AGARD-R-798 


AGARD-R-798 



ADVISORY GROUP FOR AEROSPACE RESEARCH & DEVELOPMENT 

7 RUE ANCELLE, 92200 NEUILLY-SUR-SEINE, FRANCE 


AGARD REPORT 798 

Aircraft Loads due to 


f"\ ? 'i" i # 

tJ 


JkH 09 18851 




Turbulence and their Impact on 
Design and Certification 

(Efforts avion dus a la turbulence atmospherique 
et leurs impacts sur la conception et la certification) 


This Publication was prepared at the request of the Structures 

and Materials Panel (SMP) of AGARD. Papers presented at a Workshop 

held in Lillehammer, Norway, 5 May 1994. 



NORTH ATLANTIC TREATY ORGANIZATION 

F*"fr^TRiBirnc., stat. : i i [ 

I Approved tor public release* j 
Si Distribution Unlimited i 


Published December 1994 


Distribution and Availability on Back Cover 




AGARD-R-798 



AGARD REPORT 798 


Aircraft Loads due to 
Turbulence and their Impact on 
Design and Certification 

(Efforts avion dus a la turbulence atmospherique 
et leurs impacts sur la conception et la certification) 


Accesion For 

NTIS CRA&l d 

DTIC TAB □ 

Unannounced □ 

Justification 



By.. 

Distribution/ 


Availability Codes 


Dist 


Avail and / or 
Special 


m 


This Publication was prepared at the request of the Structures 

and Materials Panel (SMP) of AGARD. Papers presented at a Workshop 

held in Lillehammer, Norway, 5 May 1994. 



North Atlantic Treaty Organization 
Organisation du traite de I’Atlantique Nord 

D'HO QuiilatiT & 


DISTRI BUTION STATEMENT A 

Approved for public release; 
Distribution Unlimited 


79950/04 070 










The Mission of AGARD 


According to its Charter, the mission of AGARD is to bring together the leading personalities of the NATO nations in the 
fields of science and technology relating to aerospace for the following purposes: 

— Recommending effective ways for the member nations to use their research and development capabilities for the 
common benefit of the NATO community; 

— Providing scientific and technical advice and assistance to the Military Committee in the field of aerospace research 
and development (with particular regard to its military application); 

— Continuously stimulating advances in the aerospace sciences relevant to strengthening the common defence posture; 

— Improving the co-operation among member nations in aerospace research and development; 

— Exchange of scientific and technical information; 

— Providing assistance to member nations for the purpose of increasing their scientific and technical potential; 

— Rendering scientific and technical assistance, as requested, to other NATO bodies and to member nations in 
connection with research and development problems in the aerospace field. 

The highest authority within AGARD is the National Delegates Board consisting of officially appointed senior 
representatives from each member nation. The mission of AGARD is carried out through the Panels which are composed of 
experts appointed by the National Delegates, the Consultant and Exchange Programme and the Aerospace Applications 
Studies Programme. The results of AGARD work are reported to the member nations and the NATO Authorities through the 
AGARD series of publications of which this is one. 

Participation in AGARD activities is by invitation only and is normally limited to citizens of the NATO nations. 


The content of this publication has been reproduced 
directly from material supplied by AGARD or the authors. 


Published December 1994 

Copyright © AGARD 1994 
All Rights Reserved 

ISBN 92-836-0002-9 


Printed by Canada Communication Group 
45 Sacre-Caeur Blvd., Hull (Quebec), Canada K1A 0S7 



Preface 


The study of the effects of atmospheric disturbances on the behaviour of aircraft started in about 1910. These disturbances 
which create discomfort for passengers and aircraft control problems for the pilots, are important dimensioning criteria in the 
field of structures. 

The first regulatory requirements in this field appeared in the 1930’s. Since this time, many other elements have intervened, 
in particular measurement of atmospheric disturbances (or their effect on aircraft), analysis of the disturbance structure, 
reduction of in-flight measurements, turbulence models, calculation of aircraft response, regulatory requirements, etc... 

The AGARD Structures and Materials Panel has always been heavily involved in this field in particular, since the mid 80’s, 
under the impetus of Dr. G. COUPRY. This renewal of interest comes from the development of transport aircraft towards 
wing configurations with high aspect ratios associated with flexible structures, the search for high-performance and optimized 
structures, the appearance of fly-by-wire controls and control laws, the development of computational power enabling more 
efficient simulation models using the results of measurements made on a large number of modem aircraft during commercial 
flights. 

In this spirit and under the impetus of J. GLASER, the Structures and Materials Panel proposed an activity which led to the 
organisation of a Workshop on the theme: “Aircraft Loads due to Turbulence and their Impact on Design and Certification”. 
This Workshop was held on 5 May 1994 as part of the 78th Meeting of the Structures and Materials Panel. 

Concentrated into a single day, this Workshop brought together the world’s leading specialists in structure calculations under 
the effect of atmospheric disturbances. It was an opportunity to link up and continue studies undertaken within other bodies 
such as the JAR Structure Study Group or the Gust Specialists’ Committee whose last meeting was organized by the FAA in 
April 1994. 

This document reproduces the papers presented. Many of them concerned the computation of the aircraft responses including 
high nonlinearities under the effect of a “continuous turbulence” type excitation. Given the recent nature of these studies, 
many complex theories and methods, which are sometimes very similar, are being developed, leading to the implementation 
of long and costly calculations. Developments for dealing with specific configurations were also presented: aircraft equipped 
with turboprops, aircraft equipped with external loads. 

The discussions held within the scope of a round table at the end of the meeting gave rise to lively and fruitful exchanges and 
additional presentations (J. HOUBOLT) which allowed us to determine the main lines for orientating the future activities of 
the Structures and Materials Panel in this field. 


Mr. D. CURBILLON 
Sub-Committee Chairman 
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Preface 


L’dtudc dcs perturbations atmospheriques sur le comportemcnt des avions a commence vers les annees 1910. Si, pour les 
passagers dcs avions de transport, elles se manifestent par unc sensation d’inconfort, et si, pour les pilotes, elles posent des 
difficulty de controle de l’avion, dans le domaine des structures, elles represented un critere de dimensionnement important. 

Les premieres exigences reglementaires en la matiere sont apparues autour des annees 30. Depuis cette epoque, de nombreux 
dldments sont intervenus, notamment la mesure des perturbations atmospheriques (ou de leurs effets sur avion), l’analysc de 
la structure des perturbations, la reduction des mesures en vol, les modeles de turbulence, le calcul de la rdponse avion, les 
exigences rdglemcntaires... 

Le Panel des Structures et Materiaux de l’AGARD a toujours eu une activite importante dans ce domaine, en particulier 
depuis 1985 sous 1’impulsion du Dr. G. COUPRY. Ce regain d’interet provient de revolution des avions de transport vers dcs 
configurations de voilure a grands allongcments associees a des structures souples, de la recherche de structures performantes 
et optimisces, de l’apparition des commandes de vo! electriques et des lois de pilotage, du developpement des moyens de 
calcul autorisant des modeles de simulation plus performants utilisant les resultats des mesures effectuees sur un grand 
nombre d’avions modenics lors de vols commerciaux. 

Dans cet esprit, et sous 1’impulsion de J. GLASER, le Panel des Structures et Materiaux a propose une activite qui a debouche 
sur 1’organisation d’un atelier sur le theme: «Les charges avion liees a la turbulence atmospherique et leur impact sur la 
conception et la certifications Cet atelier s’est tenu le 5 mai 1994 dans le cadre de la 78* mt reunion du Panel des structures et 
matdriaux. 

Concentrd sur unc joumde, cet atelier a permis de reunir les principaux specialistes mondiaux du calcul dcs structures sous 
l’effct des perturbations atmospheriques. II a etc l’occasion de faire le lien et de poursuivre les reflexions mendes dans 
d’autres organismes tels que le JAR Structure Study Group ou le Gust Specialists’ Committee dont le dernier reunion a etd 
organisd par le FAA en avril 1994. 

Ce document donne le texte des differentes communications presentees. Pour une large part, elles concement le calcul de la 
rdponse des avions comportant de fortes non-linearites sous l’effet d'une excitation du type «turbulencc continue®. Compte 
tenu du caractcrc recent de ces etudes, on assiste au developpement de theories et methodes nombreuses et complexes, parfois 
tres voisincs, aboutissant a la mise en oeuvre de calculs longs et coflteux. Ont ete egalement presentes, des ddvcloppcments 
pour traiter des configurations spdcifiques : avions equipes de turbopropulseurs, avions equipes de charges extemes. 

La discussion mende dans le cadre d’une «table rondc®, en fin de reunion, a etc l’occasion d’echanges vivants et fructueux, de 
prdsentations compldmcntaires (J. HOUBOLT) qui ont permis de determiner des axes de reflexions pour orienter les activites 
futures du Panel des Structures et Materiaux dans ce domaine. 


M. D. CURBILLON 
President du sous-comite 
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Overview of the Activities of the Ad Hoc Committee 
of Internationa] Gust Specialists 

Terence J. Barnes 

National Resource Specialist, Federal Aviation Administration 
Northwest Mountain Region Hqs 
1601 Lind Ave. (ANM - 105N), Southwest Renton 
Washington 98055 - 4056, United States 


L SUMMARY 

This paper presents the status of the work being 
undertaken by an international team of specialists to re¬ 
evaluate the gust criteria for future generations of 
commercial transports. 

The international Ad Hoc Committee of Gust Specialists 
was formed in 1986, and established the following goals: 

Criteria and means of compliance to be the same for 
FAR 25 and JAR 25. 

Reduce the total number of criteria to be met. 

Level of strength to be similar to current successful 
transports. 

Ability to handle advanced technologies such as active 
controls and load alleviation. 

Analytical work conducted by the specialists and reported 
at the annual meetings provided valuable background 
information to enable the Discrete Tuned Gust Method to 
be selected for the analysis of transport airplane dynamic 
loads due to discrete gust encounter. 

Current and future work will concentrate on improved 
ways to account for loads due to continuous turbulence, 
and handle advanced technologies such as load alleviation 
which result in non linearity. 

2. NOMENCLATURE 

AGARD - Advisory Group for Aerospace Research 
and Development 

ARAC - Aviation Rulemaking Advisory 
Committee 

CAA-UK - UK Civil Aviation Authority 
DSP - deterministic spectral procedure 

FAA - US Federal Aviation Administration 

FAR - US Federal Aviation Regulation 

JAA - European Joint Aviation Authorities 

JAR - European Joint Aviation Requirements 
NACA - US National Advisory Committee for 
Aeronautics, now NASA 
NASA - US National Aeronautics and Space 
Administration 

MFT - matched filter theory 

PSD - Power Spectral Density 


SDG - statistical discrete gust 

SSB - stochastic simulation based method 

3, BACKGROUND 

Since the recognition that turbulence produced significant 
structural loads (around 1915) there have been several 
significant steps in the development of gust criteria. The 
first gust criterion was the SHARP EDGED GUST 
formula. This was later modified to a formula specifying 
RAMP-PLATFORM GUSTS and later to ONE-MINUS- 
COSINE gusts. Finally, the criteria for CONTINUOUS 
TURBULENCE were developed. 

Additional background and history on the development of 
the current FAA criteria are given in References 1 and 2. 

In parallel with the development of criteria in the USA, 
various European Authorities developed their own criteria 
and added these as "National Variants" to the Joint 
Airworthiness Requirements of JAR 25. 

The net result was that when the gust specialist's 
committee was formed, in order to certify a transport 
airplane to both FAR 25 and JAR 25, a manufacturer was 
required to consider and possibly perform seven separate 
gust analyses, and design to the most critical loads. These 
were: 

Gust Load Formula 

PSD Design Envelope 

PSD Mission Analysis 

JAR 12.5 Chord Dynamic Gust Analysis 

CAA-UK Tuned Dynamic Gust Analysis 

RLD Negative Gust 

Round-the-clock Gust 

Further background information is to be found in 
Reference 3. 

As a result of ongoing Airworthiness Authority/ Industry 
cooperative activity and the removal of National Variants 
from JAR 25, the seven gust criteria have been reduced to 
four and modified as follows: 

Tuned Discrete Gust 
PSD Design Envelope 
PSD Mission Analysis 
Round-the-Clock Gust 


Paper presented at a Workshop on “Aircraft Loads due to Turbulence and their Impact on 
Design and Certification" held in Lillehammer, Norway on 5 May 1994. 
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Even though this significant reduction in required 
analyses was a very positive step, earlier concerns raised 
by European Industry and Airworthiness Authorities 
regarding application of the PSD criteria of FAR 25 App 
G to European transports were still unresolved. 

Furthermore, it was recognized that computer controlled 
systems such as gust alleviation would be under 
consideration on many projects, and thus the nonlinear 
effects of these systems would need to be included in any 
new gust analysis approach. 

In the absence of definitive criteria, airplanes employing 
active load alleviation concepts such as the Lockheed L- 
1011 (Reference 4) and to the Airbus A-320 (Reference 5) 
used interim design criteria and analysis procedures which 
are not discussed in this paper. 

4. DISCUSSION 

The desire to provide gust analysis loads results which 
could be used directly in a stress analysis resulted in the 
consideration of analyses in the time domain first. These 
could then provide time correlated loads directly. 

Each of the evaluations of new methods conducted and 
listed on Figure 1 will be discussed. Only key references 
are given for each method. 


New Methods Proposed to Handle Discrete Gusts 
and Continuous Turbulence 

Atmospheric 

Disturbance 

Method 

Discrete 

Gusts 

Continuous 

Turbulence 

Reference 

Statistical 
Discrete Gust 
(SDG) 

X 

X 

6 

Tuned Discrete 
Gust (TDG) 

X 


JAR 25, Ch 
13, Arndt 

91/1 

Matched Filter 
Theory (MFT) 


X 

8,9 

Stochastic 

Simulation 

Based (SSB) 
Method 


X 

11 

Deterministic 

Spectral 

Procedure 

(DSP) 


X 

12 

Spectral Gust 
Method 


X 

13 

Probability of 
Exceedance 


X 

14 


Figure 1 


Statistical Discrete Gust (SPG) : 

This method is unique in that it can handle both discrete 
gust events and continuous turbulence. The basic 
assumption of the SDG approach is that atmospheric 
turbulence can be represented by a scries of 1-cosine 
shaped ramps applied in any order. The amplitude of each 
element is adjusted to reflect equal probability of 
occurrence by scaling element amplitude according to the 
element gradient length, H, to a power: for 

continuous turbulence and for extreme amplitude 
events. 

It was first proposed by Glynn Jones in 1967, and reported 
in Reference 6. 

Although developed in parallel with the Power Spectral 
Density (PSD) method, it was not evaluated as intensely, 
especially by US manufacturers. 

Prior to proposing an in-depth evaluation of the SDG 
method by the manufacturers, the FAA requested NASA 
assistance in evaluating some of the major assumptions 
used in the SDG method. It was determined that for a 
given set of circumstances, the SDG continuous turbulence 
method produced essentially the same results as the PSD 
method, thus confirming the similarity between the H 1 ^ 
scaling law used in SDG and the von Karman spectmm 
used in PSD (Reference 7). 

An extensive evaluation by manufacturers of large and 
small transports in the US, Canada and Europe concluded 
that the method was not acceptable because of the 
laborious search procedure for the analysis of highly 
nonlinear systems, and the method was not recommended 
for consideration as a new Airworthiness Requirement. 

However, this conclusion may have been premature, since 
no other method has been found which can efficiently 
handle continuous turbulence and highly non-linear 
systems, and the SDG method is also capable of handling 
both discrete and continuous turbulence. 

Furthermore, since no single method was found to cover 
both discrete gusts and continuous turbulence, it was 
agreed that the existing situation of separate criteria would 
be continued, and each developed independently. 

Tuned Discrete Gust (TDG): 

The TDG method evolved from the method used to certify 
the Airbus A-320 (Reference 5) and includes three key 
elements. 

A gust tuning law relating gust velocity to gust 
gradient distance in the same way as for the SDG method 
using the h'^ 6 relationship. 

A reference gust velocity based on in-flight 
recordings of transport airplane loads parameters in 
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normal operations. This results in a variation of reference 
gust velocity versus altitude. 

A flight Profile Alleviation factor which accounts 
for the primary effects of the airplane operation by 
reducing the reference gust velocity at the lower altitudes 
for airplanes that predominantly cruise at high altitude. 

The TDG method is part of the JAR and has been used for 
certification of several major programs, e.g., Airbus A330 
and A340, and is part of the certification basis of the 
Boeing 777 and Cessna Citation X, among others. The 
FAA evaluates each new program individually, 
considering the equivalency of the TDG method to 
existing regulations. 

Matched Filter Theory (MFT) ; 

The MFT method identifies the gust velocity versus time 
profile which produces the maximum unit value of a 
specific load quantity. 

Unit impulses of various strengths are applied through a 
von Karman filter to the known dynamics of an airplane, 
and the impulse responses of a given load quantity are 
ascertained. For that given load quantity, the load impulse 
response waveforms are reversed in time and played back 
through the gust filter and the airplane system to 
determine the critical load quantity in question and the 
associated correlated loads. The matched excitation or 
output waveforms are a function of the initial impulse 
strengths. The whole process is repeated for impulses of 
various strengths to determine the most critical load 
response together with the associated correlated loads. 

With a nonlinear system, super position and other 
convenient features of linear MFT cannot be used, and a 
search procedure is the only practical means of finding a 
solution! 

NASA reported results for time correlated loads for linear 
systems in Reference 8, and for nonlinear systems in 
Reference 9. 

Stochastic Simulation Based fSSB) Method : 

The SSB method uses a random white noise input signal 
fed through a von Karman filter to produce a gust velocity 
time history with a von Karman power spectrum which is 
applied to an airplane with linear or nonlinear systems. 

Several simulations are performed. First, a search of time 
histories for a particular structural element load locates 
points in time where a peak load is within a specified 
range. The narrow response range is chosen for nonlinear 
systems to isolate the nonlinear nature of the airplane 
responses. Second, a sufficient amount of time prior to 
and following the peak load is selected, and the 
corresponding time history segments of the input 
excitation, the gust velocity, and the response waveforms 
are extracted. 


The average waveforms of each type from this process 
converge to provide a gust velocity description that 
maximizes a given load quantity, hence that load quantity 
and the associated correlated loads are obtained. 

Under a small contract with the FAA, NASA developed a 
SSB method to the same level as the MFT method and 
developed loads on a representative model with nonlinear 
systems, and compared the results and resources needed. 

Reference 11 summarizes the SSB method and presents 
numerical results, and concludes that the approach 
requires excessive computer resources. 

At the conclusion of this study the FAA accepted NASA's 
recommendations and released the computer code for 
MFT for industry evaluation (Reference 10). 

Deterministic Spectral Procedure (DSP) : 

The Deterministic Power Spectral Procedure, DSP, was 
introduced by Jones and has been further developed by 
Noback (Reference 12). It can be regarded as the 
translation of the continuous turbulence power spectral 
density, PSD, method from the frequency domain into the 
space or time domain. The critical gust profile, giving 
maximum response for any given load quantity and 
correlated associated loads, is obtained directly. Noback 
has compared the method with Matched Filter Theory, 
and with direct simulation. It gives identical results to the 
power spectral approach for linear systems. Results have 
been obtained for a nonlinear system using equivalent 
gains to obtain the critical gust profile, and these results, 
also, have been compared with direct simulation. 

Spectral Gust Method : 

This is based on discrete gusts which have the correct 
power spectral density. Thus, the spectral gust method 
represents a direct translation of the PSD method into the 
time domain. 

The method requires no optimization. 

However, there can be several Spectral Gusts resulting in 
different design loads for nonlinear systems, depending on 
the definition of design load. 

One solution might be to develop a few Spectral Gusts for 
the certification of each project, however, this would be a 
cumbersome process (see Reference 13). 

Probability of Exceedance : 

This method, which was proposed several years ago by 
Noback (Reference 14), determines limit loads for a 
nonlinear aircraft using the design envelope criterion. 

The method requires a time domain stochastic simulation 
of the pertinent responses to be evaluated for a number of 
input RMS gust intensities. The use of three gust 
intensities is suggested. 
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The criteria are defined so that the probability of 
exceeding the design load with the nonlinear system 
would be the same as the probability of exceeding the 
design load with a linear system. The result must be 
scaled to match linear traditional PSD results for a linear 
system. 

5. FUTURE PLANS 

Gust analysis methodology and design criteria are being 
worked by two overlapping committees. The Gust 
Specialists Committee defines, evaluates and makes 
recommendations to the ARAC FAA/JAA Loads 
Harmonization Working Group, which continues the 
process into possible rulemaking activity. 

The tuned discrete gust analysis method, which handles 
isolated gusts, has followed this procedure, and was 
incorporated in JAR 25 at Amendment 91/1 effective 12 
April 1991, and an FAA Draft NPRM has been prepared. 

Activity to re-evaluate the Power Spectral Density gust 
analysis method, which handles continuous turbulence, 
and consider other methods such as Matched Filter Theory 
which can efficiently handle nonlinear systems, is in work 
in both committees. 

The status of the various gust evaluations is given on 
Figure 2. 

The following needs to be accomplished: 

Complete re-evaluation of current PSD criteria of 
FAR 25, App G - U s jg ma and mission analysis. 

Complete evaluation of MFT one dimensional search 
method. 

Recommend method for analysis of airplanes with 
nonlinear systems flying through continuous 
turbulence. 

At the analysis level, consistent structural and 
aerodynamic modeling needs to be addressed. 

In Reference 15, John Glaser of deHavilland Inc., Canada, 
proposed an AGARD study with the following title: 
"State-or-the-Art Methods in the Computation of Aircraft 
Gust Loads". 

The proposed activity is described as follows: 

The activities to date on the Flight of Flexible Aircraft in 
Turbulence have concentrated on the models of 
atmospheric turbulence and the determination of their 
properties from flight data recordings, the most familiar 
being the discrete, (1-cos) time domain gust model and the 
von Karman frequency domain or PSD gust model. The 
importance of the gust model is that it forms the 


"excitation" side of the equations of motion and the gust 
intensity value used is essential to ensuring an adequate 
level of design strength and therefore safety. 

There are, however, other aspects in calculating gust loads 
that must be considered because of their potentially large 
influence on the design levels obtained. These have to do 
with aircraft modeling methods and the analytical 
procedures used. To date, little attention has been paid to 
these aspects within AGARD, or elsewhere for that 
matter. The following areas are offered as examples: 

Flight mechanic modeling methods and the methods 
used to achieve appropriate flight stability 
characteristics. 

Structural dynamic modeling methods. 

The methods used to combine the flight mechanic and 
the structural dynamic equations of motion. 

Computational techniques for solving the equations of 
motion. 

The calculation of response loads from the solution 
for gust response. 

A large variety of assumptions and approximations are 
made in the calculation of gust response loads that can 
strongly influence the loads results obtained. There is a 
clear need to identify where rigorous modeling methods 
and analytical procedures are required, to establish what 
those methods should be and to make those methods and 
procedures known so they can be applied uniformly 
throughout the industry. 

But in recent years, the calculation of gust loads has 
become even more difficult with the advent of active flight 
control systems with their highly nonlinear characteristics. 
This has stimulated the investigation of a variety of new 
gust modeling methods, equivalent to the von Karman 
model, labeled SDG, MFT. DSP, DPSD, etc. These new 
methods are viewed by some gust loads analysts as being 
too complex to understand, to apply in routine calculations 
or to be accepted by airworthiness authorities. 

There is a clear need to understand what effect nonlinear 
control systems (and other sources of nonlinearity) can 
have on gust design loads. If significant, it is important to 
have some simple but effective simulations of the von 
Karman model of continuous turbulence that can be 
accepted by industry and airworthiness authorities alike. 

The purpose of this activity is to organize a workshop of 
gust loads specialists from industry, research 
establishments and airw orthiness authorities to review the 
state-of the-art of gust loads computations with the view to 
identifying aircraft modeling methods and analysis 
procedures that can significantly influence design but are : 



1-5 


Inconsistently or incorrectly applied. 

Inadequate to give reliable results. 

Impractical or do not exist at all for current or 
emerging designs, in order to give framework and 
direction for a working group on gust loads methods. 

The ultimate objective of this workshop is to establish a 
terms of reference for a working group that will be 
charged with establishing recommended design gust load 
methods, that will define those methods with 
demonstrations of their accuracy and effectiveness and 
that will issue the results of their work in a manual or 
AGARDOGRAPH. 

The AGARD Working Group can work essentially 
independently of the other committees, and perform a 
useful function for the aviation community. 

6. RECOMMENDATIONS 

Incorporate TDG in FAR 25. 

Continue evaluation of methods to develop design 
loads on airplanes with highly nonlinear control systems. 

Possibly reconsider SDG if no other method can 
handle the problem of nonlinear control systems directly 
and efficiently. 

Continue review of the current PSD method of 
FAR 25, Appendix G. 

Define a limited number of patches of 
moderate/severe turbulence to be used to test the operation 
of airplane systems in the nonlinear range. These patches 
of turbulence could also be studied to determine the 
validity of the current assumptions, e.g., the 1-cosine 
shape, and to determine if it would be possible to obtain an 
improved and much simplified time history 
characterization of design gusts that would work equally 
well for nonlinear as well as linear systems. 

Proceed with the formation of an AGARD 
Working Group to review and recommend analytical 
modeling methods and procedures for dynamic gust 
analysis. 

Continue to support the efforts to collect 
worldwide operational gust statistics from the current 
commercial fleet, and collect data from high altitude 
research aircraft to support the development of HSCT gust 
design criteria. 
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SUMMARY 

During the past ten years the extensive use of electronic flight 
control systems has introduced a high level of non-linearity, 
to the behaviour of civil transport aircraft, in response to both 
pilot inputs and external influences such as turbulence. These 
systems that control, protect and, in some instances, alleviate 
the loads experienced by the aircraft have had a significant 
impact on the prediction of aircraft loads. There is an in¬ 
creasing need both to model control system non-linearity, to 
avoid designing control systems that degrade structural per¬ 
formance, and to demonstrate the effectiveness of alleviation 
systems for aircraft certification. 

The introduction of non-linear flight control systems presents 
little problem when using time-based simulations to predict 
aircraft loads due to discrete gusts. Frequency domain, 
Power-Spectral Density (PSD), analysis used to predict 
aircraft loads due to Continuous Turbulence (CT), however, 
is severely restrictive, requiring a linear mathematical model 
of both the aircraft and control systems. Over the past ten-to- 
fifteen years various analysis techniques have been devel¬ 
oped that provide time-based interpretations of the CT 
atmosphere and, therefore, means by which loads due to CT 
can be predicted for aircraft with non-linear flight control 
systems. Recent involvement in the design of Airbus aircraft 
has led to British Aerospace Airbus Limited taking an active 
role in both developing and using these non-linear analysis 
techniques. 

By drawing on experience gained during the design of recent 
Airbus aircraft, the impact of non-linear flight control systems 
on the prediction of aircraft loads due to turbulence is dis¬ 
cussed. 

LIST OF SYMBOLS 

A, inertia coefficient matrices. 

B, B^ aerodynamic damping coefficient matrices. 

C, C^ aerodynamic stiffness coefficient matrices. 

8 vector of control deflections. 

D structural damping coefficient matrix. 

E structural stiffness coefficient matrix. 

■E[x|y=y d ] expected value of load x given that load y is at 
its design value, y d . 

$>((B) normalised von Karman power spectrum. 

r| d turbulence "probability" factor { =TJ./a w }. 

H(j®) aircraft "closed-loop" frequency-response- 

function {=Y(j(o)/u g (jco)}. 

H y 0'co) frequency-response-function for load y, a 

member ofH(jco) { =y(j<o)/u g (jco) }. 

I identity matrix. 


K eq equivalent gain. 

L turbulence scale length. 

N[y d ] rate of exceedence of design load, y d . 

N 0 rate of positive going, zero crossings, for a 

linear model. 

p(o) probability density function of u. 

q modal coordinates (states), 

p air density. 

p correlation coefficient between loads x and y. 

R c , R ca aerodynamic forcing vectors due to control 
surface deployment. 

R g , R ga aerodynamic forcing vectors due to turbulence. 
ct w standard deviation (root-mean-square value) of 

turbulence, turbulence intensity. 
a y standard deviation (root-mean-square value) of 

load y. 

T turbulence time constant { =L/V }. 

U c Design Envelope Analysis peak turbulence 

intensity. 

u ce unit of computational effort. 

u g turbulence input. 

V true airspeed. 

M angular frequency. 

Y "Interesting Quantity" vector containing load 
quantities, lift coefficients, accelerations, etc. 

y a load quantity, a member of Y. 

y d the design value of load y. 

1.0 INTRODUCTION 

The task of predicting aircraft flight loads due to turbulence 
has grown considerably over the past ten years, largely due to 
the extensive use of electronic flight control systems on 
modem civil transport aircraft. 

Electronic flight control laws, with altitude and air¬ 
speed dependent gains, have meant that the use of 
extrapolation to focus on design critical flight regimes 
and aircraft configurations is not certain to identify the 
most critical loading conditions. Consequently, Design 
Envelope Analysis (DEA) requires a far more extensive 
coverage of the flight envelope and aircraft configur¬ 
ations than is required for an aircraft without electronic 
flight control systems. Non-linear flight control laws 
accentuate this problem. 

The use of load alleviation systems has served to 
eliminate the concept of a "turbulence designed" air¬ 
craft or a "manoeuvre designed" aircraft by ensuring 
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that envelope loads due to manoeuvres and envelope 
loads due to turbulence are of similar magnitude. Con¬ 
sequently, loads due to turbulence and manoeuvres must 
be analysed in equal depth throughout the design cycle. 

The above observations serve to demonstrate both the im¬ 
portance and the increasing cost associated with the pre¬ 
diction of aircraft loads due to turbulence, for both the design 
and certification of an aircraft with non-linear flight control 
systems. 

It has been the policy of Airbus Industrie to use DEA both to 
design and certificate its "fly-by-wire" aircraft. Consequently, 
all the following discussion will centre on the prediction of 
aircraft loads to establish a design envelope For the purposes 
of this discussion the aircraft model is assumed linear and, 
when combined with a linear control system, forms the 
datum against which the impact of non-linear control systems 
are judged. However, much of the discussion is equally app¬ 
licable for situations w'herc the non-linearity is exhibited by 
the aircraft model. 

2.0 THE DATUM MODEL OF THE AIRCRAFT 

The manufacturers aim is to minimise the weight of the air¬ 
craft and thereby improve its performance and marketability, 
whilst airworthiness requirements exist to provide rules that 
are not necessarily rational, but that aim to ensure a safe 
structure based on a large database of past and present 
aircraft. The loads engineer must therefore strive to satisfy 
both requirements and accordingly choose modelling tech¬ 
niques of suitable complexity, accuracy and cost 

The model of the aircraft has developed over many years to a 
standard where, for the purposes of loads prediction, it 
adequately represents the aircraft. In general, flight tests have 
confirmed this although modification to the data, feeding into 
the model, is usually required to obtain a satisfactory match 
between the measured and predicted response. Prior to 
matching it is possible that some predicted aircraft loads will 
differ from measured aircraft loads by up to 20%, however 
the predicted trends in aircraft loads are usually accurate. 

The aircraft model is not intended to "model the truth", rather 
to provide a reasonable representation of the aircraft that, in 
conjunction with other design parameters (e.g. factors of 1.5 
between limit load and ultimate load), will ensure the 
structural integrity of the aircraft without yielding over¬ 
weight designs. This does not imply that improvements to 
modelling techniques should not be sought, rather that all 
new methods must consider the accuracy and completeness of 
all model components; the purpose for w'hich the model will 
be used and the external parameters that may influence the 
level of accuracy required. 

2.1 The Mathematical Model of the Aircraft Model 

Over many years the development of the mathematical model 
of the aircraft has been guided by: 

• Its performance against flight test. 

• The computing pow'er available. 

• The airworthiness requirements. 

• The manufacturers desire to minimise weight. 

• Design costs. 

The model that has evolved combines both the structural and 
aerodynamic properties of the aircraft into a second order, 
linear, ordinary differential equation (1) described in terms of 
modal coordinates, i.e. all shapes that the aircraft can assume 
are represented by the sum of its natural vibration (mode) 


shapes, in various proportions. By eliminating the mode 
shapes that correspond to the higher frequencies of vibration 
it is possible to keep the aircraft equations of motion to a 
manageable size, with only minor degradation to accuracy. 

Aq+(pVB + D)q + (pV 2 C + E)q = pVR g u g +pV 2 R c 5 (1) 

The aircraft "Interesting Quantity" (IQ) set comprises forces, 
moments, torques, accelerations, etc. at various stations with¬ 
in the aircraft structure and are recovered through an aux¬ 
iliary equation (2) that sums the external forces acting on the 
structure. The forces and moments could also be recovered as 
the sum of internal forces, however this can prove unreliable 
when w orking with a limited number of mode shapes. 

X = pVR ga u g + pV 2 R ca 5 - A a q - pvB a q - pV 2 C a q (2) 

2.2 The Introduction of a Flight Control System 

The combination of the aircraft and control system is 
simulated by "wrapping" the control system equations around 
the aircraft dynamics represented by equations (1) and (2), 
see figure 1. 



Figure 1 . Aircraft simulation 


In the time-domain the aircraft equations are rearranged into 
state-space form and solved using a suitable integration 
method for first order, ordinary differential equations. The 
control system is modelled as a digital circuit driven by the 
feedback quantities (e.g. normal accelerations, pitch rate, 
etc.) that form a subset of the aircraft IQ set. At specified 
time-intervals within the simulation the IQ set is calculated; 
the control system activated and the commanded control de¬ 
flections updated. 

In the frequency-domain, the aircraft dynamics are described 
by Frequency-Response-Functions that relate the "Interesting 
Quantities" to the control and turbulence inputs, H. and II, 
respectively. Similarly, a linear control system is represented 
by a Frequency-Response-Function (FRF), G. 'Hie "closed- 
loop" FRF, relating the IQ set to the turbulence input, is 
calculated using equation (3). 

H( jm) = H c ( jro )£( jeo )G( jeo )H g (jeo) + H g (jeo) (3) 

where r(je») = [l-G(jeo)H c (jeo)] 

The embellishment, denotes that the FRF contains only the 
feedback quantities. 

The frequency-responses of the IQ set, hence the aircraft 
loads, are obtained via equation (4). 

Y(jw)= H(jw)u g (jw) (4) 

3.0 DESIGN ENVELOPE ANALYSIS 

The current airworthiness requirements [1][2] request that 
two types of analyses are performed to predict aircraft flight 
loads due to turbulence. These comprise the prediction of 
loads arising from both discrete gust events and CT. In turn, 
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there are two alternative approaches to determining aircraft 
limit loads due to CT: Mission Analysis and DEA. It has been 
the policy of Airbus Industrie to use DEA to design and 
certificate its "fly-by-wire" aircraft, so this forms the basis for 
further discussion. 

The concept of the Statistical Discrete Gust (SDG), intro¬ 
duced by Jones [3], attempted to combine the two analyses 
into a single analysis. However, the SDG has not been 
adopted as an airworthiness requirement and therefore offers 
little to the manufacturer, except maybe as a starting point for 
developing methods suitable for the prediction of loads for 
nonlinear aircraft models. 

The Tuned Discrete Gust (TDG), comprising a series of "1- 
cosine" discrete gusts of varying gradient length, is used to 
predict loads due to discrete gust events and is best suited to 
a time-domain analysis. CT, however, is represented by a 
stationary, Gaussian, stochastic process for which the turb¬ 
ulence is described by the isotropic power-spectral density 
function attributed to von Karman, equation (5). Consequent¬ 
ly, for a linear aircraft model, CT analysis is best performed 
in the frequency-domain. 


<£>(«) = ct>( W ) = ct’T 


1 + %(1.339Tco) 2 
{l + (1.339Tco) 2 }'^ 


(5) 


3.1 Tuned-Discrete Gusts 

For a given design condition (e.g. aircraft configuration and 
flight point) the TDG requires the prediction of aircraft loads, 
due to a family of "1 -cosine" discrete gusts of varying 
gradient length and amplitude, using a time-based analysis. 
DEA requires determination of the envelope (e.g. maximum 
and minimum) loads of all the design points considered. For 
design purposes, however, the maximum and minimum str¬ 
esses, that are a function of the forces and moments acting at 
a given point in time, are of more importance. There are two 
ways that these stresses may be recovered. 

1. Directly, by including stress quantities in the IQ set. 

2. Indirectly, as a function of time-correlated loads. 

The first method eliminates the need for time-correlated 
loads, but will require a complete turbulence analysis each 
time the structure is changed. The second method, albeit an 
approximation to the maximum and minimum stresses, en¬ 
ables structural designers to make modifications to local 
structure without recourse to further turbulence analyses. 

For logistical and financial reasons the second method is 
preferred, so it is necessary not only to recover the envelope 
loads but, also, the loads over the rest of the aircraft that 
coincide (correlate) with the envelope values. Correlated load 
sets are easily obtained for time-based analyses by sorting 
through the load response time-histories. 

The introduction of a non-linear flight control system has 
little impact on the discrete gust analysis techniques, only on 
the loads predicted. 

Note 1: computational effort. 

The calculation of a correlated load set, resulting from a 
single "1-cosine" discrete gust input, will be used as the base 
unit for a measure of computational effort. Typically, one 
unit of computational effort (lu CE ) represents the work re¬ 
quired to establish a correlated set of, say, Jive-hundred 
"Interesting Quantities" from three seconds (approximately 


200 time-points) of time-history response per IQ. For a given 
aircraft configuration and flight point, a TDG analysis com¬ 
prising a family of, say, sixteen "1-cosine" gusts will require 
16u ce (sixteen units of computational effort). 

3.2 Continuous Turbulence, Random Process Theory 

For a given design condition and a linear aircraft model, the 
prediction of aircraft loads due to CT uses a frequency- 
domain analysis, often referred to as Random Process Theory 
(RPT). The aircraft is subjected to turbulence that conforms 
to the von Karman power spectrum (5) and a statistical des¬ 
ign load is determined through a "probability" factor, r| d , that 
relates the design load, y d , to the root-mean-square value of 
tire load, a y [4], 

y d = = A y Tl d CT w = A y U a (6) 


where 


A y=7 L = J|j 0 1 H y(J co f <I)( ®) d “ (7) 

The CT problem was originally assumed to be linear, so the 
"probability" factor, q d , and the turbulence intensity, <r w , 
were combined into the single value of U a defined in the 
airworthiness requirements [1][2] as a function of altitude 
and airspeed. 

Since the design load is determined on a statistical basis, the 
correlated loads can be determined on the basis of "equal 
probability", Noback [5], The correlated load set represents 
the expected values of the other loads when the envelope load 
equals its design value, y d . For a linear aircraft model the 
response of all loads to Gaussian turbulence will also be 
Gaussian, so the loads, y j5 that correlate to the design load, 
y.=y d , are given by the conditional expectation: 

'E[yjly l = y d ] = y d Py i y i - li - ( g ) 

°yi 

which, for the linear aircraft model, yields 

‘E[y j |y i = y d ] = f H y . ( jto )H; (jeo )*D(jco )dco (9) 

Ay, Jo 

A single CT analysis, for a given aircraft configuration at a 
given flight point, will involve the calculation of tile aircraft 
"closed-loop" FRF (3) at, say, 400 frequencies. The frequency 
-responses of the "Interesting Quantities" are obtained using 
equation (4) and the correlated set formulated by integration 
using equations (6) through to (9). This requires about 10u, e 
(ten units of computational effort). 

Note 2: for a linear system, Matched Filter Theory (MFT) [6] 
offers an equivalent time-based approach to that of RPT. For 
a linear aircraft, however, RPT requires less computational 
effort and is, therefore, the preferred approach. 

The unsuitability of RPT for non-linear problems means that 
alternative methods must be sought to deal with non-linear 
flight control systems. Consequently, the prediction of aircraft 
loads due to CT is where the introduction of non-linear flight 
control systems has had most impact on the methods used. 
Before discussing this further, however, the following two 
items will be addressed: 
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• The definition of the turbulence intensity for sto¬ 
chastic analysis using the "principle of equivalent 
safety". 

• Non-linear flight control systems and the types of 
non-linear problem encountered. 


cy of exceedence variation, with turbulence intensity, shown 
in figure 2. Equation (13) would then be used to generate the 
rate of exceedence curves necessary to establish the non¬ 
linear design load. Clearly this is an impractical approach, 
therefore alternatives must be sought. 


4.0 "PRINCIPLE OF EQUIVALENT SAFETY" 

It was noted in section 3.2 that for CT DEA, due to the 
assumption of linearity, the value of the turbulence intensity, 
tr w , was combined with a "probability" factor, q d , to yield a 
single value, U B , that is defined in the airworthiness require¬ 
ments [1 ][2] as a function of aircraft speed and altitude. 

U a =T ld C w (1°) 

The choice of U a for DEA is founded on the basis of ach¬ 
ieving a constant frequency of exceedence of load at all alt¬ 
itudes [4] and by choosing a frequency of exceedence of 
design load that experience has shown to be safe. Stochastic 
based, non-linear, analyses that utilise the "principle of 
equivalent safety" (i.e. maintain a statistical level of safety 
consistent with the linear definition of design load) require 
knowledge of the turbulence intensity, a w . Turbulence intens¬ 
ities arc not defined for DEA, however they are for mission 
analysis [4]: turbulence is defined as a quasi-stationary, Gaus¬ 
sian, process with the same von Karman power spectmm, but 
with a slowly varying root-mean-square intensity with a prob¬ 
ability density function (11). 



The parameters: p,, p 2 , b ( and b 2 are specified in the air¬ 
worthiness requirements [1][2], 

The normalised rate of positive going exceedences (12) of the 
design load, y d , for a given turbulence intensity is derived 
from a formula attributed to Rice [7], 


N[y d ,tr w ] 

N 0 




( 12 ) 



Figure 2. Contribution of turbulence 
intensity to design load. 


Vinnicombe et al. [9] found that, for their particular non¬ 
linear problem, it was adequate to define a single value for 
the "probability" factor, ri max (e.g. fi max =3), using equation 
(14) for a load critical flight point. Equation (10) was then 
used to define turbulence intensities at other flight points. 

An alternative, suggested by Noback [8], is to perform a non¬ 
linear analysis for several turbulence intensities and use a 
weighted average (based on the curve in figure 2) of the 
resulting loads to determine the non-linear design load. 

The choice of turbulence intensity is not always straight 
forward and will depend very much on the aircraft and the 
nature of the non-linear problem. If the chosen turbulence 
intensity does not exercise the non-linearity then design loads 
will be biased tow'ards those for a linear model, the converse 
is also true. Clearly, this introduces an ambiguity into non¬ 
linear stochastic methods that is best addressed on a problem 
by problem basis, although there is a need for guidance 
within a future CT airworthiness requirement. 


More importantly, the frequency of exceedence of the design 
load is obtained by integrating the product of equations (11) 
and (12) over the range of turbulence intensities. 


N[y d ] 

N 0 


’ Nlyj.g, 

N n 


-p(CT w )dc7 v 


(13) 


Noback [8] and Vinnicombe et al. [9] have demonstrated that 
exceedences of the design load are dominated by the "storm 
turbulence", represented by p 2 and b 2 in equation (II), and 
that the values of turbulence intensity contributing to exceed¬ 
ence of the design load are concentrated in a narrow' band, as 
seen in figure 2. 


The value of r| d that corresponds to the turbulence intensity 
contributing most to the exceedence of design load for a given 
flight point is given by equation (14). 



A true, stochastic, non-linear analysis, based on the "principle 
of equivalent safety", would require many simulations using 
differing turbulence intensities, guided by the rate of frequen- 


5.0 NON-LINEAR CONTROL SYSTEMS 



Figure 3. Non-linear control function. 


Control system non-linearity manifests itself in an assortment 
of ways. In the context of turbulence analysis, the major non¬ 
linear problems encountered within the Airbus electronic 
flight control systems may be categorised as follows, with 
reference to figure 3. 

"Unidirectional switching" 

This type of non-linearity involves switching from, say, 
a normal flight control function to a flight protection 
function (e g. stall protection). Switching to the prot¬ 
ection function is triggered by the turbulence and re¬ 
version to the normal flight control function is only poss¬ 
ible with pilot intervention. Consequently, a turbulence 
analysis, that does not permit pilot intervention, will 
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start in the normal control function, but may end within 
the protection function. For discrete turbulence events 
the behaviour of the system is typified by equations (15) 
and (16). 


u g (t = 0) = 0 and y(t = 0)=0 j 
u g (t—>oo) = 0 and y(t-»oo)*0j 

/(u g W(-u g ) (16) 


"Additive laws" 

These functions involve the operation of additional con¬ 
trol laws in combination with the normal flight control 
laws and are typified by load alleviation systems. The 
alleviation function is activated by the turbulence and 
will cease operation when the turbulence subsides. The 
behaviour of the system, for discrete turbulence events, 
is governed by equations (16) and (17). 

u g (t = 0) = 0 and y(t = 0)=0 1 
u g (t—>oo) = 0 and y(t->oo) = 0j 

"Elemental Non-linearity" 

This is the existence of non-linear control law elements, 
such as limits and rate limits, within a control law func¬ 
tion. In general the behaviour of this non-linearity is 
governed by equations (16) and (17). 

The type of non-linear problem will influence the method 
used to predict aircraft loads. 

6.0 METHODS FOR NON-LINEAR CT ANALYSIS 

A multitude of non-linear methods have been developed that 
enable the prediction of loads. All are based on alternative 
realisations of the spectral properties of the atmosphere 
modelled by equation (5) and are categorised as: 

• RPT based stochastic methods. 

• Time-domain stochastic methods. 

• Deterministic "spectral" gust methods. 

• Deterministic "worst-case" gust methods. 

Some of the more well-known methods in each category are 
summarised in figure 4, the computational effort required to 
predict design loads, for a given aircraft configuration at a 
given flight point, is indicated by the left hand scale. The four 
categories of method are discussed in more detail in the 
subsequent sections. 

6.1 RPT Based Stochastic Methods 
These methods require linearisation of the non-linear prob¬ 
lem and the use of RPT as outlined in section 3.2. The linear¬ 
isation can take two forms: 


login (“J 



RPT Linearisation and Random Process Theory 
EG Linearisation with Equivalent Gain, Noback [10] 

SDG Statistical Discrete Gust, Jones [3] 

SGM Spectral Gust Method, Lusebrink et al. [11] 

TSS Time-plane Stochastic Simulation 

TSS* Modified TSS, Vinnicombe et aL [9] 

DSG Design Stochastic Gust, see note 3 
MFT Matched Filter Theory, Prototzky et al. [6] 

MFT 1D MFT 1-dimensional search, Prototzky et al. [12] 
MFT md MFT multi-dimensional search, Prototzky et al. [ 12] 
DSP Discrete Spectral Procedure, Jones [13] 

Figure 4. Methods for non-linear CT analysis. 


• Replacing a non-linear element by an Equivalent Gain 
(EG): In cases of "elemental non-linearity", where the 
influence of the non-linearity on the aircraft response is 
significant, it is possible to approximate the non-linear 
element by an equivalent gain (see figure 5), provided 
that the non-linearity is an anti-symmetric function. 


K„ 


°v( V Jo 


fv/(v)p(v)dv (19) 

Jo 


provided that /(v) = -/(-v) (20) 



Figure 5. "Elemental non-linearity" represented 
by an equivalent gain. 


Direct removal the non-linear function: currently this is 
the only means available for dealing with "unidirectional 
switching" non-linearities and involves replacing one CT 
analysis by two, i.e. a CT analysis is performed inde¬ 
pendently for both the normal flight control function and 
the protection function. 

Direct removal of "elemental non-linearity" is only 
possible if, in the context of a turbulence analysis, the 
non-linearity has no influence on the aircraft response. 
This is often the case for limits that restrict control 
surface movement, i.e. commanded control deflections, 
arising from turbulence, are usually well within such 
limits. 


Again the use of RPT is enabled, however, the root- 
mean-square value, ct v , of the input into the non-linear 
element is a function of the aircraft "closed-loop" FRF, 
equation (3) which, in turn, is a function of the equiva¬ 
lent gain, hence equation (19) must be solved iteratively. 
Consequently, the computational effort is greater than 
that for the linear RPT analysis. If there is more than one 
significant non-linearity, then both the complexity of the 
method and the computational effort will increase. 
Equation (19) also implies knowledge of the turbulence 
intensity and the reader is referred to the discussion in 
section 4.0. 
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The described frequency-domain techniques require amongst 
the lowest amount of computational effort of all the methods 
for CT analysis, but are only suitable for a very' limited set of 
non-linear problems. 

6.2 Time-based Stochastic Methods 

Time-based stochastic methods require the generation of a 
random time-history, or time-histories, that conform to the 
power-spectral density of the atmosphere described by 
equation (5). The Time-based Stochastic Simulation (TSS) 
requires generation of sufficient length of random turbulence 
time-history so that the spectral properties are adequately 
represented and the turbulence probability distribution func¬ 
tion is Gaussian to at least three standard deviations. This 
turbulence history is used as an input to a time-simulation 
using the techniques outlined in section 2.2. Non-linear 
design loads may then be computed either directly, using 
auto-correlation functions, or by using tire "principle of 
equivalent safety", based on frequency of exceedence of 
design load, Noback [14]. 

In an attempt to reduce computational effort, Vinnicombe et 
al. [9] used shorter simulation times and introduced the con¬ 
cept of a linear reference model to compensate for any 
statistical distortion present in the ensuing analysis. The non¬ 
linear design load is defined as that which has the same 
normalised rate of exceedence as the design load, for the 
linear reference model, given by equation (6). 

TSS methods, that use the "principle of equivalent safety" to 
determine tire design load, offer no direct means of obtaining 
correlated loads, so approximate methods must be sought, i.e. 
an averaging technique may be applied. 

The computational effort required for TSS methods borders 
on the limits of practicality within existing design parameters 
and, in this respect, the alternative offered by Vinnicombe et 
al. is more acceptable. 

6.3 Deterministic "Spectral" Gust Methods 

Deterministic, time-based, methods differ from the stochastic 
time-based methods because the spectral properties of the 
turbulence arc represented by a discrete gust In the case of 
the "spectral" gust no attempt is made to tune the gust shape 
to the dynamics of the aircraft and the method does not 
adhere to the "principle of equivalent safety". These methods 
offer a deterministic representation of equation (5), however, 
the "spectral" gust is not unique. In general, these methods 
are computationally practical within the design environment. 

6.4 Deterministic "Worst-Case" Methods 

Deterministic "Worst-case" methods derive from the "match¬ 
ed filter principle" introduced by Papoulis [15] and involve 
the evaluation of deterministic inputs that satisfy a prescribed 
energy constraint and that maximise the response of a given 
load. Like the "spectral" gust methods, the "worst-case" gust 
offers a deterministic representation of equation (5) but, in 
this case, the response is tuned to the dynamic characteristics 
of the aircraft. "Worst-case" methods require a constrained 
optimisation procedure (for non-linear problems) to deter¬ 
mine the "worst-case" gust: each gust is tuned to the load 
quantity in question. Consequently, these techniques require 
considerable amounts of computational effort, 10,000 to 
100,000 times that required for a RPT analysis, depending on 
the complexity and completeness of the search. 

Note 3: The Design Stochastic Gust (DSG), as shown in fig¬ 
ure 4.0, is an approximate "worst-case" method (based on the 


"worst-case deterministic input" developed by Jones [16]) for 
a non-linear aircraft. The "critical gust input" for a linear¬ 
ised model is determined and then used as an input to the non 
-linear model. 

7.0 NON-LINEAR CT DEA 

It was noted in the introduction that the electronic flight 
control laws, in particular non-linear flight control laws, have 
served to increase, maybe double, the computational effort 
required to predict aircraft loads for both the design and 
certification of modem civil transport aircraft. It is essential 
that the methods employed to deal with the prediction of 
aircraft loads for the non-linear aircraft ensure the structural 
integrity of the aircraft without requiring unrealistic levels of 
computational effort. 

Note 4: A DEA for an aircraft with an electronic flight 
control system may require consideration of several hundred 
design conditions (e.g. flight point, mass configuration and 
aerodynamic configuration) per design loop and require that 
both CT and TDG analyses are performed. 

There is no single prediction method for CT analysis that 
satisfies design needs for al! types of non-linear problem, i.e. 
there are no methods for dealing with "unidirectional switch¬ 
ing" (except that described in section 6.1) unless an assumpt¬ 
ion is made regarding pilot intervention This problem is not 
unique to CT analysis, but presents less of a problem when 
performing a discrete gust analysis. 

"Worst-case" methods are impractical in the design environ¬ 
ment unless the size of the IQ set is severely restricted and 
the coverage of the flight envelope significantly reduced, thus 
jeopardising structural integrity. Alternatively design time- 
scales and costs will be increased, thus jeopardising the man¬ 
ufacturer's ability to produce cost-effective designs. Given the 
"state-of-the-art" in computing, these methods are too 
complex to be justified: the model of the atmosphere is at 
best approximate and the accuracy of the mathematical model 
(refer to section 2.0) used to represent the aircraft docs not 
justify the complexity and computational effort required for a 
"worst-case" approach. Added to this there is no guarantee 
that the "worst-case" gust found by these methods is, in fact, 
the "worst-case" and the probability of encountering such a 
gust in flight is likely to be extremely remote. 

The more practical non-linear prediction methods will require 
less than 1000 units of computational effort. Of these, the 
equivalent gain method provides one of the simplest 
approaches, but is very limited its application: the non¬ 
linearity must obey equation (20). This method has been 
employed in the design of recent Airbus aircraft to account 
for the effect of a limit within the autopilot control system. 
The equivalent gain method, however, is not suitable for the 
"additive law" associated with gust load alleviation systems. 

The effectiveness of the A320 Gust Load Alleviation Funct¬ 
ion (GLAF) was demonstrated using the modification to the 
TSS method developed by Vinnicombe et al. [9], In the case 
of the A320 GLAF the loads due to CT, in general, fell below 
the design envelope, so the problems of producing correlated 
load sets, associated with the method, were not unduly 
restrictive. Where correlated loads were required it was 
possible to consider only local parts of the structure and use 
approximate methods based on either the DSG or averaging 
of the load responses resulting from the TSS analysis. 

In both the above examples, the effect of the non-linearity 
was not known prior to undertaking the non-linear analysis. 


2-7 


In fact, it is not an easy task to predict what effect a particular 
non-linearity will have on the aircraft loads, i.e. it is not 
unreasonable to assume that a GLAF will reduce wing up¬ 
loads, since this is its purpose, however the effect of the 
system on, say, tail loads is far less predictable. It is quite 
often the case that a non-linearity will both reduce and 
increase loads compared to those predicted for a linear re¬ 
presentation of the system. Consequently, it is difficult to 
judge whether or not a linearised model provides, at least, a 
conservative estimation of all loads, so it is necessary to look 
closely at the problem. Even if a linearised model is con¬ 
servative for all loads it may be in the manufacturer's interest 
to exploit the non-linearity. This is clearly the case for gust 
load alleviation systems where the clear objective is to reduce 
the wing up-loads, hence the overall aircraft weight. What¬ 
ever the situation, there is clearly a need for the manufacturer 
to establish the effect of all control system non-linearities on 
the predicted loads. The procedure adopted by British Aero¬ 
space Airbus Limited has been: 

1. Linearise the problem and use RPT to perform a DEA to 
identify critical flight conditions for the linear model. 

2. Study the results and make a judgement on which non- 
linearities, if any, may have a significant effect on the 
aircraft response. 

3. Use a selection of methods, suitable for the non-linearity 
in question, to assess its influence. 

Note 5: for recent Airbus aircraft, use of the DSG method 
(see note 3) for a reduced set of "Interesting Quantities" 
has proved to give a good indication of the influence of 
the non-linearity and of the expected change in design 
loads. Both TSS and EG methods were also used for the 
autopilot non-linearity, eventually the EG method was 
selected for DEA. In the case of the A320, the TSS method 
was de\’eloped, see Vinnicombe et al. [9], and eventually 
used to demonstrate the effectiveness of the GLAF for 
certification. In all cases the EG, TSS and DSG demon¬ 
strated a similar influence on aircraft loads for the non¬ 
linearity in question. 

Note 6: if the existing "worst-case" analysis techniques 
have any place in future turbulence analysis it would be 
at this stage of the design process where the cost may be 
tolerable. They would, however, only be used to corrob¬ 
orate the validity of other methods. 

4. Determine whether or not there is a need to model the 
non-linearity and, if so, choose the most practical method 
for use in the ensuing DEA. 

8.0 CONCLUSION 

The impact of non-linear flight control systems on the pre¬ 
diction of aircraft loads due to turbulence has been two-fold: 

• There is an ambiguity regarding interpretation of 
the CT airworthiness requirements that has led to 
the development of an assortment of methods 
based on different assumptions. 

• The cost to the manufacturer of the loads pre¬ 
diction process has risen and is likely to continue 
to rise. 

It is not the intention to address the airworthiness require¬ 
ments within this paper, however, there is an urgent need 
both to eliminate the ambiguity and ta ensure the structural 
integrity of future designs. The CT DEA requirements were 


founded on the "principle of equivalent safety", so it would 
seem that a stochastic based analysis remains the most 
suitable for dealing with non-linear aircraft. On the other 
hand, it may be that a "worst-case" approach is preferred as 
the basis for a future requirement. If this is so, practical 
methods for predicting loads must first be established. In 
either case there is a need for clarification of the CT 
airworthiness requirements, for non-linear aircraft, and for 
the development of methods that will provide practical means 
of compliance. 

Stochastic based processes currently border on the limits of 
practicality, it may be expected however that advances in 
computing will make these methods practical for the design 
of the next generation of civil transport aircraft. Even so, 
there still remains the problems of correlated loads and the 
treatment of the turbulence intensity. 

Given the current "state-of-the-art" in computing it is unlikely 
that "worst-case" gust methods will provide a practical app¬ 
roach to loads prediction for the next generation of civil 
transport aircraft. "Worst-case" approaches lose the balance 
between complexity, accuracy and cost; the accuracy of the 
mathematical model and the model of the atmosphere do not 
warrant such complexity. 

The intention of this paper is to give an industrial view of the 
impact of non-linear flight control systems on the prediction 
of aircraft loads due to turbulence. It does not recommend a 
particular method for dealing with non-linear problems, on 
the contrary it is believed that the method used should be 
tailored to the problem in hand. 

It is hoped that the reader has gained an insight into both the 
problems encountered, when predicting aircraft loads due to 
turbulence for a non-linear aircraft model, and the way that 
these problems are currently addressed by British Aerospace 
Airbus Limited. At the same time there is a need for a review 
of the CT airworthiness requirements and for the continuing 
development of non-linear analysis techniques that provide 
practical means of compliance. Such methods must balance 
both complexity and computational effort against other 
modelling assumptions and external parameters. It should be 
remembered that the analyses discussed within this paper are 
for establishing limit loads and that there is a factor of one- 
and-a-half (1.5) between these and the ultimate design loads. 
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I. RESUME 

Cet article traite du calcul des charges 
structurales des parties motrices lorsque les 
caracteristiques inertielles, les vitesses de 
rotation elevees et les grands deplacements, 
induisent un impact important des efforts 
d'origine gyroscopique. Une methode directe, 
integree aux equations de la dynamique, pour 
le calcul des efforts gyroscopiques, en 
fonction des moments d'inertie polaires, des 
vitesses de rotation et des vitesses de 
precession est presentee et comparee a la 
procedure simplifiee de superposition sur la 
base des vitesses angulaires mesurees au 
centre de gravite moteur. Les resultats 
obtenus traduisent une bonne coherence entre 
les deux methodes de calcul, pour les moments 
de tangage et de lacet. Les fonctions de 
transfert presentees mettent en evidence 
1'importance des effets gyroscopiques sur le 
modele dynamique, d'ou un impact significatif 
sur l'ensemble des efforts moteur, aussi bien 
pour le jet que pour le turbopropulseur. 
D'autre part, l'impact des forces d'origine 
aerodynamique sur les helices, induites par 
la reponse des propulseurs aux perturbations 
atmospheriques, est examine. Les hypotheses 
retenues, la methode de calcul et les 
resultats obtenus sont presentes. 

I I . ABSTRACT 

this paper deals with the problem of 
calculating loads on engines and their 
supporting structures when inertial 
characteristic, high engine speeds, large 
displacements, imposed by extreme atmospheric 
disturbances, give an important impact of 
gyroscopic effect. A direct method 
integrating the gyroscopic action in 
Lagrange's equations as a function of the 
polar moment of inertia, engine revolution 
speed and precession speed is presented and 
compared to the more common practice of 
superimposing gyroscopic effects on the basis 
of the angular speed calculated at the engine 
center of gravity. The results establish a 
good correlation between the two method on 
pitch and yaw moments on engine. On the other 
hand, tranfert funtion modifications due to 
groscopic effects, generate an increase of 
other loads for the direct method 


calculations, for both the turboprop and the 
turbojet aircraft. Effects of the 
aerodynamic force of the propeller due to 
the engine response are also examined. 

2. INTRODUCTION 

Les multiples phenomenes interferants dans 
la zone des parties motrices des avions de 
transport, phenomenes thermiques et 
dynamiques, geometrie tres perturbee, siege 
de multiples interactions aerodynamiques 
entre les elements en presence: fuselage, 
nacelle, mat, voilure, deflecteurs, surfaces 
mobiles... etc, en font une zone dont les 
modelisations sont particulierement 
diffidles a apprehender. La quantite et la 
finesse des etudes menees tant par les 
motoristes que par les avionneurs en sont la 
preuve. 

Le calcul des efforts generaux pour la 
prevision des charges limites et extremes 
est tributaire de la modelisation mise en 
place, avec une sensibilite plus ou moins 
aigue aux variations des parametres 
caracteristiques. 

Ainsi, les dimensions importantes des 
ensembles propulsifs, les vitesses de 
rotations elevees, les structures de plus en 
plus souples, notamment la voilure, et les 
deplacements importants qu'elles impriment 
aux elements tournants, sont autant de 
parametres qui font que les moteurs et leurs 
elements de liaison avec la voilure sont 
amenes a subir des efforts gyroscopiques et 
aerodynamiques non negligeables pour leur 
dimensionnement. De ce fait, les 
reglementations imposent la prise en compte 
de ces effets dans de nombreux cas de vols, 
dont les cas de turbulences atmospheriques 
(§25.371 pour la FAR et le JAR). 

Dans ce cadre, cette etude presente 
1'influence des efforts inertiels 
gyroscopiques et des effets aerodynamiques 
des propulseurs, sur les charges 
structurales des parties motrices de deux 
avions de transport civil, un turbojet de 
type Airbus et un turbopropulseur de type 
ATR. 


Paper presented at a Workshop on "Aircraft Loads due to Turbulence and their Impact on 
Design and Certification” held in Lillehammer, Norway on 5 May 1994. 
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3. AVANT-PROPOS - MODELISATION 

Les avions traites etant de conception et de 
physionomie trbs differentes, les hypotheses 
et les modelisations adoptees le sont 
egalement. Dans le cas du turbopropulseur, le 
moteur etant en partie integre a la voilure, 
ce dernier ne fait l'objet d'aucune 
modelisation aerodynamique particuliere. II 
est represente par plusieurs masses 
comportant 6 degres de libertes chacune, dont 
le centre de gravite et le moyeu de l'helice, 
que l'on utilisera pour 1'integration des 
forces d'origines gyroscopique et 
aerodynamique provenant du propulseur. A 
1'oppose, les parties motrices du turbojet 
etant separees de la voilure par un mat, les 
nacelles sont modeiisees par un prisme de 
base hexagonale forme de 36 paves, permetta.nt 
la definition de forces aerodynamiques 
instationnaires suivant les axes vertical et 
lateral. Dans ce cas, les efforts 
gyroscopiques sont evalues au niveau du 
centre de gravite du moteur. Pour les deux 
avions traites, les forces aerodynamiques sur 
les surfaces modeiisees (voilure 
empennages,...) font l'objet d'un calcul 
instationnaire par la methode dite des 
rbseaux de doublets (Dubblet Lattice Method). 

Cette modelisation suppose done, pour les 
caracteristiques geometriques des motours, 
que le mouvement du point de reduction cos 
efforts (moyeu ou centre de gravite) soit 
representatif du mouvement des elements en 
rotation, e'est a dire que les differents 
composants du moteur soient indoformabies et 
que les axes de rotation restent aligr.es au 
cours du temps (pas de balourds). 

D'autre part, la schematisation utilisec est 
limitee au demi avion et done fondee sur une 
representation symetrique do la structure. De 
meme, seules les perturbations atmospheriques 
verticales et constantes en envergure sent 
traitees. Ceci represente une simplification 
importante compte tenu du probleme etudie. En 
effet, le sens de rotation identique cos 
moteurs gauche et droit, implioue a travers 
les effets gyroscopiques, une dissymetrie 
dans la reponse dynamique de la structure, 
dissymetrie non modelisee dans notre cas. 

Enfin, la representation modale do la 
structure prend en compte les moteurs a 
l'arret, e'est a dire leur masse inerte. Les 
effets inertiels de gyroscope sont done 
superposes & cette representation. 

La vitesse avion et le regime de rotation des 
moteurs (il) sont consideres commc constants 
lors du deroulement d'un calcul. 

4. CALCUL DES EFFORTS GYROSCOPIQUES 

4.1 M6thodes 

II est commun, pour le calcul des charges 
structurales, d'utiliser une representation 
modale de la structure, theoriaue ou vnlidce 
a partir d'essais de vibration au sol, les 
moteurs etant a l'arret. Les equations 


caracterisant la reponse dynamique de 
1'avion sont done resolues sans prise en 
compte des effets gyroscopiques, qui sont 
simplement superposes a partir de la 
relation classique suivante, fournissant le 
moment gyroscopique de solides en rotation : 
-* -* -+ 

MCyr (t) = Ip ft A fip(t) 

Ip etant le moment d'inertie polaire du 
solide en rotation, 

Qsa vitesse angulaire (regimes moteur), 
£2p(t) la vitesse angulaire de precession 

Une methode plus precise, dont la theorie 
utilisee est presentee par Beatrix et Dat 
[2], permet la prise en compte directe des 
termes gyroscopiques dans les equations de 
Lagrange. L'integration de ces moments au 
calcul de la reponse dynamique de 1'avion 
souple a une rafale atmospherique est 
brievement rappeie dans ce qui suit. 

Compte tenu des hypotheses formuiees, les 
parties tournantes sont uniquement 
caracterisees par leur energie cinetique, 
dont la derivee s'ecrit : 



La representation modale de la structure 
etant la meme que pour la methode initiale, 
les termes d'inertie classiques des parties 
tournantes sont deja integres et il reste 
done a ajouter a la modelisation le terme 
lie a la rotation et done aux effets 

gyroscopiques : [g] 

[ G 1 est une matrice carree antisymetrique 
de coefficients constants sur les lignes et 
les colonnes correspondant aux degres de 
libertes en tangage et lacet des points de 
reductions des deplacements de chaque 
element tournant. 
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La formulation de la derivee par rapport au 
temps du vecteur des deplacements en 
fonction de la matrice des deform^es 
unitaires modales [O] permet l'expression 
du terme gyroscopique dans la base modale et 
i'ecriture matricielle du mouvement de 
l'avion dans la rafale est done la suivante: 

[ K ] { } + [ [°] + [<t>] [g] [<d] ] {q} 

+ [ R ]{q} = [<t>] { Fraf } 

(les matrices [M],[C] et [K] contiennent les 
termes aerodynamiques de mouvement) 
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4.2 R6sultats 

Ces deux methodes, dites "par superposition 
effets gyroscopiques rajoutes" et "methode a 
calcul integre", ont ete utilisees pour 
devaluation des efforts gyroscopiques de 
deux avions, representatifs de leur 
categories respectives, un turbopropulseur 
court-courrier de type ATR et un turbojet 
moyen courrier double couloir de type Airbus. 
D'autre part, les deux formes de 
perturbations atmospheriques exigees par les 
reglementations, la rafale discrete accordee 
et la turbulence continue ont ete appliquees. 

4.2.1 Resultats obtenus pour le turbojet 
Une etude complete a ete menee, fondee sur 
les cas definis lors des phases de conception 
et de certification de cet avion. Les charges 
moteur, au niveau du centre de gravite, sont 
presentees, dans les cas de poussee maximale 
continue, de poussee nulle et moteur a 
1'arret. Pour chaque effort et acceleration, 
sont comparees les valeurs critiques obtenues 
dans chacun des cas de poussee, les effets 
gyroscopiques pouvant modifier la hierarchie 


des cas traites. Seuls 1'increment d'effort 
du a la perturbation atmospherique est pris 
en compte a 1'exclusion des charges dues au 
vol stationnaire en palier a Nz=l. 

Les graphiques suivants presentent les 
resultats obtenus pour la rafale discrete par 
la methode a calcul integre. Les moteurs 
gauche et droit sont uniquement differencies 
par 1'inversion du sens de rotation, les 
interactions entre les deux cotes n'etant pas 
consideree. 

L'effet du couplage de la reponse dynamique 
avec le sens de 1'effort gyroscopique 
apparalt a travers la comparaison de ces deux 
graphiques. On constate dans la majorite des 
cas, que les effets gyroscopiques modifient 
assez peu le comportement dynamique des 
moteurs de cet avion, et se superpose (>0 ou 
<0) au calcul moteur a 1'arret. Seuls les 
efforts Fz, Mx et Mz mettent en evidence un 
couplage tres different entre les deux cotes, 
notamment a travers 1"apparition d'un nouveau 
cas critique (n°l, 2 ou 3 dans les figures). 
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D'autre part, pres de la moitie des criteres 
(efforts ou accelerations) voient leur cas 
critiques modifies par 1' integration des 
efforts gyroscopiques. 

Une analyse identique a ete effectuee dans le 
cas de l'excitation par la turbulence 
continue. Les graphiques suivants mettent en 
Evidence un impact globalement moins 
important que pour la rafale discrete, limite 
dans ce cas a environ +5% sur le moteur 
gauche. On notera que la turbulence continue 
reprfesente les cas dimensionnant pour la 


structure des parties motrices de cet avion. 

Les resultats moteur droit sont compares 
pour le cas a poussee maximale, aux valeurs 
pour le moteur gauche. 

On remarque egalement le faible impact des 
cas a poussee nulle, ceci dO notamment aux 
regimes moteurs relativement faibles. 

D'autre part, 1' influence des effets 
gyroscopiques moteur sur les charges subies 
par la voilure et le fuselage a 6te evalu^e. 
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On note que cet impact, resultant du couplage 
entre les modes souples structuraux des 
parties motrices et ceux de la voilure et du 
fuselage s' avere tres faible. On constate 
une variation des charges inferieure a 0,2% 
sur la quasi- totalite de la voilure et du 
fuselage. Seul le fuselage avant fait 
exception, ou le facteur de charge vertical 
augmente de 0,9%. La verification de cette 
hypothese nous a permis de travailler a 
partir d'un modele symetrique limite au demi- 
avion, chaque cote etant caract6ris6 par le 
sens donne a la rotation des moteurs. 

Un point essentiel mis en evidence par ces 
resultats, est l'apport procure par 
1 1 amplification des couplages modaux, sur les 
charges non soumises directement aux effets 
gyroscopiques. L'effort et 1'acceleration 
laterale (Fy et Ny) ainsi que le moment de 
roulie Mx pour la turbulence et l'effort et 
1'acceleration verticale (Fz et Nz) ainsi que 
le moment de roulis Mx pour la rafale en sont 
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les meilleurs exemples. Non pris en compte 
par la methode dite "par superposition 
efforts gyroscopiques rajoutes", on trouve 
des augmentations de l'ordre de 5 a 10%. 

En ce qui concerne les charges directement 
associees aux efforts gyroscopiques, les 
moments de tangage My et de lacet Mz, les 
graphiques suivants presentent des tendances 
variables. Neanmoins, le faible impact sur le 
moment My et des niveaux enveloppes voisins 
pour le moment Mz, enregistres pour la 
methode "directe" et pour la methode "par 
superposition", traduisent une assez bonne 
coherence entre ces deux methodes de calcul 


pour les avions traites, a la fois pour la 
turbulence continue et la rafale discrete. 
Dans le cas de la turbulence continue, le 
calcul des efforts gyroscopiques par la 
methode dite "par superposition" A ete 
effectue a partir de la relation suivante: 

■» 4 4 

CMGyr = Ip Pm£2 p n A on? 

G represente l'ecart type de la grandeur et 

Pm£2p le coeff icient de correlation entre la 
composante M de moment considere et la 
vitesse angulaire f2p correspondante. 
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4.2.2 ResOltats obtenus pour le 
turbopropulseur 

L'6tude men6e sur cet avion, est fondee sur 
les cas determines comme critiques au niveau 
du centre de gravite moteur, lors de la 
certification de cet avion. Les efforts 
g6n6raux au CdG moteur sont presentes dans 
les deux cas : moteur a l'arret, c'est a dire 
sans effets gyroscopiques et pour le regime 
100% (environ 1200 tours/mn), legerement 
sup6rieur au regime de croisiere egal a 86%. 
Les effets gyroscopiques induits par le 
propulseur sont integres dans le modele, au 
niveau du moyeu de l'helice. 

De meme que pour le turbojet, les 
comparaisons presentees portent uniquement 
sur 1'increment de charge induit par la 
perturbation atmospherique, a l'exclusion des 
charges dues au vol en palier. 

Pour le cas de la rafale discrete, les 
r6sultats obtenus par la methode dite "a 
calcul integre" pour les 4 premiers efforts 
(Tx, Ty, Tz, et Mx), et les deux sens de 
rotation des moteurs, seul l'effort lateral 
Fy met en evidence une sensibilite 
significative aux effets gyroscopiques, avec 
plus de 15% d 1 augmentation sur le moteur 
gauche, et un comportement symetrique pour le 
cot6 droit, les autres efforts ne presentant 
pas d'impact superieur a 2%. 

Les comparaisons sur les moments de tangage 
et de lacet sont presentes par les deux 
graphiques suivants, pour les deux sens de 
rotation et pour les deux methodes de calcul 
"integre” et "effets gyroscopiques rajoutes". 
On ne constate aucun impact significatif sur 
le moment de tangage My, et a 1'oppose, un 
impact majeur sur le moment de lacet Mz. Les 
deux methodes de calcul sont coherentes aussi 
bien sur la tendance generate de 1'impact que 
sur les niveaux rencontres. 


Les r£sultats obtenus pour cet avion, sont 
tres differents de ceux presentes pour le 
turbojet, principalement du fait des 
differences de conception de 1'assemblage 
moteur-mat-voilure. Le turbopropulseur se 
caracterise, comme la plupart des avions de 
ce type, par 1■integration du moteur a la 
voilure. Les mouvements du moteur induits par 
les deformations de la voilure sans flfeche, 
sont done essentiellement dans le plan 
vertical (deplacement suivant z et rotation 
en tangage). De ce fait, le mouvement 
coplanaire de la voilure, associ£ aux 
deformations du mat-moteur sont les sources 
essentielles de rotation de lacet au niveau 
du moteur. La difference d'impact entre le 
deux moments My et Mz, s'expliquent done de 
fagon assez classique, par le fait que la 
part gyroscopique My g induite par rotation de 

lacet a z du moteur, d 1 amplitude tres faible, 
est negligeable devant 1'importance des 
efforts Tz et My resultant de la reponse de 
la voilure en flexion/torsion et de l'avion 
en pompage-tangage. A 1'inverse, le moment 
gyroscopique Mz g induit par ces mouvements 
de tangage est preponderant devant le efforts 
Fy et surtout Mz d'origine principalement 
inertielle, resultant du mouvement de la 
voilure en coplanaire ou des deformation du 
mat-moteur en lacet. 

La mise en evidence du contenu fr6quentiel de 
la reponse dynamique est realisee par les 
figures en pages suivantes, sur la base de 
reponses a plusieurs rafales discretes de 
longueurs d’onde differentes mais de memes 
intensites, pour les efforts Tx max, Ty max, 
et Mz min. On constate des modifications dans 
le couplage des modes structuraux, se 
traduisant par des reductions (cas Tx) ou des 
augmentations de niveaux des charges (cas Ty 
et Mz) et par une modification de la 
frequence critique (cas Mz). 
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ENGINE YAW MOMENT FOR DIFFERENT DISCRETE GUST GRADIENT DISTANCES 


I 


(%) 



De meme, les resultats obtenus sous une 
excitation par la turbulence continue sont 
present6s. D'une maniere generate, les 
tendances ressenties pour la rafale discrete 
se confirment, notamment pour l'effort Ty 
(+8% cot6 gauche) et les efforts Tx, Tz, Mx, 
et My pour lesquels on rencontre peu ou pas 
d'impact, comme le montre la figure suivante. 
La principale difference se situe pour le 
moment Mz, pour lequel on ne retrouve pas 
1’amplitude de 1'impact mis en evidence par 
les r6ponses a la rafale discrete. Le 
caractere spectral de la turbulence continue 
a certainement une influence importante dans 


ce resultat, et les fonctions de transfert 
pour les efforts Ty et Mz sont presentees. 
On remarque notamment la presence des modes 
de mouvement lateral du moteur (lateral 
engine~5Hz), de coplanaire voilure (l st wing 
fore-aft~7,8Hz), et l'apparition du mode de 
tangage moteur (engine pitch~10Hz) dans les 
cas avec effets gyroscopiques pour l'effort 
Mz. De meme la participation du mode de 
lacet moteur (engine yaw~12Hz) est remplac6 
par le mode suivant, tres similaire, mais 
comportant ses deformees laterales (lacet 
moteur et coplanaire voilure) en opposition 
de phase. 



fore/aft lateral vertical roll pitch yaw 

load load load moment moment moment 


Hleft engine max thrust □ no gyro effect 0 right engine max thrust 
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5. FORCES AERODYNAMIQUES DES 
PROPOLSEURS 

La theorie utilisee, est celle presentee par 
Houbolt et Reed [1], permettant le calcul 
quasi-stationnaire des efforts aerodynamiques 
sur les pales. Cette formulation est enrichie 
de l'effet de compressibility de l'air et 
d'un effet instationnaire, modelis^ par la 
fonction de Theodorsen, qui traduit le 
dephasage entre la variation d'angle 
d'attaque et la force aerodynamique 
correspondante sur la pale. Ces efforts 
rdsultent de la rotation du plan de l'helice 
en tangage et en lacet, induite par les 
mouvements provenants de la reponse dynamique 
de la voilure, du mat ou du moteur aux 
perturbations exterieures. 

Les hypotheses et la formulation retenues 
permettent 1'integration de ces termes aux 
equations de la dynamique de 1'avion en vol 
exprimees dans la base modale. 

5.1 Mdthode 

Pratiquement, les efforts sont calcules a 
partir de 1'evaluation en ecoulement 
compressible de l'effort local sur une 
tranche de pale en fonction de la vitesse de 
l'air et du point de vol. La vitesse de l'air 
sur la pale, consideree comme rigide des lors 
qu'elle n'apparait pas dans le schema modal 
utilise, est la composition de la rotation de 
l'helice, de la vitesse avion et des vitesses 
induites par les reponses dynamiques de la 
voilure, du mat ou du moteur. 

La figure suivante propose un schema de cette 
modelisation. 

De nombreux effets sont done negliges, tels 
que l'impact des rafales atmospheriques sur 
les pales, les dissymetries, la variation du 
regime de rotation, de 1'angle d'attaque des 
pales et des efforts de traction, les 
balourds, etc... 


[ K' ' ' ] =S 


Le probleme se pose done d'intdgrer ces 
efforts aerodynamiques au calcul de la 
reponse dynamique de 1'avion souple a une 
rafale atmospherique. 

La vitesse de l'air sur une tranche de pale, 
etant exprimee en fonction du rayon local, 
des parametres de position du moyeu et de 
leurs derivees, la limitation aux termes du 
premier ordre permet la mise en place d'une 
expression relativement directe des efforts 
et moments elementaires sur la tranche. La 
prise en compte des differentes pales 
(notamment les effets de symetrie autour de 
1'axe moteur) et 1'integration des efforts 
locaux sur le rayon des pales, permettent 
l'expression de l'effort global sous la forme 
d'une relation matricielle directement 
integrable au second membre des equations de 
Lagrange (travaux de Beatrix et Dat [2]) : 



Les matrices de coefficients [K'''3 et [K'3 
contiennent, pour les degres de libertes du 
moyeu de l'helice y, z, 0, 4* (respectivement 
lateral, vertical, tangage et lacet), des 
termes que 1'on peut rapporter aux gradients 
aerodynamiques correspondants de l'helice : 

degres de liberte du moyeu moteur 
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q etant la pression dynamique (DaN/m 2 ), 
D le diametre de l'helice (m) et 
S sa surface (m 2 ) 


(X, Y, Z) repere lie a l'avion 
(U, V, Z) repere lie au moteur, plan 
de l'helice sous un deplacement 0 et 
l|/ respectivement tangage et lacet. 


Tangle d'htflice 
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Les coefficients de ces deux matrices 
s'expriment de la fagon suivante : 

Czq = -(4 £2 Cr/V) II ; Cz q = (4 £2 Cr/V) J2 

Cyg = -(4 £2 Cr/V) J1 ; Crrg = -(2 £2 Cr/V) J2 

Cy q = -(4 £2 Cr /V) 12 . Cm q = -(2 £2 Cr/V) 13 

Cng = -(2 £2 Cr /V) 12 . Cn q = -(2 £2 Cr/V) J3 

De plus, les symetries des helices autour de 
leur axe de rotation, permettent d'ecrire les 
relations suivantes : 

c y>y = - Cz g 

CZy = Cyg 
C Y r = ~ Cz q 

Cz r — Cy q 

Les termes en gras dans les matrices [K*] et 
[K'"] proviennent de l'effet purement 
instationnaire, partie imaginaire de la 
fonction complexe de Theodorsen : F(k)+iG(k). 
Dans l'expression des gradients 
aerodynamiques, Cr est la corde dite "de 
reference" de la pale et V la vitesse de 
l'avion. Les termes II, Jl, 12, J2, 13 et J3, 
realisent 1■integration sur l'ensemble de la 
pale des efforts elementaires a partir des 
rayons courants adimensionnels (Tl=r/R) et de 
la corde locale C (T|) (m) : 

N£ C1CC ^2 Ar r cm) F(k) ^ 

4 271 Cr H (Tl, Mach, Ar) 

■Ho 

1 5 

N£ Clot ^ Ar r T\ z ccrp F(k) dT] 

4 271 Cr J H (tl, Mach, Ar) 

^0 
1 4 

N£ Cla Ar r Tp cap FOO dT) 

4 271 < - r H (Tl , Mach, Ar) 

■Ho 


il = 


12 = 


13 = 


Cmy - Cng 
Cn^ = Cmg 
Cm r = - Cn q 
Cn r = Cm q 


Np etant le nombre de pales de l'helice. 
L'effet de compressibilite est pris en 
compte a travers la fonction H (TJ, Mach, Ar) 
dependant du nombre de Mach, de 
1 1 allongement de la pale (Ar) et de la 
position en envergure (T|) de la tranche 
consideree : 

H (Tl, Mach, Ar) = ~~ a/ ( v 2 +£2 2 r 2 t| 2 ) * 

[ 2 + Ar VI - M 2 /V 2 (V 2 + £2 2 R 2 T| 2 ) ] 

en remarquant que le terme entre parantheses 
represente le carre de la vitesse locale 
tangente a la pale dans le repere de 
l'helice, et V 2 /M 2 le carre de la celerite 
du son. 

5.2 R6sultats 

L"evaluation de 1'impact des effets 
aerodynamiques a ete realisee dans les memes 
conditions que celles decrites pour les 
effets gyroscopiques, limitees au cas du 
turbopropulseur, la methode utilisee n'etant 
pas adaptee aux fans ou aux helices 
carenees. Les resultats obtenus pour la 
turbulence continue sont presentes figure 
suivante. On remarque la variation 
importante des efforts Ty et Mz. La 
modification de comportement vis a vis du 
mode lateral moteur (lateral engine~5Hz) est 
mise en evidence sur les fonctions de 
transfert page suivante. Cet effet a ete 
enregistre pour la plupart des cas traites 
(voilure vide ou pleine de carburant). A 
1'oppose, un impact tres faible, inferieure 
a 3%, a ete releve pour la rafale discrete. 

Les effets de compressibilite de l'air et 
d'instationnarite, ne modifient que tres 
faiblement ces resultats. 



fore/aft lateral vertical roll pitch yaw 

load load load moment moment moment 
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6. CONCLUSION 

Les resultats obtenus sur l'impact des effets 
gyroscopiques sont significatifs, tant par les 
variations induites sur les charges moteur que 
par la modification du contenu en frequence 
des reponses dynamiques calculees, ceci dans 
des proportions tres liee a l'avion considere. 
En ce qui concerne les moments de roulis et de 
lacet. My et Mz, malgre les differences dans 
la prise en compte de la reponse dynamique des 
moteurs, la comparaison entre la methode dite 
"par superposition" et la methode integrant 
ces effets aux equations de la dynamique met 
en Evidence des niveaux enveloppes globalement 
voisins ce qui traduit neanmoins, pour les 
avions studies, une assez bonne coherence 
entre ces deux methodes. On note egalement un 
impact induit parfois important sur les autres 
efforts moteur (notamment pour 1'effort 
lateral Ty). La prevision des charges reste 
delicate, du fait de la grande sensibility des 
resultats aux configurations massiques 
traitees, a 1'integration motrice et aux 
conditions de couplage. 

En ce qui concerne la seconde partie, 
traitant des forces aerodynamiques sur les 
pales d'helices, l'impact mis en evidence est 
g6n6ralement assez faible. On remarque 
neanmoins la modification de comportement du 
mode de deplacement lateral du moteur du 
turbopropulseur impliquant une augmentation 
significative des charges Ty et Mz dans le 
cas de la turbulence continue, augmentation 
non retrouvee dans pour la rafale discrete. 

On prolongement a cette etude serait la mise 
en ceuvre des points suivants. Tout d'abord. 


un modele asymetrique avion complet 
apporterait la prise en compte des effets 
antisymetriques entre les deux moteurs, done 
des couplages differents negliges ici. De 
plus, les variations d'incidence induites par 
les rafales atmospheriques sont egalement 
generatrices de charges sur les pales, et 
1'evaluation de leur impact reste a 
effectuer. 

Au dela de ces aspects, il doit etre 
souligne, que la prevision des charges sur 
les parties motrices reste d'une complexity 
considerable, tant du fait de la modelisation 
de 1 1 aerodynamiques locale, que de 
1'hypersensibilite des charges au modele 
structural utilise. II apparait egalement 
que, compte tenu de la complexity 
suplementaire que les effets gyroscopiques ou 
aerodynamiques moteur, ajoutent au modele de 
calcul, et a l'analyse des resultats, leur 
integration doit etre effectuee avec prudence 
a partir d'une connaissance suffisante des 
caracteristiques dynamiques du nouvel avion. 
La maitrise de ces problemes par la 
modelisation utilisee actuellement a 
1'Aerospatiale, associee aux procedures de 
dimensionnement, a montrer assurer un bon 
niveau de security aux elements ryalisant la 
liaison des parties motrices sur les avions. 
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Treatment of Non—Linear Systems by Timeplane—Transformed C.T. Methods 

The Spectral Gust Method 

H. Lusebrink, J. Brink—Spalink 
Deutsche Aerospace Airbus GmbH — Structural Dynamics 
P.O. Box 95 01 09, D—21111 Hamburg, Germany 


1. SUMMARY 

The present paper discusses the problems, which occur when 
PSD/Continuous Thrbulence (C.T.) gust methods, originally 
developed for linear dynamic gust analysis in the frequency 
plane, are applied to modern transport A/C with flight control 
systems, containing a variety of nonlinearities for signal 
shaping. 

Power spectral density (PSD) gust velocity intensities are 
defined by the von Karman power spectrum in the frequency 
plane, scaled by the design gust intensity U a (RMS). This 
definition of gust in the frequency plane causes the basic 
problems for nonlinear dynamic analysis, in contrary to the 
timeplane defined Timed Discrete Gusts. 

Discrete Gust calculations show that aircraft nonlinearities 
cannot be neglected, lb treat nonlinearities for C.T. there are 
basically three approaches (apart from linearization methods): 


(SSB) Stochastic Simulation Based Methods 

(timeplane stochastic simulation), usually in combination 
with exceedance rate based design load definition 

(MFB) Matched Filter Based Methods (’Worst Case Gusts’): 

For each load station, a search within a class of shape 
functions of prescribed ’spectral energy’ is necessary 
to find that shape giving the highest load. 

(SG) Spectral Gust Method: 

Uses discrete gusts having exactly the von Karman 
spectrum, and computes the energy norm of the load time 
history. 

A practical method for nonlinear C.T. method should have the 
following characteristics: 

-give loads consistent to existing requirements (for 

linear aircraft), 

— give correlated (balanced) loads 

-straight-forward and economical, when applied for 

a large number of design cases to be calculated. 

Since the SSB methods need long simulation times and give no 
balanced loads (together with exceedance rate counting), and the 
MFB methods are very expensive, the present paper proposes 
the SG method. It needs only one step to derive design loads, is 
consistent to Design Envelope Analysis and provides correlated 
loads. 


2. INTRODUCTION 

Continuous Thrbulence (C.T.) gust design methods were 
introduced in the requirements ( FAR/App. G, JAR/ACJ 
25.305(d)), because the Discrete Gusts (DG) represented a too 
artificial description of atmosphere, not taking into account its 


randomness. The C.T. methods, first the Mission analysis 
(MA), later the simpler Design Envelope Analysis (DEA) are 
based on (von Karman) power spectral density description of 
gust velocity, and used frequency domain techniques to compute 
the loads. 

The first problem with C.T. methods was to adjust the design 
gust intensity parameters: bj, b 2 , P 2 (for MA) or U„ 
(for DEA) respectively, based on A/C designs with good 
structural safety records. 

The next problem was the computation of balanced (correlated) 
loads, a demand from stress offices. Meantime this was solved, 
but only for DEA. 

The third problem with C.T. methods since many years is how to 
treat aircraft nonlinearities. While aerodynamic nonlinearities 
can be linearized in a conservative manner, the control system 
nonlinearities (limiters, thresholds) of modern transport aircraft 
are more difficult to handle. Statistical linearization methods are 
not straight-forward. They need a two-step analysis 
procedure: First calculate Equivalent Gains, second: calculate 
design loads.The method is expensive and time consuming, if 
several such nonlinearities exist. The most adequate way to treat 
nonlinearities seems to be simulation in time domain. But time 
domain methods for C.T. introduce three new analytical 
problems: 

-Unsteady aerodynamics have to be translated into time 

domain. 

-Random gusts have to be realized in time domain or 

to be replaced by ’equivalent’ discrete gust time 
histories. 

-The specification of ’design load’ from the 

corresponding load time history consistent with linear 
methods. 


3. PROBLEM DESCRIPTION 

The problem is the definition and computation of design loads 
for nonlinear aircraft subjected to random gusts. 

Note that gust models, defined by Spectral Densities and 
intensity parameters and their time—plane realizations, are the 
most representative description of the atmosphere. 

The nonlinearities can be found in: 

-Controls (limiters,thresholds, actuator dynamics), 

-Flight Mechanics (large rigid body motions) 

-Aerodynamics (high angles of attack) 

linearization usually gives conservative loads 

— Structure (free play, material,...); 
usually not important for loads 


Paper presented at a Workshop on “Aircraft Loads due to Turbulence and their Impact on 
Design and Certification ” held in Lillehammer, Norway on 5 May 1994. 
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Resulting problems 

Correct treatment of nonlinearities requires time domain 
methods, leading to following analytical difficulties 

-Unsteady aerodynamics have to be translated into 

time domain. 

Solutions: [14] 

-rational function approximation (s—plane) 

-z—transform approach 

— Random gusts have to be realized in time domain. 
Solutions: 

— Stochastic simulation 

-Worst Case Gusts (MFB, SDG,...) 

-Spectral Gusts 

-The definition of ’design load ’ for nonlinear 

systems needs further discussion. 

Solutions: 

-- RMS load 
-MFB load 

— Exceedance probability based load 
-Spectral Gust based load 


4. DISCUSSION OF POSSIBLE SOLUTIONS 
(Lin) Linearization 

This is the usual method today, when PSD methods ( DEA or 
Mission Analysis) are applied. Linearization is achieved in case 
of saturation nonlinearities in control laws (yaw damper, load 
alleviation systems) by reducing the gains (conservative). In 
more general cases with isolated nonlinearities statistical 
linearization (’’Equivalent Gain”—Method) is applied, which is 
appropriate for random gust . 

Disadvantage: Not straight-forward, time consuming and 
expensive 


(TSS) Timeplane Stochastic Simulation 

A random gust shape u(t) is generated using 

-filtered (Gaussian) white noise, or 

-superposition of sinewaves with random phase 

and von Karman amplitudes and subsequent 
Inverse Fourier transform. 

Then the nonlinear system is simulated in time domain (using 
rational function approximation of unsteady aerodynamics). 

Disadvantage: high cost (long simulation times); balance ? 


(MFB) Matched Filter Based Methods (NASA) 

For each load quantity a ’’worst case gust” shape u(t) has to be 
found by evaluating 

y* s = % max max, | y (t)| 

HLt ^ 2 

where the ’’spectral energy” 

II u ||* = \ j I G(im) | 2 /<J>((o) dm 
0 


(which is a measure for the ’’relative probability” e [4]) 
is constrained. 

The MFB methods are based on a variational principle 
(Cauchy—Schwarz Inequality), which for linear systems 
translates the RMS —norm of the DEA method into the 
max—norm of a load time history. 

The MFB—1 —Dimensional Search method [8] and the 
Statistical Discrete Gust method (SDG) are also based on this 
maximum principle. 

Disadvantage: high cost (optimization for each load station is 
required) 

(e.g. up to : 60 stations, each 12 responses; 

to be done for : 40 mass cases, 10 velocities/altitudes) 


(FPE) Fokker Planck Equation 

The FPE represents a partial differential equation in the state 
variables (and t) for statistical properties of a stochastic 
differential equation [ 5, 6]. 

Disadvantage: extremely high cost 


(Est) Estimation methods 


Determine scaling factors = f rorn discrete gusts 

and apply them to linear PSD results. 

Disadvantage: unknown accuracy, balance ? 


(SG) Spectral Gust Method (DA proposal) 

By a ’’Spectral Gust” we understand a discrete gust u G (t) with 
correct ’spectral density’, i.e. simply the scaled impulse 
response of a gust filter G(s): 

u G (t) = yitT L 1 G(i(o), . 

T (= L/V) being a reference time. Then we define the design 
gust load by an energy norm of the load response y(t) due to a 
design Spectral Gust u = U Q u G (t): 

y des.sG =11 y lit / </t , y(t) = y[U 0 u G ](t), 


where ||y||, 


i 


1/2 


I y(01 2 dt 


and U 0 (= o w ri d ) 


Lo J 

the design gust velocity from the DEA method. 


is 


Then for linear systems it can be shown that the SG method is 
totally equivalent to the (balanced) DEA method in the 
frequency plane. Hence the SG method is a direct translation of 
the DEA method into time domain, which requires no 
optimization. We discuss the SG method in more detail in the 
next section. 

Remark: 

We have || u c ||, = /f, and the natural selection T = L/V, 
t = Vt/L, implies || U o u 0 || t = U G and y des SG =|| y || t . 
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5. THE SPECTRAL GUST METHOD 

It is well known that the discrete (1 - cos) gust does not reflect 
correctly the power spectral density of atmosphere, see figure 1. 



Here, for reasons of comparison the discrete gust 

u(t) = j (1 - cos(Qt)), £2 = jtV/H, H = GustGradient, 
was periodically continued to get a nonvanishing power spectral 
density , which is defined for u(t) by [2] : 

+T 

<E>(co) = lim 4p | I u(t) e _imt dt | 2 

T=>« 31 A I 
= 2 * 

This inconsistency of the spectral densities is one reason why 
Discrete Gust calculations and PSD methods cannot give 
comparable results. 

In the following we derive ’’Spectral Gusts” u(t), i.e. discrete 
gusts having correct PSD (and the Laplace transform should 
exist). 

LetT be the corresponding reference time (e.g. T=L/V, L the 
scale of turbulence). Then we require 

=^I U MI 2 . 

00 

where <b 1 (to) is the given PSD, normalized by j dijdw = 1. 

o 

Now let G(s) be a normalized gust filter for d>, i.e. a transfer 
function which generates random gusts with prescribed PSD 
from white noise. G(s) is a holomorphic function satisfying 
4>i(tn) = | G(ico) | 2 , GOO = G ( s )- 

Hence a normalized Spectral Gust is given by 
U G (itt>) = ,/irT G(iu>), and in time domain by 

u G ( t ) = m l -1 g(s) . 

It is simply a scaled impulse response of the gust filter G(s). Its 
energy norm is || u G || t = /f, || u G || t/T = 1, whereas the gust 
filter is normalized in frequency domain : || G || m = 1 

It should be remarked that there exist many gust filters; they are 
obtained from a special one by multiplication with a real, 
odd, analytic phase function rji(a)) : 

G^(ico) = G(ito) eW"). 


Von Karman Spectral Gusts: 


We apply this procedure to the normalized von Karman power 
spectrum 

* k((0) = JL ^ + f (1-339 coL/V) 2 
nV [1 + (1.339o)L/V) 2 ] U/6 

Here L= 2500 ft is the scale of turbulence. We found by means 
of Beta-Functions that the constant 1.339 can be replaced by 

itwnU ' 1 - 3389852W - 

[ Brink-Spalink, Sept. 1992], 


Von Karman gust filters are given by 

fT 1 ± JH7$ sT 0 j 

T *- 1 ' 339 'fc 

considering the different signs. Note that a minus sign in the 
denominator would give a time reversed impulse response (filter 
is not realizable). It is possible to transform these gust filters 
analytically into time domain [ Lusebrink, Aug. 1990 ] : The 
denominator of form (s — a) -n has an inverse Laplace 
transform, and the numerator represents a differential operator. 
The resulting Spectral Gusts are quite simple: 

u g(‘) = ± c o . e-VTo • (^-‘Z 6 • (1 - c^), 

with constants c 0 = |^/| / T(^), c t = ± 

r(ll/6) = 0.940655858, andT 0 = 1.339L/V, whereas T 
(=L/V) is the Spectral Gust reference time. 

The singularity t _I/6 is very mild; it can be approximated with 
arbitrary accuracy. Note that the impulse response of a Dryden 
filter or a rational function approximated von Karman filter 
does not have a singularity. 

In case of a von Karman gust filter with an additional phase 
factor e nt ’(“ ) , an analytical Inverse Fourier transform will not be 
possible in general. On the other hand a numerical transform will 
converge extremely slow, because |G(ico)| a <d~ 5/6 for the 
higher frequency band. But convergence can be accelerated by 
low pass filtering and numerical differentiation in time domain: 

g(.)-(i + k !) L _1 {ttt? °< s )} • 


Mathematical Basis of the Spectral Gust Method 

For linear systems with transfer function H y (s) the design load 
as used in DEA is defined by 

r -i 1/2 


y<ies,RMS — U D A y , 


r -i 1 / 2 

oo 

j | H y j 2 4> K do) . 

_o 


If G(ito) denotes a normalized gust filter, | G | 2 = <f>i(= <1> K ), 
then is the energy norm in frequency domain of the combined 


transfer function H y G: 


*;= |J |H y G| 2 dmJ =|| H y G ||„. 

To translate this into time domain we remember Parseval’s 
theorem; its general form reads 


u(t) • v(t) dt 


-xfw 


U(iro) V(ico) dco , 


where U(s), V(s) are the Laplace transforms of u(t), v(t). We get 


|HyG|| ra = 


HyG HyGdcn 


I hy * g 1 2 dt 


= yit II hy * g III 

where (h y * g)(t) denotes the convolution of the impulse 
responses of H y (s) and G(s) . Hence 

y des,RMS = D a || H y G || M = U 0 || h y * g || t . 
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For linear systems we may introduce a scaling factor c : 

YdeSjRMS = l|hy*(c g) !|. - 

c corresponds to the reference time T as derived above: 

c = U 0V /jtT. NowfornonlinearsystemswedefinetheSpectral 
Gust design load by 

ydes,SC> =lly[UoU G ](t)||, / /f , 

where u G (t) = /itT L _ ^ G(s) is a Spectral Gust and T 
(= L / V ) is a reference time . 

Note that for linear systems this definition is independent on 
U 0 , and independent on the selected gust filter G(s). For linear 
systems it is totally equivalent to the DEA method in the 
frequency plane. Moreover, also the correlated loads from the 
DEA method can be defined in time domain, using Parseval’s 
theorem once more: 

Zy : = Zd's Qy: , wh CrC 
00 

j y(t) z(t) dt 

Qyz = ° II y lit INI, 

is a correlation coefficient in time domain. 

C. T. Design Load Definitions for nonlinear systems 

In case of nonlinear systems the above definition will be 
dependent on the selected Spectral Gust. On the other hand no 
agreed definition exists in this case for random gusts. Possible 
definitions are: 

— the RMS gust load (used in the DEA method) 
y dcs = RMS ( y [U D u K (t)] ] 

with u K (t) = von Karman — filtered Gaussian unit 
white noise . 


the MFB gust load (prescribed ’gust probability’) 
y des max max, |y(t)| , c = U 0 , 

IHill . < 


load with prescribed exceedance probability 
Nyfyda) = N des (used in Mission Analysis) , 
where y(t) = y [u K (t) • U 0 /rj d ], 

the Spectral Gust based design load : 
yd CS ,so =1 y III / v/T, y(t) = y[U o u 0 ](t), 


6. COMPARISON OF LOADS FROM 
VARIOUS METHODS 

The nonlinear test model (open loop) 

The best way to demonstrate the behaviour of the various 
nonlinear gust loads methods is to apply them to a nonlinear 
system, which is simple enough to compute results analytically, 
but complex enough to represent a broad class of systems 
(eigenfrequencies, dampings, nonlinearities). Our test model 
consists of a linear transfer function (harmonic oscillator, 
including load recovery) followed by a static quadratic 
nonlinearity (progressive) : 


V 2 + b o 5 + c o 


v . n 4- ivh 

+ 2 PqWq s + (i)q 



x l 1 ^ l x U 


Hence the load response y(t) due to a general gust u(t) is 

i 

y(t) = f ( I hx(t-t) u(t) dr ) = f (h x * u) 

0 

where h x (t) = L ^ H x (s) and f(x) = x • (1 + |x|). 

Using this 5 -parameter model we can demonstrate the effects 
of different 

-gust inputs ( random gust , Worst Case Gust, 

Spectral Gust), each with varying intensities 

-design load definitions: 

ymax , Oy , II y || 2 , N(y) = N des , and 
P(q > y) = P d „ [3]. 

-nonlinearities: f(x) = x • (1 + |x|)* 1 , 

-system damping! p 0 = 0.10, 0.02 


Computation of RMS based design load 
(assuming Gaussian distribution) 

The RMS —definition of design load as defined in JAR/ACJ 
25.305(d) and based on FAA-ADS-53 , is given by 

,, t nr- _ 2l - RMS W 

ydes,RMS U ° % * % O w RMS(u) ’ 

where the design gust velocity U 0 (= o w q d ) is prescribed as a 
function of altitude, but the gust velocity o w itself is not 
prescribed. But for nonlinear systems the A value depends on 
the gust velocity o w . 

For our simple nonlinear model we can compute the 
RMS-based design load analytically, assuming Gaussian 
probability distribution of the gust velocity u(t): 

p - < “ ) " M&je; 

Then the linear response x = II X u has Gaussian distribution 
p x (x), with RMS o x = o w || H X U l^ . Now the nonlinear 
response y(t) = f(x(t)) has a distorted Gaussian distribution: 

Py(y) = P*W / ^ > where y = f ( x )> 

and f(x) is strictly increasing. We obtain 


oo 

Oy = J y 2 Py(y) dy = J IW P*(x) dx = 

— oo — oo 

= o x [l + /32/jt o x + 3o x ]. 

Here we have used the formulae 
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00 

J X k e _Xx2 dx = Ax °2 1 r (-- t 1 ), X > 0,k > - 1 and 
o 

T(x + 1) = xT(x), T(l/2) = it 1 / 2 . Hence the RMS—design 
load is 

y<Jes,RMS “ U„ IIH X G 1 [l + JW* o x + 3o 2 ] 1/2 



which depends on gust intensity and nonlinearity. 


Computation of MFB design load 

The MFB design load for a nonlinear system H y can be written 
as 

y des.MFB = max max, |y(t)| , 

where y = y [u] = y [G w ]. For our model y=f(H*u). Now 
since f(x) is a monotonic function, the worst case gust fory(t) is 
the same as for x(t), hence 

ydes,MFB ^ (Xdes) X des (^ “h X des ), where 

X des x deSj MFB U„ JJ H X G ||(,;. 


Computation of MFB —IDS design load 

This simplified MFB—method was recently published by NASA 
including software [8]. Here the search for worst case gust is 
restricted to a 1—parameter family. The design loads obtained 
numerically will be slightly smaller than from the full MFB 
method. 


Computation of Spectral Gust based design loads 

The discrete gust used as input is a scaled Spectral Gust 

U a u G (t) = U 0 /iff L -1 G(s) 

(seethe analytical solutions in section 5). The corresponding 
Spectral Gust design load of our nonlinear system is 

y d es,SG =|| f (hx * (U 0 U G ) II, / Jr, (T = L/V). 


Computation of exceedance rate based 
design load (for Gaussian distribution). 

Based on timeplane stochastic simulation of the nonlinar aircraft 
with input process u(t) assumed to be stationary Gaussian with 
RMS(u) = o g (= U a /3 e.g.) given and output process 
y(t), the exceedance rate based design load y des EXC is defined 
by the equation 

N y (ydes,Exc) = N des (= l/50000h e.g.), 
where N des is a prescribed design frequency, and the function 
N y (y,) denotes the exceedance rate of level y, for signal y(t) 
(number of crossings per unit time of level y, with positive 
slope). 

For the linear subsystem x = H x (u) the exceedance rate is given 
by Rice’s formula [1] 

N x (x) = N 0x exp(-^£), where 

N = ± °* II ia)H * G II 

0jX 2jc o x 2k || H X G || 

is the ’’characteristic system frequency” for given PSD 
= | G(ico) | 2 ; solving for x we get 



For our nonlinear response, y(t) = f(x(t)), and from 
monotonicity of f(x) each x—crossing of level xl corresponds to 
an y—crossing of level yl = f(xl): 

Ny(yi) = N x (xj), where y, = f(x,) . 

We obtain 


where x desRMS „ x , ia 

f \V /2 

_ X des.EXC 1 / N 0,x\ 

’Iexc/rms - x des _ RMS “ r) d 

is a design load ratio for linear systems and Gaussian processes. 


Note that N 0 x is infinite if | H x (iw) | co 2 / 3 is unbounded, since 

for the von Karman Spectrum |G(ico)| - to -5/6 ; hence in our 
model x des EXC = oo 5 if a 0 ^ 0 . 


Graphical presentation of design loads from variousmethods 


This can be computed only numerically. The intensity U„ and 
default system parameters 
P 0 = 0.10, (i) 0 = 4 ir, a 0 = c 0 = 0 , b 0 = 1, 
were varied. 



Time [sec] 



Time [sec] 

Time response due lo Spectral Gust ul(t) 


Cont. Turb. Loads for Nonlinear System y=f(H*u) : Open Loop 



gust velocity: k * Usigma [ft/sec] 
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7. CONCLUSION 

Based on discrete gusts having the ’correct’ power spectral 
density,the Spectral Gust method represents a direct translation 
of the PSD method into time domain. 

-It requires no optimization (which is inherent 

in Worst Case Gust methods by definition) , 

-The gust shape is independent of the load stations and 

their response characteristics. 

....’’the airplane does not generate the turbulence”.. [1] 

—It produces loads consistent to DEA for linear A/C 

-It gives also correlated (balanced) loads 

-It is straight forward and very economical 

The problem that several Spectral Gusts exist, leading to slightly 
different loads in case of nonlinear systems, reflects the lack of 
agreed definition of ’design load’. The selection of a Spectral 
Gust and a reference time T (varying the scale of turbulence) 
could be adopted, e.g. to the RMS definition, which was used in 
DEA . 

The final Spectral Gust parameters are then obtained from 
general studies during the means—of—compliance process. 
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A STUDY OF THE EFFECT OF STORE UNSTEADY 
AERODYNAMICS ON GUST AND TURBULENCE LOADS. 
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H. Climent 
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CONSTRUCCIONES AERONAUTICAS S.A. 
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SUMMARY 

Gust and turbulence loads are the most severe con¬ 
ditions for some parts of the aircraft, specially in 
the outboard wing region. The addition of external 
stores usually alleviates the gust/turbulence loads on 
the wing but becomes critical for other components 
as pylons and wing/pylon attachments. 

Safety implications and all-weather envelope of 
current fighters force an accurate study of the 
gust/turbulence response of the aircraft for each 
store configuration. Inertia is the most important 
contribution of the store to the gust/turbulence 
loads. Nevertheless, store unsteady aerodynamics 
should be included if significantly changes the results 
of the gust/turbulence response. 

This paper is devoted to the analytical study 
of the effect of store unsteady aerodynamics on 
gust/turbulence response. A large number of con¬ 
figurations have been included in the study. Several 
regulations to define the gust/turbulence excitation 
have been used. 

1 INTRODUCTION 

The effect of store aerodynamics in aeroelastic prob¬ 
lems has been devoted mainly to study its influ¬ 
ence in wing/store flutter. Published investigations 
started as early as 1949 [1], and continue over the 
time. A complete multivariate analysis can be found 
in [2], 

On the other hand, very little published work con¬ 
cerning the effect of store aerodynamics in gust and 
turbulence loads is available. The reason may be 
in the fact that the addition of underwing external 
stores usually alleviates the gust/turbulence loads on 
the wing, while the flutter speed of the aircraft may 
be adversely affected because of the inertia, elastic 
and aerodynamic coupling between the wing and its 
stores. 

Nevertheless, some parts of the airframe as wing- 
pylons attachments or the pylons themselves become 
critical with the addition of external stores. If the 


load envelope due to gust and turbulence overshoots 
the rest of loads envelope (manoeuvre, buffet, re¬ 
lease...) as is the case using some current regulations, 
some pylon sizing cases will be product of gust and 
turbulence analyses. Therefore the accurate compu¬ 
tation of the gust and turbulence responses is also a 
matter of safety in an all-weather flight envelope. 

In the subsonic regime, current analysis standards 
[3]-[4] make possible to account for the wing/store 
unsteady aerodynamic coupling in gust and turbu¬ 
lence response while in the supersonic regime it is 
still an open question. The computation of the un¬ 
steady aerodynamic contribution of the store and 
the store/wing interference effect will increase signif¬ 
icantly the cost of each single run and the associated 
engineer cost of data preparation, model validation 
and results evaluation. The economics of including 
store aerodynamics must be considered together with 
other technical aspects. 

Due to the variety of stores that can be carried by 
a typical modern fighter and their derived configura¬ 
tions considering the release sequence (including fail¬ 
ure cases), the total number of possible aircraft/store 
configurations is in the order of thousands. This con¬ 
stitutes a large amount of work for flutter analyses, 
because even using engineering judgement criteria 
several hundreds should be selected. In gust and tur¬ 
bulence response inertia is the most important con¬ 
tribution of the store to the pylon loads and therefore 
a significant reduction in the number of configura¬ 
tions to analyse can be made based on that param¬ 
eter. It means that the inclusion of store unsteady 
aerodynamics in gust and turbulence loads calcula¬ 
tion is feasible with current computers, and should 
be included if significantly changes the results of the 
gust and turbulence response, specially when changes 
are nonconservative. 

Multivariate analysis of similar studies [2] has con¬ 
cluded that general guidelines can not be provided, 
only particular advices for use in a particular aircraft. 
Therefore, it has been preferred to deal directly with 
an actual set of external stores in an actual aircraft 


Paper presented at a Workshop on “Aircraft Loads due to Turbulence and their Impact on 
Design and Certification" held in Lillehammer, Norway on 5 May 1994. 
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model. 

The paper presents the general procedure followed to 
calculate the gust and turbulence load envelopes at 
monitoring stations. Some validations and sensitiv¬ 
ity analysis of store aerodynamics are briefly out¬ 
lined. The main part of the paper is devoted to 
present the effect of the store unsteady aerodynam¬ 
ics on gust/turbulence response by covering a large 
number of actual wing/store configurations first in 
single carriage and then in multiple carriage. The 
paper ends with conclusions and guidelines derived 
from the results of the analysis and description of 
some limitations of current methods. 

2 GENERAL PROCEDURE 

Loads acting on the aircraft due to the atmospheric 
turbulence are computed as response either to a 
tuned discrete gust or to continuous turbulence. A 
complete survey of vertical and lateral, “up” and 
“down”, gust conditions, flight envelope points and 
mass configurations covering all possible airplane 
characteristics and current regulation requirements 
[9]-[12] have to be performed in the certification pro¬ 
cess. Nevertheless, for the intention of the present 
study it has been considered enough to analyse only 
tuned discrete gust and continuous turbulence ver¬ 
tical excitation (up and down) in one point in the 
subsonic regime. 

2.1 Tuned discrete gust response 
The tuned discrete gust definition used for the anal¬ 
ysis corresponds to the requirements of the DEF 
STAN 00-970 regulations [9] which define the gust 
intensity without attenuation with gust length. Fig¬ 
ure 1 shows the shape of the five different gust lengths 
used in the analyses, from 100.0 ft to 473.5 ft (25 c, 
c = mean aerodynamic chord) with 66 f.p.s. of gust 
intensity. The loads obtained with this excitation are 
LIMIT loads. 



0.1 0 2 0 1 0 4 0 * O • T »•) 


Linear behaviour of the structure is assumed. There¬ 
fore, the total loads are obtained adding up the tran¬ 
sient loads due to the gust or turbulence response to 
the stationary loads due to the flight condition anal¬ 
ysed. The flight condition considered in the results 
presented herein corresponds to the 1G level flight. 
MSC/NASTRAN aeroelastic capability [4] running 
in a CRAY XM-P computer and some internal CASA 
in-house codes have been used to obtain the tran¬ 
sient response to the gust or turbulence excitation. 
A modal solution has been selected to reduce the 
computer time needed for the analysis. The dynamic 
equation defining the gust problem in the frequency 
domain is: 


{-w 2 [M} + iw[B] + [K] + [<?]}{£} - {P(o;)} (1) 


where: 


UJ 

Frequency 

[M] 

Generalized mass matrix 

[B] 

Generalized damping matrix 

[K] 

Generalized stiffness matrix 

[Q\ 

Generalized aerodynamic coefficients 

(O 

Generalized modal coordinates 


Generalized aero gust loads 


The elements of [M], [B] and [Tf] are obtained from 
the structural model. The unsteady aerodynamic 
theories (Doublet Lattice and Slender and Interfer¬ 
ence bodies method [5]-[8]) give the [Q] matrix, as 
a function of Mach number and reduced frequency, 
in the frequency domain. A Fourier transform of the 
time histories defining the gust excitation is needed 
to obtain the aerodynamic loads {P(w)}. (£(u>)} is 
obtained as the solution of the system, and an Inverse 
Fourier Transform provides {£(<)}. 

Back substitution gives the nodal displacement vec¬ 
tor as: 

= mm ( 2 ) 

where [<f>\ is the modal matrix. 

From this displacement vector, a quasi-static set of 
nodal gust loads can be obtained, at each time in¬ 
stant t, using the stiffness matrix [Kaa]'- 

= [iYHWO} (3) 

This nodal loads distribution is summed up for each 
previously selected monitoring station (j) to obtain 
the resulting loads. Total loads are obtained sum¬ 
ming up positive and negative incremental gust loads 
(to account for up and down gusts) to the stationary 
loads (1G loads). 

F T] (t) = F Sj ± F gi (t) 


Figure 1: Tuned discrete (1 — cos) gust 


(4) 
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Single gust load envelopes at each monitoring station 
are obtained from these time-histories for each gust 
length and aircraft configuration as shown in figure 2. 



Figure 2: Typical single gust load envelope 


2.2 Continuous turbulence response 
For a fighter aircraft, current military regulations 
[10]-[11] require continuous turbulence analyses in 
the supersonic regime. This study has been focused 
in the subsonic regime, therefore the definition of the 
continuous turbulence model has been obtained fol¬ 
lowing a civil regulation [12]. The continuous tur¬ 
bulence model used corresponds to the Von Karman 
Spectrum, with scale of turbulence L=2500 ft and 
112.5 f.p.s. of gust intensity. Loads obtained with 
this excitation are LIMIT loads. 

The airplane response to continuous turbulence is 
calculated by means of typical statistical approach of 
random analyses [13]. The unit Power Spectral Den¬ 
sity (PSD) of airplane response (Sj (w)) can be com¬ 
puted from the unit PSD of the atmospheric turbu¬ 
lence (S„(w)) using techniques of frequency response: 

Sj(u) = \H ja (u>)\ 2 x S a (u>) (5) 

where Hj a (w) is the frequency response of any phys¬ 
ical variable Uj , due to an excitation source, Q a that 
verifies the transfer function theorem: 


Global gust load envelopes at each monitoring sta¬ 
tion are obtained ’’enveloping” single load envelopes 
for all gust lengths and airplane configurations. 
These gust loads are then compared with other en¬ 
velopes (i.e. static loads) in order to select points to 
provide sets of nodal loads for checkstress. Figure 3 
shows a typical pylon global gust load envelope com¬ 
pared with a a typical static load envelope. The gust 
envelope overshoots the static load envelope at some 
points. In this example, eight sets of nodal loads 
were provided. 
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Uj(u) - Hja(u) X Qa{w) (6) 

where Uj(u ) and Q a (ui) are the Fourier transforms 
of uj and Q a . The rms value ( A ) of each response 
is computed as the squared root of the PSD area of 
the response: 



The limit loads are obtained as: 


L — u a x A (8) 

Where L is the limit load and u a the gust inten¬ 
sity (112.5 fps), value that corresponds to the r.m.s. 
value of the gust velocity times a factor to account 
for peaks. In-house CASA programs are used to com¬ 
pute the rms values of the internal loads (shear force, 
bending moment and torsion moment) at the selected 
monitoring stations. 


Figure 3: Typical pylon global gust load envelope 
with points for nodal loads delivery 
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3 ANALYSIS CONDITIONS 

All the results presented in this paper have been com¬ 
puted in the subsonic regime at Mach=0.8, sea level 
conditions and considering only vertical excitation. 

3.1 Structural Models 

The dynamic model used corresponds to an actual 
fighter aircraft model, with delta wing, foreplane, 
wing tip pod and four underwing stations to carry 
external stores. The model has been obtained by the 
assembly of several substructure models representing 
the main components of the aircraft. Normal mode 
frequencies and mode shapes have been validated us¬ 
ing ground resonance test results. 

The pylon models consist of a mass joined to the wing 
structure through a stiffness matrix with a reduced 
set of degrees of freedom, because the dynamic be¬ 
haviour of the store is adequately represented by a 
few modes (pitch, yaw, lateral, roll etc.), that have 
also been validated with ground resonance test. 

3.2 Aircraft Unsteady Aerodynamic Model 
Unsteady aerodynamics has been calculated using 
Doublet Lattice method for lifting surfaces combined 
with Slender and Interference Body Methods for 
fuselage and tip pod [4]-[8j. 

Figure 4 gives a representation of the aerodynamic 
model. The wing aerodynamic model consists of a 
flat surface with 182 boxes and a body for the tip 
pod with 6 slender and interference elements. The 
number of boxes for the foreplane is 30. The fuselage 
model consists of 9 slender and interference elements. 
Values of C[ ja and position of the wing pressure cen¬ 
tre obtained for zero reduced frequency have been 
compared with the values obtained using more ac¬ 
curate steady aerodynamics methods and available 
wind tunnel results showing good correlation. 


Four pylons located in four different spanwise posi¬ 
tions of the wing structure have been considered: 


Pylon 

ID. 

Spanwise 
position (%) 

Inboard 

(I/B) 

32.0 

Center 

(CTR) 

53.0 

Outboard 

(O/B) 

71.0 

Stub 

(STUB) 

87.0 


For each pylon, several monitoring stations are de¬ 
fined, representing the loads acting on the store, store 
-f launcher or store + launcher + pylon. The com¬ 
puted loads at each monitoring station are: 

• F y (Lateral Force) 

• F z (Vertical Force) 

• M x (Bending Moment) 

• M y (Pitch Moment) 

• M . (Yaw Moment) 

Figure 5 shov. s a pylon monitoring station. 
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3.3 Monitoring stations definition 

Monitoring stations have been defined covering all Figure 5: Pylon monitoring station, 

aircraft components. In gust/turbulence response, 
wing and pylon monitoring stations have received a 



Figure 4: Aircraft Unsteady Aerodynamic Model 
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3.4 Sensitivity analysis to unsteady aerody¬ 
namic modelling of store and pylons 
The loads convergence with unsteady aerodynamics 
model mesh has been verified with sensitivity analy¬ 
ses to number of: 


Table 1 describes the different stores considered in 
the analysis, the pylon or pylons where it has been 
assumed that these stores can be carried, weight and 
maximum cross section. 


• Slender body elements. 

• Store fin boxes. 

• Store canard boxes. 

• Pylon boxes. 

A quick convergence of loads has been observed with 
the number of slender body elements. Using as refer¬ 
ence 25 slender body elements, 3 used to be enough 
in most of the cases to reproduce 98 % of body loads. 
In addition, no variation of loads has been observed 
above 10 slender body elements. For lifting surfaces, 
results have shown that relatively coarse meshes can 
be used for I/B, CTR and O/B pylon modelling while 
finer meshes should be used for STUB pylon, store 
fins and store canards. 

3.5 External Stores Considered 
The different type of external stores that can be car¬ 
ried in a fighter aircraft can be grouped, from the 
gust loads point of view, in three classes: fuel tanks, 
missiles and bombs. In general, heavy stores corre¬ 
spond to full fuel tanks while light stores are missiles. 

The aerodynamic shape of the external fuel tanks 
corresponds to a slender body incorporating in some 
cases stabilizing fins. Missiles include both types of 
lifting surfaces: canards and fins to allow missile con¬ 
trol. The third type of stores (bombs, equipment 
pods etc.) corresponds, in general, to stores with 
slender body without fins or control surfaces, and 
masses and cross-sectional areas that are between 
those of fuel tanks and missiles. 

For this study, two types of external fuel tanks (one 
including stabilizing fins) have been considered. For 
missiles, a typical short range air-to-air missile has 
been considered in single or twin carrier configura¬ 
tion. For bombs, two typical heavy bombs have been 
analysed. The second bomb has been included in the 
study to allow the analysis of a twin bomb carrier. 

Some functional assumptions have been performed 
to allocate the stores on the pylons. Usually, fuel 
tanks are carried at a unique underwing position with 
the adequate fuel system connections to allow fuel 
pumping from the fuel tanks. Missiles are located at 
the outboard regions of the wings. Bombs and pods 
can be carried at several span wise positions, but in 
general, the heavier the store the more inboard is 
carried. 


TYPE 

ID. 

UNDER WING STORES 

w 

SECT. 



I/B 

CTR 

m 

o 

STUB 

(kg) 

(m 2 ) 

FUEL 

T1 


X 



1300. 

0.430 

TANK 

T2 


X 



920. 

0.283 

MIS- 

M 




X 

157. 

0.014 

SILE. 

2M 



X 

X 

297. 

0.028 


B1 

X 

X 

X 


496. 

0.139 

BOMB 

B2 



X 


263. 

0.138 


2B2 

X 




591. 

0.276 


Table 1: External stores considered 


Five of the seven stores are always carried in the same 
pylon. The twin missile carrier can be carried in two 
pylons and the B1 bomb in three pylons. This gives 
a total of 10 “single store” configurations. In addi¬ 
tion large number of “multiple store” configurations 
can be arranged combining simultaneously the single 
ones. 

3.6 Unsteady Aerodynamic Models for External 
Stores 

Unsteady aerodynamics has been calculated using 
Doublet Lattice method for the pylon surfaces, store 
fins and canards, combined with Slender and Inter¬ 
ference Body Methods (Images) for the underwing 
stores [4]-[8]. Full interference effects with the air¬ 
craft unsteady aerodynamic model have been taken 
into account. 

The final unsteady aerodynamic models of pylons 
and stores have been obtained after the sensitivity 
analyses representing a: balanced trade-off between 
economy and accuracy. Pylon lifting surfaces models 
are shown in figure 6 and store models in figure 7. 
The store models are not at the same scale to allow 
visibility of small stores details. Only slender models 
have been represented. 

NOTE: current methods do not allow easy modelling 
of stores for which control surfaces are inside of in¬ 
terference body. Interference body section should 
match slender body maximum section, which in turn, 
should match store maximum section. Control sur¬ 
faces close to longitudinal axis should be modelled 
translated outside the interference body with the as¬ 
sociated problem of structural-aerodynamic interpo¬ 
lation. Figure 8 shows two typical stores that present 
this modelling problem, a laser guided bomb with ca¬ 
nards and a balistic bomb with fins. 



I/B PYLON 


O/B PYLON 



Figure 6: Pylons unsteady aerodynamic models 



FUEL TANK (Tl) 



BOMB (Bl) 





FUEL TANK (T2) 


BOMB (B2) 



TWIN MISSILE CARRIER (2M) 


Figure 7: Stores unsteady aerodynamic models (not to scale) 
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Figure 8: Stores difficult to model with current methods 
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4 RESULTS OF “SINGLE STORE” CONFIG¬ 
URATIONS 

The 10 single store configurations have been consid¬ 
ered in a first step to present results in a system¬ 
atic way. Although the loads have been obtained 
for all pylon and aircraft monitoring stations, results 
are presented only for the pylon monitoring station 
that takes into account the contribution of store plus 
launcher. Loads have been obtained with and with¬ 
out store aerodynamics. 

The loads envelopes displayed in figure 9 show a typ¬ 
ical effect of the store unsteady aerodynamics on py¬ 
lon loads. Continuous lines correspond to the case 
without store unsteady aerodynamics while dotted 
lines correspond to the case with aerodynamics. The 
larger envelopes correspond to all gust lengths of the 
tuned discrete gust. The small envelopes correspond 
to only considering the 25c gust length case. 

Some regulations [11]-[12] request for compliance 
with 25c gust length only. Figure 9 shows that due 
to the small loads obtained in this particular case, to 
consider or not the store unsteady aerodynamics is a 
matter of minor importance. 


-WITHOUT STORE UNST. AERO. 

-WITH STORE UNST. AERO. 



Figure 9: Effect of store unsteady aerodynamics on 
gust load envelopes 

Figure 10 shows the effect of the store unsteady aero¬ 
dynamics in the response to continuous turbulence. 
Continuous line is the case without store unsteady 
aerodynamics and dotted line is the case with aero¬ 
dynamics. 


_ WITHOUT STORE UNST. BEROOYNAM1CS 

.WITH STORE UNST. AERODYNAMICS 




Figure 10: Effect of store unsteady aerodynamics on 
continuous turbulence response 


From the envelopes (gust) or limit loads (turbu¬ 
lence), the effect of the unsteady aerodynamic mod¬ 
elling has been computed for each store-pylon com¬ 
bination and each load:. 


AiOA£>(%) = 100 x 


LOAD A - LOADq 
LOAD 0 



Where: 

LOAD Peak values (gust) or 

Limit loads (eq. 8) (turbulence) 

A ’’With Aerodynamics” 

0 ’’Without Aerodynamics” 

Table 2 displays the values obtained for the 10 sin¬ 
gle store configurations. The two figures in each box 
correspond to the influence in the tuned discrete gust 
response (above) and the influence in the continuous 
turbulence response (below). Negative values rep¬ 
resent alleviating effect of considering aerodynamics 
while positive values represent an increase of loads. 
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TYPE OF STORES 

| FUEL TANKS 

MISSILES 

| BOMBS j 

Tl 

CTR 

T2 

CTR 

M 

STUB 

2M 

1 _B1 

B2 

O/B 

2B2 

I/B 

O/B 

STUB 

I/B 

CTR 

O/B 

A F y (%) 

4.4 

4.6 

-8.8 

45.3 

23.8 

11.1 

1.7 

-8.5 

-0.6 

43.0 


0.6 

2.6 

-10.0 

15.3 

17.1 

10.2 

1.6 

-10.1 

-5.3 

37.1 

A F t (%) 

-7.9 

-1.5 

7.2 

16.6 

58.7 

-1.1 

0.9 

-2.3 

-1.3 

-7.3 


-8.1 

-3.9 

2.9 

5.9 

35.5 

-1.3 

0.5 

-2.3 

-1.6 

-7.1 

A M x (%) 

4.0 

4.8 

6.1 

60.2 

17.0 

9.6 

1.7 

-8.8 

-1.0 

38.9 


0.5 

2.7 

-8.4 

20.2 

-3.5 

8.3 

1.3 

-10.1 

-5.4 

30.2 

A M y (%) 

-17.3 

-5.8 

-5.0 

-2.4 

-20.9 

-6.3 

-6.9 

-5.2 

-6.0 

-7.3 


-24.7 

-4.2 

-7.1 

-7.9 

-9.3 

-1.3 

-6.2 

-3.3 

-4.2 

-6.8 

A M z (%) 

17.0 

6.6 

-6.4 

159.1 

-21.7 

6.5 

30.9 

27.1 

2.8 

7.0 


17.4 

2.8 

-13.8 

44.6 

-15.4 

1.3 

29.1 

28.5 

15.0 

1.9 


Table 2: Variation of gust/turbulence loads due to store unsteady aerodynamics for “single store” 

configurations 


In general, results show that similar effects are ob¬ 
tained for both excitations tuned gust and contin¬ 
uous turbulence. It is remarkable that the effect is 
always alleviating for the M y load. The largest store 
unsteady aerodynamics effects are obtained for the 
missiles and twin carriers, as was expected because 
of the presence of lifting surfaces embodied in the 
models of those stores. On the other hand, a small 
effect is reported for stores consisting basically on 
slender bodies. To get a better insight into the dif¬ 
ferent aerodynamic contributions, the following sen¬ 
sitivity analyses have been performed: 

• PRESENCE OF FINS IN FUEL TANKS. 

A tuned gust run was performed removing the 
fins in the fuel tank T2 model. Comparison of 
results with the run including the fins shows that 
the influence of fins is not significative. 

• PRESENCE OF FINS AND CANARDS IN 
THE MISSILE. 

Two tuned gust runs were performed either re¬ 
moving the fins or the canards in the Missile M. 
model. Comparison of results with the run in¬ 
cluding fins and canards show that all the effect 
is basically due to the canards and the influence 
of fins is not significative. 


• PRESENCE OF CARRIER 

Two tuned gust runs were performed removing 
the aerodynamic model of the carrier but main- 
tainig the close bodies in the twin missile 2M 
and twin store 2B2 carriers. Results show that 
although the influence of the close bodies is very 
important in both cases (bombs and missiles) 
the influence of the carrier itself should also be 
taken into account. 

• SPANWISE POSITION OF SLENDER BOD¬ 
IES 

The store B1 is carried in three spanwise posi¬ 
tions. Results from table 2 show the effect of 
store unsteady aerodynamics: 

- Increases on M z with spanwise position. 

- Decrases on F y and M x with spanwise po¬ 
sition. 

- Remains almost constant on F t and M y . 
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4.1 Effect on Wing Monitoring Stations 
Usually, inclusion of stores produces an alleviation 
on wing gust/turbulence loads. Nevertheless, an ex¬ 
ample has been selected to illustrate that depending 
on the configuration and monitoring station, a wide 
variety of cases can occur, including increase of loads 
after store inclusion. The consideration of the store 
unsteady aerodynamics can contribute to expand the 
load envelopes. 

The following plots show the comparison of results 
for three configurations: 

• Clean wing (without missile). 

• Missile M on stub pylon without store aerody¬ 
namics. 

• Missile M on stub pylon with store aerodynam¬ 
ics. 

Figure 11 shows the comparison for a typical tuned 
gust load envelope in the wing root monitoring sta¬ 
tion, far inboard the missile section. The influence 
of the missile inertia is small and the effect of the 
unsteady missile aerodynamics very reduced. 

Figure 12 shows the comparison in the wing mon¬ 
itoring station located at the same missile position. 
The presence of the missile inertia substantially mod¬ 
ifies the load envelope shape, producing additional 
critical points. The unsteady missile aerodynamics 
significantly contributes to enlarge the load envelope 
and therefore, if not included in the calculations, non 
conservative results can be obtained. 


WITHOUT MISSILE 


WITH MISSILE 



•WITHOUT MISSILE 


■WITH MISSILE 



Figure 12: Wing 89% spanwise monitoring station 


WITHOUT MISSILE 
WITH MISSILE 



Figure 13: Wing 99% spanwise monitoring station 

Figure 13 shows the comparison at a wing monitor 
station located outboard from the missile position. 
Comparison shows the clear alleviating effect of the 
missile inertia, only slightly modified by the missile 
unsteady aerodynamics. 


Figure 11: Wing root monitoring station 
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5 RESULTS OF “MULTIPLE STORE” CON¬ 
FIGURATIONS 

In principle, it is assumed that the inertia added in 
the “multiple stores” configurations should alleviate 
the loads in each particular pylon already obtained 
with the “single store” configurations. Although in 
general the assumption is quite appropriate, some 
critical load cases coming from multiple stores con¬ 
figurations have been obtained. 

Configurations with multiple stores can been ar¬ 
ranged combining the stores of table 1. The max¬ 
imum number of possible multiple store configura¬ 
tions including all possible combinations is in the or¬ 
der of hundreds. Nevertheless, a large reduction in 
the number of configurations can be achieved based 
on particular functional considerations (non-flight 
configurations, two or more failures in the release se¬ 
quence, etc.). An' additional reduction in the number 
of configurations to analyse can be obtained remov¬ 
ing the configurations which loads are clearly cov¬ 
ered by other configurations (lightest stores at each 
pylon). 

The remaining “multiple stores” configurations have 
been analysed with and without store unsteady aero¬ 
dynamics and the loads due to gust/turbulence re¬ 
sponse obtained. At the end, 6 “multiple stores crit¬ 
ical configurations” have been obtained. A configu¬ 
ration is considered critical when it produces critical 
loads in at least one monitoring station. Loads are 
critical when overshoot the envelopes already pro¬ 
duced by “single store” configurations. 



UNDER WING STORES 

I/B 

CTR 

O/B 

STUB 

1 

Bl 

Bl 

Bl 

M 

2 

2B2 

T1 

B2 

M 

3 


T1 

2M 

M 

4 


T1 


M 

5 

2B2 



M 

6 



Bl 

M 


Table 3: “Multiple Stores” Critical Configurations 

An attempt has been made to reduce the amount 
of computation, trying to derive guidelines on which 
store unsteady aerodynamics should be considered 
and which can be removed in a multiple store config¬ 
urations. Several runs have been performed with the 
6 critical configurations considering a partial num¬ 
ber of stores with unsteady aerodynamics instead of 
all of them. Table 4 shows a typical example of ob¬ 
tained results, corresponding to the first critical con¬ 
figuration. Again, percentages have been obtained 
following equation 9. In each box, values above rep¬ 
resent the effect in tuned discrete gust and values 
below represent the effect in continuous turbulence. 
Negative values mean alleviating effect and positive 


AERO 

A 

UNST. AERO.EFFECT j 



I/B 

CTR 

O/B 

STUB 


Fy (%) 

-7.8 

13.6 

-10.9 

-21.3 



-4.6 

6.0 

-9.9 

-10.8 


Fz (%) 

8.9 

15.9 

4.7 

1.5 



2.8 

6.4 

-0.4 

-2.7 

3 Slen. 

Mr (%) 

-9.8 

13.0 

-10.9 

25.9 

+ miss. 


-4.6 

5.6 

-10.0 

-19.6 


My (%) 

4.4 

-0.6 

-6.2 

-41.7 



0.3 

-13.7 

-3.8 

-26.6 


Mz (%) 

11.3 

103.3 

29.0 

18.8 



16.8 

21.5 

-32.3 

-26.9 


Fy (%) 

-9.8 

-3.3 

-11.1 

-2.2 



-5.7 

-7.8 

-12.1 

-4.7 


Fz (%) 

10.1 

13.7 

-2.8 

-2.4 



7.5 

5.6 

-3.7 

-3.1 

3 Slen. 

Mr {%) 

-12.2 

-3.7 

-11.2 

-1.3 



-5.5 

-8.1 

-12.3 

-2.4 


My (%) 

4.1 

-8.4 

-13.3 

-4.8 



0.4 

-8.8 

-7.0 

-3.6 


Mz (%) 

2.2 

104.9 

28.5 

-5.4 



17.6 

22.5 

-29.5 

-3.0 


Fy (%) 

1.1 

18.8 

-0.8 

-10.8 



0.6 

18.5 

1.1 

-5.6 


Fz (%) 

-1.5 

3.5 

8.1 

2.9 



-4.6 

1.5 

3.7 

0.2 


Mr (%) 

1.1 

19.2 

-0.8 

21.9 

miss. 


0.6 

18.4 

1.1 

-18.3 


My {%) 

0.3 

9.4 

8.9 

-36.1 



0.1 

-4.6 

3.6 

-24.1 


Mz (%) 

-6.5 

4.6 

2.6 

11.9 



1.2 

-1.2 

-4.8 

-24.1 


Table 4: Variation of gust/turbulence loads 

due to store unsteady aerodynamics for a “multiple 
store” critical configuration 

values mean load increases. First block of the table 
presents the influence of store unsteady aerodynam¬ 
ics considering all the stores. Second block presents 
the results considering only the aerodynamics of the 
3 bombs slender bodies. Third block presents the 
results considering only the missile aerodynamics. 

Inboard pylon results are relatively similar between 
first and second block indicating that the presence of 
missile aerodynamics at the stub pylon has a minor 
effect on the inboard pylon. Stub pylon results are 
relatively similar between first and third block indi¬ 
cating that the presence of missile aerodynamics is 
the main contribution to stub pylon loads. 

In spite of the above expected considerations, ta¬ 
ble 4 shows no other clear guidelines. Therefore, 
the conclusion (confirmed with similar results with 
other configurations) should be that the store un¬ 
steady aerodynamics should be considered for all the 
stores in critical multiple store configurations. 
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6 SUMMARY OF RESULTS 

So far, the effect of the store unsteady aerodynam¬ 
ics has been considered by comparing with the same 
situation without aerodynamics. When different con¬ 
figurations have been analysed, a range of values rep¬ 
resenting that effect are obtained, one for each con¬ 
figuration. When application of results have to be 
made to an actual aircraft to obtain the global effect 
per pylon basis, it is better to normalize the effects 
to the same value in each pylon. Figure 14 shows 
a histogram of the effects per pylon basis obtained 
with the expresion: 


&load(%) 


100 x 


[LOAD A ) - ( LOAD 0 ) 
[MAX(LOAD 0 )] 


( 10 ) 


Where: 


LOAD 

A 

0 

MAX(LOAD) 


Peak values (gust) or 

Limit loads (eq. 8) (turbulence) 

’’With Aerodynamics” 

’’Without Aerodynamics” 
Maximum Load in that pylon 


The histogram represents the range of values of the 
store unsteady aerodynamics effect for all the con¬ 
figurations analysed with single or multiple stores. 
A table like table 5 with the maximum positive val¬ 
ues (in form of factors) for each load on each pylon 
can be derived for this particular aircraft and could 
help to decide whether a new configuration should 
be analysed including store unsteady aerodynamics 
given the gust/turbulence loads without such aero¬ 
dynamics. 



I/B 

CTR 

O/B 

STUB 

Fy 

1.37 

1.07 

1.09 

1.00 

F z 

1.26 

1.17 

1.09 

1.11 

M x 

1.36 

1.06 

1.14 

1.28 

M y 

1.07 

1.00 

1.13 

1.00 

M z 

1.16 

1.21 

2.24 

1.79 


Table 5: Factors for pylon loads 


The histogram shows larger influence of store un¬ 
steady aerodynamic than expected at the begining 
of this study, in both discrete tuned gust and contin¬ 
uous turbulence excitation. In general, comparison 
between the two types of excitation shows similar ef¬ 
fects with the only exceptions of STUB pylon F z , M x 
and M z . 
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Figure 14: Histogram of Global Effects per Pylon 
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7 CONCLUDING REMARKS 

An extensive analytical study of the effect of store 
unsteady aerodynamics on gust/turbulence response 
has been presented. The influence has been larger 
than expected in both discrete tuned gust and con¬ 
tinuous turbulence. The effect in single store con¬ 
figurations is very similar for both excitations and 
also for multiple store configurations with the only 
exceptions of STUB pylon F z , M x and M,. 

In general, considering store unsteady aerodynamics 
alleviates the pylon M y loads but there are not clear 
trends for the rest of loads. 

“Single store” configurations results indicate that un¬ 
steady aerodynamics of stores with lifting surfaces, 
(missiles, twin carriers) have to be considered. Oth¬ 
erwise, non-conservative results can be obtained even 
for some wing monitoring stations close to the store 
position. For missiles models, the contribution of 
canards is the most significative with only limited 
influence of fins. 

“Multiple store” configurations results indicate that 
it is necessary to include unsteady aerodynamics of 
all stores in critical configurations. Apart from this, 
no other clear guidelines can be derived from the re¬ 
sults. 

Current methods do not allow easy modelling of 
stores which control surfaces that are inside of in¬ 
terference body. Interference body section should 
match slender body maximum section, which in turn, 
should match store maximum section. Control sur¬ 
faces close to longitudinal axis should be modelled 
translated outside the interference body with the as¬ 
sociated problem of structural-aerodynamic interpo¬ 
lation. According with results presented in this pa¬ 
per, it will not be a serious problem with balistic 
bombs with fins, but it is likely to be a major prob¬ 
lem with laser guided bombs with canards. 

Extensions of the work performed herein can be ap¬ 
plication to engine nacelles or stores in supersonic 
regime. 
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Current continuous turbulence Power Spectral Density (PSD) 
gust analysis methods are valid only for linear aircraft. In the 
past, many methods for the analysis of nonlinear aircraft to 
meet the current PSD gust criteria have been proposed. In this 
paper, three stochastic simulation based and one deterministic 
function based method are compared and evaluated for the 
compliance of nonlinear aircraft to the current continuous tur¬ 
bulence gust criteria. The aircraft configurations analyzed in 
this paper include a symmetric aircraft model coupled to a 
nonlinear gust load alleviation (GLA) system and an anti¬ 
symmetric aircraft model coupled to a nonlinear yaw damper 
system. Results from these four nonlinear methods are com¬ 
pared to linear closed and open-loop results as well as currently 
used linear approximation methods. Computing performance 
issues are addressed to provide the reader with a complete pic¬ 
ture of the trade-off between the analysis accuracy and com 
puling cost. 

Nomenclature 

Scalars 

A, Unit gust root-mean-square value of response y 

from linear PSD analysis 

A,R,C Constants used in the approximation of the cumu¬ 
lative probability distribution integral for the prob¬ 
ability of exceedence criteria 
Non-storm and storm scale parameters for proba¬ 
bility distribution of root-mcan-squarc gust velocity 
Ct.Cj Coefficients for linear regression fit 
h y (l) Impulse response of the first system in Matched 

Filler Theory 

I! y (ia>) Fourier transform of h y (l) 

K Strength of impulse function used in the nonlinear 

Matched Filter Theory 

kMFT Matched Filter Theory waveform normalization 

factor 

No Number of zero response crossings 

Y(y) Number of crossings of response level y 

Pt,Pi Proportion of flight time spent in non-storm and 
storm turbulence 

P(y>yd) Cumulative probability distribution function for y 
greater than y d 

p(y) Probability density function for y 

Te Maximum simulation time 

to Time shift for the Matched Filter Theory waveform 

U. Design gust velocity U„ = 

IVg Input gust time history 

x(t) Matched excitation waveform in Matched Filter 

Theory analysis 

X(ico) Fourier transform of x(l) 

yd Design value of response y 

y(t) Time history of response y in Matched f iller 

Theory when excited by the matched excitation 
waveform. 

Y(iw) Fourier transform of y(t) 

*? Peak to root-mean-square ratio 


a> Circular frequency (rad/scc) 

Uniform random phase angle from 0 to 2n 
4>(«>) Input turbulence power spectral density 

Root-mcan-squarc gust intensity 
iTy Rool-mcan-square value of response y 

Subscripts 

L Linear 

Nl, Nonlinear 

Superscripts 

* Complex conjugate 

Introduction 

All modern aircraft today utilize active control systems of some 
type during normal everyday operation. These include Yaw 
Dampers and Stability and Control Augmentation Systems 
(SCAS). Future aircraft will also include Maneuver Load 
Alleviation (MLA), Gust Load Alleviation (GLA), and Elastic 
Mode Suppression (EMS) systems. To properly account for die 
effects and benefits of such systems, methods which include the 
systems nonlinear behavior, such as saturation (control surface 
rate and deflection limitations) and nonlinear aerodynamics, 
must be used. 

Certifying agencies require limit load analysis of the aircraft 
responding to continuous turbulence (CT) as well as discrete 
gust events (1-3). When the current continuous turbulence cri¬ 
teria was developed, one of the main assumptions was that the 
aircraft model was linear. This is demonstrated by the defi¬ 
nition of A in the requirements as die unit gust root-mean- 
squarc (RMS) response level, as well as the specification for 
U„ as combination of the RMS design gust intensity and a 
factor which transforms the response from a RMS to peak 
value. 

On current aircraft programs, because of this linear assump¬ 
tion, gust load analyses arc performed in the frequency domain. 
Typically when design or mission Power Spectral Density 
(PSD) atid I-cosine discrete gusts are analyzed, the nonlinear 
effects of the system (such as saturation) arc side-stepped by 
reducing the gain of the control system until die saturation is 
eliminated. I his is a conservative approach which reduces the 
full benefit of the active control system. Currently also, the 
effects of nonlinear aerodynamics cannot be properly modeled 
using frequency domain approaches, hut correction factors are 
instead applied to the linear aerodynamics. 

Actual aircraft designs arc in fact nonlinear and for some of the 
new systems introduced on modern day aircraft, they are highly 
nonlinear. 'The analysis of these nonlinear systems for the cur¬ 
rent continuous turbulence criteria requires additional assump- 
dons in order to determine the limit load or response levels. 
Satisfying the requirements of these certifying agencies when 
analyzing aircraft with highly nonlinear control systems will 
include determining the effects of these nonlinearities. 


Paper presented at a Workshop on “Aircraft Loads due to Turbulence and their Impact on 
Design and Certification” held in Lillehammer, Norway on 5 May 1994. 
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Gould(4), Nohack (5), Vinnccomhc, llockcnhull and Dudman 
(f>), and Scot!, Polotzky, and Perry (7-9) have published pro¬ 
cedures to comply with the design envelope requirement for 
nonlinear aircraft in the current continuous turbulence criteria. 
Goggin (10) includes some preliminary work performed to 
correlate stochastic linear and nonlinear time domain limit load 
results with those from traditional frequency domain PSD gust 
analyses. 

This paper provides a brief overview of the procedures listed 
above, as well as a comparison of the results from these 
approaches for a large transport aircraft. Results for a vertical 
gust ease with a nonlinear active gust load alleviation system 
and a lateral gust ease with a nonlinear yaw damper control law 
will be compared and correlated with the linear results from the 
traditional frequency domain analyses. 

Mathematical Model 

Goggin (II) developed an extremely accurate state-space for¬ 
mulation of the flexible aircralt/activc control system model. 
In this formulation, the unsteady aerodynamics arc generated 
using Doublet Lattice theory and Slender Body theory with 
Interference effects developed by Giesing, Kalman and Roddett 
(12). The motion dependent aerodynamics arc then approxi¬ 
mated with a Rational Function Approximation (RFA) devel¬ 
oped by Sevart (13) and Roger (14). 

Hie unsteady gust forces arc approximated by a RFA that was 
developed by Dykman (15). This form of RFA was demon¬ 
strated to be a much more accurate description of the unsteady 
gust forces, especially at higher frequencies and for components 
that are a considerable distance from the gust reference point. 
This would be true for tail forces on large transport aircraft. 

The final equations of motion arc formulated in modal coordi¬ 
nates and include the coupling of the flexible aircraft equations 
and the control law equations of motion. Loads are determined 
using a summation of force method to combine the inertial 
forces with the aerodynamic gust and motion induced forces. 


Leak Response 
** “ RMS Response 

The peak gust velocity, (/„, for that response is derived from the 
pcak-to-RMS ratio and the simulation RMS gust intensity <r„. 

V„ = rpr,, (2) 


After processing the results from many simulations using dif¬ 
ferent RMS gust intensities, a table of peak load vs. peak gust 
intensity is established. This table then provides the design gust 
response when interpolated from using the design gust velocity. 

For the analysis shown here in this paper, a modification of the 
above method was used due to the asymmetric nature of a 
spoiler based GI.A system. When determining the peak and 
RMS load for a positive response quantity, only positive 
bending peaks were used, and only that part of the response 
time history that was positive was used when calculating the 
RMS response. Similarly, the negative peaks and negative 
parts of (lie time history were used when generating negative 
response quantities. 


Proba bilil y of Fxcccdencc Criteria 


Noback (5) developed and evaluated another method for 
determining limit loads from a nonlinear aircraft using the 
design envelope criterion. In this paper this method will be 
referred to as the Probability of F.xcccdcncc Criteria (PF.C). 
This method also requires a lime domain stochastic simulation 
of the pertinent responses to be evaluated for a number of input 
RMS gust intensities. Noback suggests the use of three gust 
intensities which may be considerably less that the APC 
method. lie based the criteria so that the probability of 
exceeding the design load in the nonlinear system would be the 
same as the probability of exceeding the design load in a linear 
system. For a linear Gaussian system, the probability density 
of the output y of that system is defined as: 


(3) 


Analysis capabilities for nonlinear systems include rale and 
deflection limitations for control surface deflection degrees oT 
freedom, as well as two types of nonlinear time marching cor¬ 
rections for modifying the linear aerodynamic theory to match 
steady experimental nonlinear aerodynamic data. For further 
information on the formulation of the math model, nonlinear 
capabilities, or correlation of this math model with linear and 
nonlinear maneuver and 1-cosine gust excitations, see Refer¬ 
ence 11. 


Average Peak Criteria 


Early analyses of nonlinear active control systems (ACS) were 
performed using what will be called the Average Peak Criteria 
(APC). This method was initially used by Gould(4) for the 
evaluation of the nonlinear characteristics of a GLA system on 
a version of the L-1011 aircraft. This procedure requires the 
time domain stochastic simulation of the output response 
quantities to be evaluated for many input root-mean square 
(RMS) gust intensities. These response quantities were then 
processed to determine the average of a given number of peaks 
and the RMS of the response. Gould showed that the results 
from using one peak in this analysis and using three peaks 
produced essentially the same results. 


When the peak response quantities arc known for all of the 
responses, the pcak-to-RMS ratio tj is calculated for each 
response. 


Where e, is the standard deviation of y. 


Substituting rr } = /ijrr,, and integrating the probability density 
function above some design load yj will result in the probability 
distribution function for a given RMS gust intensity. 



This is the probability that a given design load, y will be 
exceeded in a linear system when excited by a Gaussian input 
with a standard deviation of n w . 


Now in the design envelope criteria, only (/„ is specified and not 
n„. Noback introduces probability density of turbulence that is 
used in the mission analysis criteria and was initially developed 
by Press and Slciner(lb). 




_2_ - 1 /J(» w /a 2 > 


Where P\ and /’2 are the proportion oT (light lime spent in non- 
storm and storm turbulence and fr| and hi arc the non-storm 
and storm scale parameters in the probability distribution of 
root-inenn-squnrc gust velocity. 
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Neglecting the non-storm terms of turbulence and combining 
with Equation 4: 



This is the probability that the design load yd, for a linear system 
will be exceeded. Noback proposes that this may be approxi¬ 
mated to be within 0.6% using: 


I’(X>yd)- r 2 


g b 2 


2 

1 


+ J2*ujb 2 


V) 


For nonlinear systems, the probability distribution function, 
P(y>yd,a w ). is not known and must be determined numer¬ 
ically. 


J .oo p oo 

p(<%) p (y > Yd, °«) drT w = e(%M<z„ 

0 J 0 


(«) 


Noback suggests that for an estimate of the design load, y„ t , the 
probability distribution function P(y>y es ,,a w ) is evaluated for 3 
input gust intensities. The integrand of Equation 8 is then 
approximated using the form: 

QlP J = At~ *** " cV = P(c w ) P(y >)•„,' a w ) (9) 

and this equation is integrated analytically resulting in: 

p (y>y' !t )±^ 4 ir e - 1 * C (10> 


If this is performed for a few estimates of the design load, the 
design load can be found by interpolation using the probability 
of exceeding the design load in a linear system from Equation 
7. 


Due to the asymmetric nature of some types of aeroscrvoclaslic 
systems, the probability distribution function is evaluated 
numerically for both positive and negative estimates of the 
design response level. For more details on tins method sec 
Reference 5. 


Threshold Exceedence Count Criteria 


Vinnccomb, Hockenhull and Dudman (6) developed another 
stochastic simulation based method which utilizes normalized 
exceedence counts as the defining criteria for the limit load 
level. In this paper this method will be referred to as the 
Threshold Exceedence Count Criteria (TEC). From Rice (17) 
the exceedence curve for a Gaussian system can be defined as: 


N(y) = W 0 exp 



(ID 


At the design load level, the number of exceedences is: 

111 


Niy d ) = /V 0 exp 


2(7 ^ 


( 12 ) 


Vinnecombe et. al. propose that for shorter simulations, the 
oulput exceedence curves will be biased on the statistics of the 
input gust time history and that the biasing of a linear reference 


model and the nonlinear model will be similarly distorted. An 
estimate of the limit load is (hen dclcrmincd for these short 
simulations as: 


^N,(yd NL ) Nd* n v c ) 


(13) 


Note that if the ouiput response remains Gaussian, then this 
remains consistent with the previous results derived by Rice in 
the limit. 



In this analysis, a linear regression fit was used to smooth the 
data between linear and nonlinear reference models. Al 
common normalized exceedence levels, the linear load from the 
linear reference model and the nonlinear load from the non¬ 
linear model arc plotted againsL each other. This data is then 
approximated using linear regression. Vinnecombe et. al. found 
for the type of notdincaritics under evaluation, that the linear 
regression fit accurately represented the data. As will be shown 
later in this paper, the linear regression accurately represents 
the type of nonlinearity used here as well. The form of the fit 
is: 

Yw, = ClYl. + c 2 ( 15 ) 

Where yni. is the nonlinear load at a given normalized 
exceedence level and yi, is the linear load at the same normal¬ 
ized exceedence level. Vinnecombe et. al. generally used a 
pcak-to-RMS ratio of 3, but also investigated the effect of other 
ratios. For more details on this method see Reference 6. 

This method might have some advantages over the previous 
two described. First, if the results from various input gust 
intensities arc very similar, then the nonlinear simulation may 
need only be performed for one input gust intensity on both the 
nonlinear and linear reference model. Secondly, if the input 
gust time history lends to bias the linear reference and nonlinear 
model similarly, especially near the limit load, then the simu¬ 
lation lime might be able to be reduced. 

Matched Filter Theory Criteria 


Matched Filler Theory (MFT) has been under development 
during the last decade at NASA Langley for computing maxi¬ 
mized and correlated gust loads and responses for both linear 
and nonlinear aircraft systems (7 9,18-19). For a linear aircraft 
model, the process behind MFT is shown in Figure 1 and is as 
follows: 

A unit impulse function is used as input to a gust filter modeled 
after the atmospheric turbulence spectrum. 'Flic gust filter used 
for the results within this paper is documented in Reference 20. 
This is then used as the excitation of the linear aircraft model. 
The output will be called the impulse response and is designated 
by h y (t). Flic combination of the gust filler and the aircraft 
model here will be referred to as the first system. 

Flic impulse response h y {() is then normalized by a factor k*frr, 
reversed, and shifted ill time an amount i 0 . This scaled, shifted, 
and reversed impulse response function x(J) is referred to as the 
matched excitation waveform and is then used as input to the 
gust filter/aircraft model a second time. The combination of the 
gust filler and (he aircraft model here will be referred to as the 
second system. 
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The output from Ihc second syslcm is then dcsignnlcd hy y(i) 
Note that for a linear syslcm, the lime response y(i) will reach 
a maximum at exactly t= t 0 This is because for a linear system: 




Y(i<o)e iuit doo 


X(io >): 


-iVl. 
//>«) , 


K MTT 

Now because this is an even function: 
1 


X(la,)n y (io,)e"°'d(o (17) 

(18) 
(19) 
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even function: 
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MFT 


( 20 ) 


So the maximum occurs at time f=i 0 . If k S m is selected as 
unity then the maximum ofy(t) will occur at t 0 and will be equal 
to Aj/i r. The selection of fc Mr7 - for the linear analysis is not very- 
important, because for a linear system, the response y(t) will be 
linearly proportional to the input x(i). For the nonlinear analysis 
however, the selection of kMrr is crucial. 


For the MFT analysis, k M rr is selected so as to create a 
waveform x(l ) that will be of unit energy. Jones (21) shows that 
this constraint on the input waveform enforces that all or the 
input waveforms will have an equal probability of occurrence. 
If k M rr is selected such that: 

X 

*Mrr="2~|j H y (iw)I!*(iai)doX (21) 

then for a linear system y m , x = >>(fo) = A y 

For nonlinear systems, Scott et. al. (7-9) has been actively 
researching methods to optimize the MFT waveform to obtain 
peak load levels. Two types of nonlinear MFT processes arc 
currently documented in this research. The first type oT opti¬ 
mization is what is referred to as Ihc multi dimensional search. 
This procedure uses Chebyshev polynomials for the approxi¬ 
mation of the input MFT waveform. The coefficients or the 
polynomials arc then used as the design variables in the opti¬ 
mization process with the energy of the MFT waveform con¬ 
strained to unity. This nonlinear analysis is computationally 
expensive and will not be addressed further in this paper but the 
reader is encouraged to look into the references for more details 
of this procedure. 

The second type of optimization is actually more of a trial and 
error search technique. This technique is referred to as the 
one-dimensional search method in the references and will be the 
only nonlinear MFT method documented in this paper. The 
one-dimensional search method is depicted in Figure 2 and is 
almost identical to the linear MFT process. The main differ 
cnees are: 1) the strength of the impulse response in the first 
syslcm is no longer unity but a variable (K), 2) a loop around 
the entire linear MIT analysis is created, and 3) in the the first 
and second systems for this analysis, the aircraff model is non¬ 
linear. 


If the nonlinear Mi l analysis is executed for many impulse 
strengths (AT), the resulting responses from each of these analyses 
can be searched through to generate a table of the maximum 
and minimum responses. Note that for a nonlinear analysis the 
maximum load does not necessarily occur at 


Scott et. al. showed that the results from the one-dimensional 
search closely approach those of (he multi-dimensional opti¬ 
mized results at a fraction of the computing lime (7). lie 
recommends the use of Ihc one dimensional search when using 
the MFT nonlinear analysis method. 

Aircraft Model 


l lie model used in this analysis is that of a wide body transport 
aircraft. This aircraft has 3 engines including one that is 
mounted in Ihc vertical tail of the aircraft. The aircraft also has 
both upper and lower winglcls. Seventeen elastic and three 
rigid body modes arc used in the gust analyses as well as one 
rigid mode for each of the control surfaces that are utilized in 
the active control system. The vibration modes used in the 
analysis arc generated using a beam stick representation of the 
aircraft. I hc Doublet l.atticc aerodynamic model of the air¬ 
craft is shown in Figure 3. All significant components of the 
aircraft arc included. Wing engines and the tail engine arc 
modeled using 'ring wings' and the fuselage is modeled with 
Slender Body theory. I hc methods of Images and Interference 
Bodies arc included to account for the wing-body interference. 
The steady and unsteady aerodynamics arc modified to match 
the experimental values using weighting factors that remain 
constant with frequency. 

The gust load alleviation (Gl.A) syslcm that is included in the 
vertical gust analysis utilizes sensors al the center of gravity 
(CG) of the aircraft as well as accelerometers at the winglips. 
Hie control law drives Ihc outboard ailerons and spoilers. A 
pitch compensation loop also drives the inboard and outboard 
elevator to cancel Ihc pitching forces that arc generated by the 
wing control surfaces. T he control law model used in this 
analysis has not been optimized to achieve the greatest 
reduction in loads but is included instead to provide a com¬ 
parison of the nonlinear gust analysis methods. 

Nonlincaritics for the Gl.A system include aileron and spoiler 
rate and deflection limitations. The aileron in this Gl.A anal¬ 
ysis is constrained to deflect between + 20 degrees while the 
spoiler control surfaces arc limited to 0.0 degrees trailing edge 
down (TFD) and 30.0 degrees trailing edge up (TKU). The rale 
limitations included in these analyses are slightly more compli¬ 
cated and arc not constant. Ihc rate limitations arc a function 
of the load on the actuator or the control surface. For the 
aileron, Ihc no-load rate is + 63.5 dcg'scc and the half load rate 
is + 40.0 dcg 'scc. I hc spoiler no-load rate is + 95.5 dcg/scc 
and the half load rate is + 66.4 dcg/scc. All of the deflection 
and rate limitations arc maintained at the limit load level. 

'Ihc nonlinear yaw damper model that is included in the lateral 
gust analysis utilizes roll attitude, lateral acceleration and yaw 
rale sensors tip near Ihc pilot station of the aircraft. This con¬ 
trol law model drives Ihc upper rudder control surface only and 
is typical of a failed yaw damper configuration. The analysis 
includes rate and deflection limitations similar to those used in 
the nonlinear Gl.A analysis. I hc rudder deflection was limited 
to + 2.5 degrees in this analysis and the rudder rate was limited 
to + 44.5 dcg/scc at no load and + 30.0 dcg/scc at half-load. 
As in the symmetric gust analysis, all of the deflection and rate 
limitations arc maintained at the limit load level. 

Analyses 

The first aircraft configuration that will be examined in this 
paper is the symmetric case with the Gl.A active control 
system. Results will be shown for this case in both the open 
loop (Gl.A off) and closed loop (Gl.A on) configurations. 
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Unless otherwise noted in this paper, Ihe stochaslic simulation 
results use a maximum simulation lime of 2500 seconds, the 
TFC analysis uses simulations excited by a RMS gust intensity 
derived from a pcak-to-RMS ratio of 3.0, and the ARC method 
uses the average of the highest 5 peaks. The peak gust velocity 
U„ for tilts analysis was 112.2 fps and this corresponds to a Vr 
flight condition within the 85,0 fps design envelope crilcriori. 
In the analysis method used here, Ihe wing bending load sign 
convention is positive wing tip down. 

The results from linear stochaslic simulations and linear MFT 
analyses using the mathematical model described above are 
shown in tabular form for several loads and responses in Table 
1. This table lists the differences in the unit gust RMS responses 
between the traditional linear frequency domain PSD analyses 
and the linear time domain analyses. As shown in the table, 
all of the responses for the stochastic simulation methods arc 
within 4.5% and the results from the linear MFT analysis are 
within 2%. 

An illustration of the accuracy of Ihe state-space mathematical 
model is also shown in the comparison of frequency responses 
due to a harmonically oscillating one dimensional stationary 
gust field for two integrated loads and two physical Dcgrcc-Of- 
Freedom (DOF) accelerations. Figure 4 illustrates the fre¬ 
quency response of the wing root bending moment. Figure 5 is 
the horizontal tail root shear, Figure 6 is the C.G. vertical 
acceleration, and Figure 7 is the C.G. pitching acceleration. 
The frequency response data generated from the state-space 
model is shown with the solid line, and the dala from the tra¬ 
ditional frequency domain method is shown by the solid circles 
in these figures. As shown, the mathematical model is 
extremely accurate throughout the entire frequency range of 
interest. The higher frequency clastic modes in the analysis arc 
captured with this accurate mathematical model as well as the 
lower frequency, rigid body and clastic modes. 

A couple of representative probability functions that were used 
in the PEG analysis are shown in the next two figures. Figure 
8 illustrates the probability functions for the wing root bending 
moment in the up-bending direction and Figure 9 shows the 
down-bending direction. Included in Ihcse figures arc proba¬ 
bility functions for the linear open loop analysis (shown with the 
X's), the linear closed loop analysis (shown with the triangles), 
and the results from three nonlinear simulations with different 
input RMS gust intensities (a w = .32.0, 37.4, and 44.9 fps). These 
figures also include markings that show where the probability 
functions were evaluated for use in the PEC analysis method. 
Within the range that these probability functions were used, the 
functions are smooth. This indicates that 25D0 seconds of 
response time should be an adequate amount to be used in this 
method. 

The next two figures show the corresponding numerically 
determined values of the integrand of Equation 8 in the PEC 
analysis method. Figure 10 is the integrand for the wing root 
up-bending load and Figure 11 is the integrand for the wing 
root down-bending load. Included in these figures arc numer¬ 
ically evaluated values of the integrand for seven input RMS 
gust intensities as well as a curve representing the approxi¬ 
mation of the integrand. The numerically evaluated values of 
the integrand are shown in these figures by the solid black cir¬ 
cles and the approximation of the integrand is shown by the 
solid continuous line. The formula for the approximation of the 
integrand is shown by Equation 9. T his approximation of the 
integrand was determined using the numerically evaluated 
integrand at three of the input RMS gust intensities 
(<r w ~ .32.0, 37.4, and 44.9 Tps). The value of *2 for this symmetric 
flight, condition is 11.25 so this corresponds to <values of 


2.84, 3.32, and 4.0. The remaining four numerically evaluated 
points in this figure show that Ibis integrand is properly repres¬ 
ented with the use of this approximation. 

Intermediate results from the TEC analysis are shown in the 
next four figures. These figures illustrate the threshold 
exceedence curves for Ihe same wing root bending load. Figure 
12 is for wing root up-bending, and Figure 13 is for wing root 
down-bending. Once again included in these figures are the 
linear open loop rcsulls (shown by the X's), the linear closed 
loop rcsulls (shown by (lie triangles) and the results from the 
same three nonlinear simulations that were described above. 
Markings arc also shown in these figures to illustrate the regions 
within these curves which correspond to the limit load level. 
These exceedence curves arc not as well behaved as the prob¬ 
ability functions that were used earlier. As mentioned earlier, 
a linear regression fit was used in conjunction with these curves 
to smooth the exceedence data. 

The next two figures illustrate the accuracy of the linear 
regression fit in smoothing the data from these threshold 
cxcccdcncc curves. Figure 14 illustrates a cross-plot of Ihe wing 
root up-bending loads at. common normalized exceedence levels 
and Figure 15 illustrates Ihe same data for the wing root 
down-bending load. Ihcse curves were generated using the 
linear closed loop rcsulls and the results from the nonlinear 
simulation with a peak-lo-RMS ratio of 3.0. In these figures, 
the solid circles represent Ihe numerical rcsulls interpolated 
directly from the threshold exceedence curves. The solid line 
here reprcscnls the linear regression fit of these numerical 
rcsulls. Also shown in these figures arc markings to show the 
area of the cross-plot that is used for the limit load determi¬ 
nation. An equation is also included in these figures to repre¬ 
sent relationship between the linear and nonlinear load results. 

Intermediate rcsulls from the MFT analysis are shown next. 
Figure 16 illustrates the effect of the maximum wing up-bending 
load with the change in the impulse strength. Figure 17 shows 
the same information for the wing down-bending moment. In 
the rcsulls shown here for (he MFT analysis, 21 impulse 
strengths were used Tor the wing up-bending analysis and 9 
impulse strengths for the wing down-bending analysis. These 
impulse strengths were logarithmically spaced as shown by the 
solid circles in these two figures. Additional impulse strengths 
were used in Ihe up-bending analysis to be assured that a peak 
was not overlooked in the one dimensional search process. 
T his concern for accuracy will be discussed further later within 
this paper. 

T he next two figures illustrate a comparison of these four non¬ 
linear OF methods and include two linear methods for contrast 
as well. Figure 18 illustrates the reduction of wing bending 
moment across Ihe span of the wing for the up-bending loads. 
Figure 19 includes similar dala for the wing down-bending 
loads. As shown by Ihcse figures, the linear closed loop analysis 
(solid line, no symbols) provides the lowest (unconscrvativc) 
load levels across the span of the wing for both up and down- 
bending loads. T he linear reduced gain method shown here is 
the currently used linear approximation method. Simply put, 
the control law gain is reduced in this method until both the 
control surface rates and deflections arc within the prescribed 
limits. For the rale limitations at. the limit load level, the no- 
load control .surface rates have been used. As we can sec in 
these figures, for this symmetric case, this method results in Ihe 
highest (most conservative) loads. 

flic PIT’ and T EC results are very similar for both the wing 
up and down-bending loads and are within 3.5% across the 
span of the wing. T he A PC method for the wing up-bending 
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loads is within a couple percent of the PF.C method, but for the 
down-bending loads, the A PC method produced loads that arc 
up to 12% lower (less conservative). 

The MFT results arc up to 5.5% lower than the PFC results for 
the wing up-bending analysis and up to 3.5% lower for the 
down-bending analysis. Many points were added to the ori¬ 
ginal nine impulse strengths in the wing up-bending MFT 
analysis to provide assurance that the critical matched 
excitation waveform was found. The final number of impulse 
strengths used was 21. The MFT analysis will always approach 
the nonlinear limit load from the unconscrvativc side, due to the 
constraint of the energy of the matched excitation waveform. 
In the nonlinear MFT one-dimensional search method, there 
is no guarantee that the critical matched excitation waveform 
will ever be found. It is believed in the MFT analysis shown 
here, that the critical matched excitation has not yet been 
obtained. Further use of another nonlinear MFT procedure 
such as the multi dimensional optimization should show that 
both the wing up and down-bending loads will increase. 

The A PC method used here resulted in unconscrvativc loads for 
tiic wing down-bending moments and further investigation into 
this behavior was needed. In the APC method, the 
pcak-to-RMS ratio is recalculated for each input gust intensity. 
For the analysis shown here, the pcak-to-RMS ratio for the 
wing up-bending loads remained fairly constant with changes 
in the input gust intensity. For the wing down bending loads this 
was not the case. In the wing down-bending analysis, the 
pcak-to-RMS ratios increased as the RMS gust intensity 
increased. 

A modification to the APC criteria was included that generated 
the peak gust velocity U„ from the peak to RMS ratio calcu¬ 
lated in the linear closed loop analysis. This modification to the 
APC method can, in certain special circumstances, be identical 
to the TFC method. The results from this modification are 
shown for the wing up-bending loads in Figure 2D and for the 
wing down-bending loads in Figure 21. As shown in these fig¬ 
ures, the results from the APC method wing up-bending loads 
arc relatively unchanged. The results from the APC method 
wing down-bending loads have increased and arc now within 
1.1% of the loads from the PFC method. 

An investigation into the effect of the simulation gust intensity 
used in the TFC method is illustrated in the following two fig¬ 
ures. In these two figures, the TFC method is analyzed using 
three input gust intensities and is compared to the PFC method. 
For the up-bending loads of Figure 22, the results from the 
TFC method with a pcak-to-RMS ratio of 3.5 almost exactly 
matched the results from the PFC method. For the TFC results 
from the other two pcak-to-RMS ratios, the results were about 
3-5% higher (more conservative) across the span of the wing 
For the down-bending loads shown in Figure 23, the results 
from the TFC method using a pcak-to-RMS ratio of 3.0 were 
almost identical to those from the PFC method. The results 
from the analyses using the remaining two pcak-to-RMS ratios 
were with ± 1-2% of the PFC results. 

An investigation into the amount of stochastic simulation time 
needed to obtain reasonable results from the TFC method was 
also performed. The results of this investigation arc shown in 
the following two figures. In Figure 24, the efFcct of (he max¬ 
imum simulation Lime on the wing up-hending loads is illus¬ 
trated. As shown in this figure, results from the 2500, 2000, and 
1500 second simulations provide extremely similar results. 
When the maximum simulation time is reduced to 100 or 500 


seconds, the results diverge and become about 5-6% less con¬ 
servative. 

A similar investigation is shown in Figure 25 for the wing 
down-bending loads. Once again, for the maximum simulation 
time of 2500, 2000, and I 50(5 seconds, the results are indistin¬ 
guishable. For simulation limes of 1000 and 500 seconds, the 
results also begin to diverge and become slightly less conserva¬ 
tive. 

Tile second aircraft configuration that is examined in this paper 
is an anti symmetric case with a yaw damper (YD) active con¬ 
trol system. Results arc shown Tor this case in the closed loop 
(YD off) configuration only. Once again, unless otherwise 
noted in this paper, the stochastic simulation results use a 
maximum simulation time of 2500 seconds, the TEC analysis 
uses simulations excited by a RMS gust intensity derived from 
a pcak-to-RMS ratio of 3.0, and the APC method uses the 
average of the highest 5 peaks. The peak gust velocity U„ for 
this analysis was also 112.2 fps corresponding to a Vs flight 
condition within 85.0 fps design envelope criterion. 

Figure 26 illustrates a comparison of the four nonlinear CT 
analyses along with two linear analyses for the vertical tail 
integrated normal shear forces. Once again, the linear closed 
loop results arc the lowest and most unconscrvativc over most 
of the vertical tail span. Flic results from the linear reduced 
gain method result in integrated loads that arc deemed to be the 
most conservative of the methods shown here. All of the non¬ 
linear methods remain within about 3% across the span of the 
vertical tail except for the results from the APC method which 
arc 8-1 2% lower than the other nonlinear methods in the region 
of the lower vertical tail. 

Another investigation into the discrepancy between the APC 
and the other (I'FC.PFC, and MFT) nonlinear CT analyses 
was performed. It was found, once again, that the 
pcak-to-RMS ratio of the nonlinear simulations increased at the 
higher gust intensities. The APC method was modified again 
to use the peak lo-RMS ratio from the linear closed loop anal¬ 
ysis and the results arc shown in Figure 27. As shown in this 
figure, this modification to the APC method brought the results 
in line with the other three nonlinear CF methods. 

Figure 28 illustrates the effect of different input RMS gust 
intensities used in the T FC method on the vertical tail normal 
shear forces. Figure 29 shows tfic effect of reducing the amount 
of simulation lime used in the TFC method on the vertical tail 
normal shear forces. As shown in these two figures, the effect 
of input RMS gust intensity and simulation time for the anti¬ 
symmetric case are minimal within tficsc reasonable bounds. 

Computing Performance Issues 

One of the major issues when discussing the use of these non¬ 
linear gust analysis methods for use with 'real world' problems 
is that of computing time (and effort) needed to analyze these 
problems. It is important to point out that for the APC, PF.C 
and TFC' methods that the computing lime is relatively constant 
for any number of loads or responses that arc evaluated, 
because the loads and responses can typically be determined 
at a much coarser time interval than the integration of the 
equations of motion. For the MFF analysis, the opposite is 
true. The computing time needed for the MFT analysis is 
almost linear to the number of output loads or responses that 
arc requested in the analysis. This is because in the MFT 
analysis, the aircraft model is excited by a series of matched 
excitation waveforms that are derived from the impulse 
response functions of each load or response. 
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Included below is a brief description of (he computing lime 
needed to produce the results shown here. This description is 
not meant to necessarily be conclusive but instead should be 
used as a tool to show the relative amount of computing time 
that could be needed for real problems. The amount of simu¬ 
lation time required for the A PC method is highly subjective 
and will not be discussed here. 

For the symmetric gust analysis, Ihe analysis included 15 wing 
load stations. 3 load components were determined at each of 
these load stations. 7 physical DOF accelerations were also 
included within this analysis. In the MFT analysis, Ihe wing 
up-bending load analysis was performed using 21 input impulse 
strengths and the wing down bending analysis used 9 input 
impulse strengths. The simulation of Ihe first system in the 
MFT analysis was performed for 10 seconds using the impulse 
function and the simulation of the second system was per¬ 
formed for 15 seconds using each of the matched excitation 
waveforms. When we add up all of the simulation lime needed 
for this analysis we can see that this equals (21 + 9)*1*I0 + 
(21+9)*(45 + 7)*15 = 23,700 simulation seconds. The first 
term in the equation is the amount of simulation time needed 
in the first system and the second term is the amount of simu¬ 
lation time needed in Ihe second system. 

For the PFC method, four simulations were needed. Three of 
the simulations were needed for the nonlinear analysis of the 
PFC method and one linear baseline simulation was required 
to correlate the time domain method to Ihe results from the 
linear frequency domain methods. The amount of lime for the 
PFC method was 2500 seconds per simulation. Thus the total 
amount of simulation time required for the PFC method shown 
here was 10,000 seconds. 

For die TFC method, two simulations were needed. One linear 
baseline and one nonlinear simulation. The results from the 
TFC method showed that the reasonable results were obtained 
when simulation times of 1500 to 2500 seconds were used. 
Thus for the TEC method, the amount of simulation time 
ranges from 3000 to 5000 seconds. 

For the anti-symmetric case shown in this paper, the computing 
time used for the PEC and TEC analyses was the same as that 
used in the symmetric analysis. For the MFT analysis, the 
amount of simulation time was much lower. This was due to 
the fact that the number of loads and responses under analysis 
was reduced (7 vertical tail load stations, 3 components of load 
per station, and 8 physical DOF accelerations), only I I impulse 
strengths were used in the one-dimensional search process 
(instead of 21 and 9 for the symmetric analysis), and also due 
to the symmetry of this nonlinearity. When we add up Ihe 
simulation time needed for this analysis we can see that this 
equals 11*1*10 + 11*(21+8)*I5 = 4,895 simulation seconds. 

A tabular comparison of (he simulation time required for these 
analyses is shown in Table 2. It is important to note here that 
the analyses performed for this paper included only a small 
subset of the loads and responses on the aircraft. Typical pro 
duclion gust analyses will require the analysis of many more 
responses and the computing time for the MFT analysis will 
increase accordingly. 

Conclusion 

The nonlinear continuous turbulence analyses arc currently 
extremely expensive to use. More research is needed into 
methods that reduce the cost of these analyses. The use of the 
nonlinear gust analysis methods should only be after per¬ 


forming an exhaustive filtering of Ihe flight conditions using 
linear methods and approximation techniques. 

For the analyses shown here in this paper, the TEC method 
illustrated a slightly conservative behavior while requiring the 
least amount of simulation lime. It is believed that the PEC 
method is the most accurate of those shown in this paper but 
ihe computing cost of Ibis method is at least double that of the 
TFC method. Comparisons of the two methods should be 
performed for a representative condition to be assured of the 
conservative nature of the T FC method. A sensitivity analysis 
should also be performed both with gust intensity and amount 
of simulation time that is required in the T FC method. 

Correlated loads arc difficult if not impossible to determine 
using methods based on the T FC and PEC methods. Corre¬ 
lated load conditions may be determined using die MFT anal¬ 
ysis methods for a small group of loads or responses. Some 
scaling of the results from Ihe MFT analyses will need to be 
performed to eliminate Ihe unconscrvativc nature of the anal¬ 
ysis. 
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Figure 1. Flowchart of the Linear MFT Process Figure 2. Flowchart of the Non-Linear MFT Process 
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Table 1. Comparison of Time Domain and Frequency 
Domain Results 
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Table 2. Comparison of Simulation Time Needed for 
Nonlinear Gust Analyses 



Figure 3. Doublet Lattice Model of a Large Transport 
Aircraft 
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Figure 5. Frequency Response of Horizontal Tail Root 
Shear from a Vertical Gust 
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Figure 4. Frequency Response of Wing Root Bending Figure 6. Frequency Response of C.G, Vertical 
Moment from a Vertical Gust Acceleration from a Vertical Gust 
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Figure 10. Approximation of the PEC Probability 
Integral for Wing Root Up-bending 
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Figure 7. Frequency Response of C.G. Pitching 
Acceleration from a Vertical Gust 
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Figure 11. Approximation of the PEC Probability 
Integral for Wing Root Down-bending 
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Figure 8. Cumulative Probability Functions for Wing Exceedence Curves for Wing Root Up- 
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Figure 9. Cumulative Probability Functions for Wing Figure 13. Exceedence Curves for Wing Root Down- 
Root Down-bending Bending 
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Figure 14. Linear Regression Fit of the Exceedence 
Data for Wing Root Up-bending 
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Figure 18. Comparison of Four Nonlinear Gust Ana¬ 
lyses for Wing Up-Bending Loads 



Figure 15. Linear Regression Fit of the Exceedence 
Data for Wing Root Down-bending 



Figure 16. Effect of Impulse Strength on Nonlinear MFT 
Wing Root Up-Bending Load 



Figure 17. Effect of Impulse Strength on Nonlinear MFT 
Wing Root Down-Bending Load 
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Figure 19. Comparison of Four Nonlinear Gust Ana¬ 
lyses for Wing Down-Bending Loads 
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Figure 20. Comparison of Four Nonlinear Gust Ana¬ 
lyses for Wing Up-Bending Loads (with 
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Figure 21. Comparison of Four Nonlinear Gust Ana¬ 
lyses for Wing Down-Bending Loads (with 
Modified APC Method) 
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Figure 22. Effect of Input RMS Gust Intensity on TEC 
Results for Wing Up-Bending 
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Figure 23. Effect of Input RMS Gust Intensity on TEC 
Results for Wing Down-Bending 


Figure 


S — Unew O Nonlinear (TEC) £ 

g Q Nonlinear (PEC) A Nonlinear (MET) g 

ja O Nonlinear (Modified APC) 96 Linear Reduced Cain 5- 
1 r i-1-i-i-i-i i i i 

50 100 150 200 250 300 350 400 450 500 

Vertical Load Station Coordinate 
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lyses for Vertical Tail Shear Forces (with 
Modified APC Method) 
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Figure 24. Effect of Simulation Time on TEC Results 
for Wing Up-Bending 
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Figure 25. Effect of Simulation Time on TEC Results 
for Wing Down-Bending 
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Figure 28. Effect of Input RMS Gust Intensify on 
Results from TEC Method for Vertical 
Tail Shear Forces 
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1. SUMMARY 

Statistical Discrete Gust (SDG) methodology for extreme tur¬ 
bulence provides the basis of a design criterion that accounts 
for: (1) the statistical nature of design level turbulence; (2) the 
effects of lightly damped aircraft responses; (3) critical re¬ 
sponses that may have significant frequency separation; and 
(4) the effects of nonlinear control systems. The criterion 
would eliminate the need for multiple design criteria and 
could increase the structural efficiency of new aircraft. The 
increased efficiency would be the result of having equally 
probable design limit loads for all gust critical components of 
the aircraft. 

2. LIST OF SYMBOLS 

A Standard deviation 

E- i th Extreme value of a response 
EAS Equivalent airspeed 

F Flight profile alleviation factor 
if Gust development length, ft 

P k SDG probability factor 

RB Rigid Body 

U Maximum velocity for SDG ramp gusts, ft/sec (TAS) 
Uj Intensity factor for JAR gusts , ft/sec(TAS)/ft 1/6 
U ref Reference gust velocity, ft/sec (EAS) 

U 0 SDG gust intensity factor, 

ft/sec(TAS)/ft 1/6 for extreme turbulence 
ft/sec(TAS)/ft 1/3 for continuous turbulence 
U c Continuous turbulence design intensity, ft/sec 

TAS True airspeed 

y SDG design amplitude 

O air density ratio 

3. INTRODUCTION 

Analyses of a swept wing business jet with fuselage mounted 
engines has raised questions about existing design gust crite¬ 
ria. Areas have been identified where existing criteria may 
not provide levels of safety comparable to those for aircraft 
having frequency response characteristics that differ from 
those of the business jet. To alleviate this concern a design 
criterion has been developed based upon the SDG model of 
extreme turbulence. This development has followed verifica¬ 
tion of the equivalency of power spectral density (PSD) and 
SDG results for linear systems in continuous turbulence 1 and 
the extension of the SDG model to include extreme turbu¬ 
lence 2,3 . This paper concentrates on the results obtained from 
SDG methodology for a particular aircraft and comparison of 
the results with those from present criteria. However, a brief 
discussion of SDG-PSD equivalency is included since the 
continuous turbulence methodology is directly applicable to 
extreme turbulence. 


3.1 Continuous Turbulence 

Present criteria specify use of the von Karman turbulence 
spectrum with 2,500 feet as the characteristic scale of turbu¬ 
lence. With this model for PSD analyses and unity for the 
turbulence intensity in each analysis, the relationship between 
PSD and SDG results is given by: 

y = 10.4 A (1) 

Equation 1 is consistent with units of feet and seconds. Other 
units of distance will require adjustment of the constant, 10.4. 
The constant must be adjusted by the cube root of the conver¬ 
sion factor between feet and the desired unit of measure. This 
relationship has been satisfied to a reasonable degree of 
accuracy by deleting one assumption included in references 1 
and 4. The assumption that a gust profile consisting of either 
one, two, four, or eight (1-cos) ramp gusts would adequately 
define the critical profile for a given load or response has been 
found to be inappropriate. The number of ramps in a critical 
profile may take on any value, and may even be greater than 
eight for lightly damped modes. 

The basic SDG model represents continuous turbulence as a 
series of non-overlapping (1-cos) ramp gusts of alternating 
sign. The maximum velocity for the j th ramp gust is a func¬ 
tion of gust intensity, U 0 , and the gust development length, Hj. 

Uj = U 0 Hf (2) 

Using this relationship and the principle of superposition for 
linear systems, aircraft responses are calculated for a family of 
ramp gusts of varying development length. For each desired 
load and response the extreme values of alternating sign, Ej, 
resulting from the family of ramp gusts are identified. The 
identification of these extreme values requires a peak tracking 
routine. However, that routine will not be described here 
since peak tracking can be accomplished by several methods 
and it is not necessary when using the SDG model of extreme 
turbulence. 

Once the extreme values have been determined, they are se¬ 
lected in descending order of magnitude with alternating sign. 
The times of occurrence and the gust development lengths 
associated with each extreme value are retained to form the 
critical profile and to assure the ramps do not overlap. The 
design amplitude, y , for a given load or response is the 
maximum of the peak responses, y k , resulting from profiles 
consisting of k ramps. 

Yk= P klM for k = 1 t0 n ( 3 ) 

i=l 


Paper presented at a Workshop on “Aircraft Loads due to Turbulence and their Impact on 
Design and Certification” held in Lillehammer, Norway on 5 May 1994. 
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The probability factors, P k , take on values of 1.0 for k=l and 
l/(0.88Vk) for k>l. The constant, n, is the number of ex¬ 
treme values retained from the response analyses. The prob¬ 
ability factors maintain a constant probability of exceedance 
for gust patterns consisting of multiple ramps. The design 
amplitude is then: 

y = max(Y 1 ,y 2 , y 3 „..,Y n ) (4) 

Development lengths for the family of gusts must cover a suf¬ 
ficient range to assure finding maximum values of y for each 
of the loads and responses. For vertical gusts a development 
length range of 30 to 500 feet is usually adequate. Lateral 
gust analyses may require development lengths greater than 
1000 feet to accommodate lightly damped dutch roll modes. 
Each gust response analysis must extend over a sufficient pe¬ 
riod of time to provide the necessary number of extreme val¬ 
ues for calculation of y . The duration of each analysis may 
vary from two to three seconds for vertical gusts; while it may 
be 15 seconds or more for lateral gust analyses. 

Equation 1 is satisfied within two percent when unit values of 
U 0 and U a arc used and gust development lengths arc chosen 
such that the gust velocity increases approximately two per¬ 
cent for each succeeding ramp gust in the SDG analysis. Thus 
the validity of SDG methodology for linear systems in con¬ 
tinuous turbulence has been verified. 

Verification of the equivalency between SDG and PSD analy¬ 
ses of linear systems in continuous turbulence and the exten¬ 
sion of the SDG model to extreme turbulence has led to a 
practical means of predicting aircraft limit loads reflecting the 
statistical nature of design level atmospheric turbulence. 

3.2 Extreme Turbulence 

Extension of the SDG model to extreme turbulence requires an 
alternate power law for the maximum gust velocity and a set 
of probability factors that reflect the intermittent nature of the 
turbulence. For extreme turbulence the relationship between 
the maximum gust velocity and development length for the j th 
ramp gust is given by: 

Uj = Uo-Hj /6 (5) 

At present, provisional values (obtained from J .G. Jones, 
DRA Farnborough) arc being used for the factors that 
maintain a constant probability of exceedance for multiple 
ramp profiles. The values for P,, P 2 and P 3 are 1.0, 0.75 and 
0.48, respectively. The number of ramp gusts forming a criti¬ 
cal profile is restricted to cither one or two; since 

0-75 i|E;| >0.48 X|Ej| (6) 

i=l i=l 

for the worst case, Ej=E 2 =E v This restriction to the gust 
shape leads to simplification of the SDG methodology de¬ 
scribed above for continuous turbulence. A maximum of two 
extreme values is required to define critical gust profiles and 
these values must be of opposite sign to satisfy the require¬ 
ment for ramp gusts of alternating sign. Thus, overall maxi¬ 
mum and minimum values from the complete family of gusts 
analyzed arc all that arc required. This eliminates the need for 


peak tracking routines and reduces the magnitude of the com¬ 
putations for SDG gust analyses to a level similar to that re¬ 
quired for the JAR 25.341, Amendment 91/1 discrete gust. 
Computation time is increased slightly for lateral gust analy¬ 
ses to accommodate the dutch roll response of the aircraft. 
Details of SDG methodology for extreme turbulence are out¬ 
lined in the Appendix. The following discussion of gust load 
analyses based upon present discrete gust and continuous tur¬ 
bulence criteria and a SDG criterion illustrates the benefits 
provided by a statistical model of extreme turbulence. 

4. DISCUSSION 

Gust load analyses using continuous turbulence methods 
(either SDG or PSD) maximize contributions to structural 
loads from rigid body and elastic responses of the aircraft. 
However, there are two undesirable characteristics resulting 
from using a model of continuous turbulence for design limit 
loads. First, for relatively high frequency responses that are 
critical for short gust development lengths, the design values 
for U c yield gust velocities that are much lower than what is 
required to represent design levels of extreme turbulence. 
Again, it is the extreme turbulence that has the potential to 
cause structural overload and failure. Secondly, for lightly 
damped low frequency responses that arc critical for long gust 
development lengths the continuous turbulence model yields 
extremely high gust velocities for design values of U 0 . The 
combination of light damping and high gust velocities yield 
exaggerated values for limit loads. 

Scaling the continuous turbulence model to design levels does 
not account for changes in the nature of turbulence with in¬ 
creasing intensity. Incremental gust velocities associated with 
the continuous turbulence model increase rather slowly with 
increasing development length (H 1A1 power law). On the other 
hand, velocities associated with extreme turbulence increase 
rapidly at short development lengths but at a much lower rate 
for longer development lengths (H l/6 power law). Unless the 
character of extreme turbulence is accounted for, gust induced 
limit loads will not have equal probability of exceedance. 

Results presented in this section are based upon gust loads 
analyses of a swept wing business jet with fuselage mounted 
engines. Two other aircraft (a straight wing business jet and a 
transport with four wing-mounted engines) have been 
analyzed. When the differing response characteristics arc 
taken into consideration, SDG analysis results for those 
aircraft are consistent with the results presented below. 
Analyses are based upon JAR 25.341 Amendment 91/1 
discrete gust and FAR Part 25 Appendix G continuous 
turbulence requirements, and a SDG model with the gust 
intensity factor derived from U ref , the JAR 25.341 reference 
gust velocity. SDG modeling provides a statistical 
representation of the turbulence that caused the extreme events 
from which values of U ref were derived. The SDG extreme 
turbulence model is applied for design speeds V 0 , V c and V,,, 
with the gust intensity reduced 50 percent for speed V,,. 
Pending further study, design speed V R is retained to assure 
that the effects of very lightly damped rigid body modes are 
taken into account. However, the design gust intensity is 
unchanged from that for speed V c . This will be discussed fur¬ 
ther in the paragraphs on loads due to lateral turbulence. 

Design gust intensities for discrete gust loads analyses arc 
compared in terms of gust intensity factors U Q and U, for the 
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SDG and JAR discrete gusts, respectively. The factors are 


defined in terms of true airspeed by: 


U 0 = Uref / 350 1/6 /Vo 

(7) 

and 


Uj = F g • U ref / 350 1/6 ! 

(8) 


The design gust intensity factors are illustrated in Figure 1 
with the probability factors applied to the SDG intensities. 
Even though intensities for JAR gusts are greater than those 
for the two-ramp SDG profile, the resulting loads are consis¬ 
tently greater for the SDG method. The shape of the two- 
ramp SDG profile has a significant influence on the magnitude 
of the gust induced structural loads. The altitude variation of 
design gust intensities for the two discrete gust models also 
contributes to the differences in the design limit loads pre¬ 
sented in the following paragraphs. 

Comparison of design level extreme and continuous turbu¬ 
lence is facilitated by using SDG models for both types of tur¬ 
bulence. Gust intensity factors equivalent to the design values 
of U a for continuous turbulence PSD analyses are given by: 

U 0 =iyi0.4 (9) 

The gust intensity factors for extreme and continuous turbu¬ 
lence have different units. Therefore, the following compari¬ 
sons will be in terms of true gust velocity. Present design 
level velocities for continuous turbulence and those selected 
for use with the SDG extreme turbulence model are compared 
in Figure 2. The design gust velocities are for single ramp 
gusts having development lengths of 350 feet. The variation 
of design velocities with altitude illustrated in Figure 2 is a 
major contributor to differences in limit loads resulting from 
SDG extreme and continuous turbulence models. 



Figure 1 Design Gust Intensity Factors for Discrete Gust 
Analyses 

For vertical turbulence, an upper limit for critical ramp 
lengths is approximately 350 feet. For lateral turbulence 
critical ramp lengths tend to be somewhat greater and may ex¬ 
ceed 1000 feet for both extreme turbulence and continuous 
turbulence. The actual critical development lengths will de¬ 
pend upon the amount of damping present in the dutch roll 
mode. For a development length of 1000 feet, extreme turbu¬ 
lence maximum gust velocities are 19 percent greater than 
those shown in Figure 2; while for continuous turbulence 


maximum velocities are 42 percent greater than those shown. 
The long gust development lengths for continuous turbulence, 
the associated high gust velocities and light damping of the 
dutch roll mode yield very conservative limit loads for com¬ 
ponents such as the aft fuselage and vertical fin. 
Representative loads for the wing and horizontal stabilizer due 
to vertical turbulence and vertical fin loads due to lateral 
turbulence are compared in the following paragraphs. 



Figure 2 Example Design Gust Velocities for Extreme 
and Continuous Turbulence 

4.1 Vertical Turbulence 

Discrete gust analyses yield maximum positive wing bending 
(wing tips up) for flight at maximum dynamic pressure and 
maximum mach for speed V c . For this flight condition, SDG 
methodology was used for continuous turbulence analyses. 
This provides results in a form comparable to those from dis¬ 
crete gust analyses. Critical gust profiles for this flight condi¬ 
tion are shown in Figure 3 for the two discrete gust models 
and for continuous turbulence. The same maximum gust 
velocity was chosen for the continuous turbulence analysis as 
that for the SDG extreme turbulence analysis. This provides a 
direct comparison of the two analyses without introducing 
variations due to criteria. Using the same maximum gust 
velocity for both analyses serves to illustrate the similarities 
and difference between the two models. 



Figure 3 Critical Gust Profiles for Wing Root Vertical 
Bending Moment, Positive Tip Up 
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It can be seen from Figure 3 that critical SDG profiles for 
continuous turbulence and extreme turbulence have similar 
shapes. Each profile begins with a relatively long ramp and 
ends with a shorter ramp. The long initial ramp and the 
overall duration of the gust provide maximum excitation of 
the short period mode; while the shorter ramp provides maxi¬ 
mum excitation of the elastic modes. Gust velocity envelopes 
corresponding to these profiles are illustrated in Figure 4. 



Gust Development Length - Feet 

Figure 4 Gust Velocity Envelopes at the Altitude for 
Critical Discrete Gust Loads 

Wing root bending moments due to the critical SDG profiles 
arc illustrated in Figure 5. The variation of the bending mo¬ 
ment with gust penetration distance illustrates the differences 
between the extreme and continuous turbulence models. 
Compared to the extreme turbulence results the continuous 
turbulence model produces lower bending moment contribu¬ 
tions from both the short period mode and the elastic modes. 
The loads from the extreme turbulence analysis are 14 percent 
greater than those from the continuous turbulence analysis. 
This, even though the maximum gust velocities are the same. 
The continuous turbulence intensity in this example corre¬ 
sponds to a U a value of 86.2 ft/scc (which is 27 percent 
greater than that required for design speed V c by FAR Part 25 
Appendix G). 



Figure 5 Normalized Incremental Wing Bending Moment 
Due to Critical SDG Profiles 


The bending moments shown in Figure 6 for the JAR analysis 
show no elastic mode contribution resulting from the profile 
yielding the maximum bending moment. This is due to the 
significant separation between the frequencies of the elastic 
modes and the rigid body short period mode. The bending 
moment response for the maximum contribution from the 
elastic modes yields only 91 percent of the moment resulting 
from the critical JAR profile. SDG methodology, on the other 
hand, maximizes the contributions of both rigid body and 
elastic modes of the aircraft. As a result the maximum 
bending moment from the SDG analysis is 16 percent greater 
than that from the JAR analysis. The SDG loads are greater 
even though the JAR gust velocity is 19 percent greater than 
that for the SDG analysis. These results illustrate the 
importance of gust shape in the determination of design limit 
loads. 



Figure 6 Normalized Incremental Wing Bending Moment 
Due to SDG and JAR Discrete Gust Profiles 

For aircraft with lower frequency wing bending modes the 
JAR profile may yield significant elastic contribution to criti¬ 
cal incremental loads. For example, JAR analyses yielded in¬ 
cremental wing bending moments for the four-engine transport 
that were slightly greater than those from the SDG analysis. 
However, the maximum velocity for the JAR gust was ap¬ 
proximately 20 percent greater than that for the SDG critical 
profile. The JAR profile will usually yield wing bending mo¬ 
ments greater than those from the SDG extreme turbulence 
model only if the maximum gust velocity is significantly 
greater than the velocity for the two-ramp SDG profile. 
However, depending upon the response characteristics of 
particular aircraft, some wing bending moments arc critical 
for a one-ramp SDG profile. In that case the velocity for the 
SDG profile will generally be greater than that for the JAR 
profile. 

4.1.1 Design Limit Loads 

Normalized design limit loads at the root of the wing and hori¬ 
zontal stabilizer arc shown in Figures 7 through 10 for the 
three analysis methods. All loads have been normalized by 
the results from SDG analyses. The majority of the SDG loads 
due to vertical gusts arc greater than those from JAR and PSD 
analyses. Exceptions to this occur because the design levels 
of turbulence for the three turbulence models vary sig¬ 
nificantly with altitude. Because of these altitude variations, 
the JAR and PSD results for torsion, shown in Figure 7, ex¬ 
ceed those from the SDG analysis. The critical value of tor- 
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sion from the JAR analyses occurs at high altitude where the 
maximum gust velocity is much greater than that for the criti¬ 
cal two-component SDG profile. The critical positive torsion 
from PSD analyses occurs at low altitude and speed V B where, 
according to FAR Part 25 Appendix G, U a has a value of 99 
ft/sec. 
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Figure 7 Normalized Positive Wing Root Limit Loads 
Due to Vertical Turbulence 



The negative wing loads shown in Figure 8 are also influenced 
by the distribution of design gust intensities versus altitude. 
Down loading on the wing tends to be critical for light weight 
configurations flying at maximum dynamic pressure (sea level 
and design speed V c ). The SDG gust intensity decreases very 
little at low altitude as can be seen in Figure 1. However, at 
low altitude, the JAR gust velocities are reduced significantly 
by the flight profile alleviation factor. For this particular 
aircraft, down loading of the wing due to the SDG extreme 
turbulence model tends to be critical for single-ramp gusts. 
As can be seen from the data in Figure 1, SDG intensities for 
single-ramp profiles are greater than the JAR intensities at all 
altitudes. At low altitude the JAR design gust intensities are 
relatively low which causes down loading of the wing to 
become critical at a higher altitude. Because of the critical 
single-ramp SDG profile and the low JAR gust intensities at 
the lower altitudes, JAR loads are somewhat less than those 
from SDG analyses. The light weight configurations have 
minimum fuel in the wing. The reduced mass yields increased 
frequencies for the elastic modes. PSD loads are reduced 
because of these increased frequencies. The reduced critical 
gust development lengths associated with the increased 
frequencies yield relatively low design gust velocities . 

The negative horizontal stabilizer loads (downward loading) 
illustrated in Figure 9 follow the same trend as the positive 
wing loads. Here again the torsional moments from the JAR 
analysis exceed the SDG results. The critical condition for 
JAR torsion is at high altitude, where the JAR gust velocity is 
as much as 33 percent greater than that for the critical two- 
ramp SDG profile. PSD loads are relatively low because of 
the relatively high frequencies of the elastic modes. 

The variation of steady-state balancing tail loads with flight 
condition has a significant influence on the positive (upward 
loading) horizontal stabilizer loads illustrated in Figure 10. In 
this case the net load represents the difference between steady- 


state and incremental loads. The JAR loads are critical for 
higher altitudes where the gust intensities are greater, while 
the SDG and PSD loads are critical for lower altitudes. 
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Figure 8 Normalized Negative Wing Root Limit Loads 
Due to Vertical Turbulence 
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Figure 9 Normalized Negative Horizontal Stabilizer Root 
Limit Loads Due to Vertical Turbulence 
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Figure 10 Normalized Positive Horizontal Stabilizer Root 
Limit Loads Due to Vertical Turbulence 
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4.2 Lateral Turbulence 

SDG continuous turbulence analyses were performed for ver¬ 
tical fin lateral bending moment to provide data comparable to 
extreme turbulence results. The critical condition for bending 
moment at the root of the vertical fin due to continuous turbu¬ 
lence is at low altitude and design speed V B . Typical critical 
profiles for the two types of turbulence arc illustrated in 
Figure 11. The continuous turbulence profile is made up of 
six ramps having development lengths of 1,363 feet and one 
ramp with a length of 912 feet. On the basis of a U 0 value of 
99 ft/scc the maximum gust velocity for the longer ramp is 
105 ft/scc and 92 ft/sec for the short ramp. For a gust profile 
consisting of seven ramps, the probability factor, P ? , reduces 
the maximum incremental velocities to 45 and 40 ft/sec 
respectively, for the long and short ramps. In contrast, the 
extreme turbulence critical profile is formed from one ramp 
with a development length of 997 feet, followed by a second 
ramp with a length of 89 feet. The maximum velocity asso¬ 
ciated with this profile is 53 feet per second. 



Figure 11 Typical SDG Critical Profiles for Fin Lateral 
Bending Moment 

The bending moments resulting from the critical profiles are 
illustrated in Figure 12. For continuous turbulence, maximum 
bending moment occurs at a distance of 9,528 feet into the 
gust pattern. The long development lengths in the critical 
continuous turbulence profile cause the critical bending mo¬ 
ments to be entirely due to rigid body response characteristics 
of the aircraft. In contrast, the extreme turbulence profile 
provides excitation of both the elastic modes and rigid body 
modes. The maximum bending moment for extreme turbu¬ 
lence occurs at a gust penetration distance of 2,213 feet, com¬ 
pared to the distance of 9,528 feet for the critical continuous 
turbulence profile. The JAR gust was not included in this 
comparison because the maximum development length of 350 
feet fails to provide significant excitation of the dutch roll 
mode. 

Considering the bending moment data for the wing (Figure 5) 
and the vertical fin data in Figure 12; the question must be 
raised as to whether or not gust loads analyses based on a 
model of continuous turbulence should be required for design 
limit loads. While the continuous turbulence model (at inten¬ 
sity levels well below those being used for design limit loads) 
is appropriate for defining fatigue loads, it clearly docs not 


represent the nature of the extreme turbulence that is likely to 
cause structural overload and failure. 
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Figure 12 Normalized Vertical Fin Lateral Bending 
Moment Due to Critical SDG Profiles 

4.2.1 Design Limit Loads 

The normalized design limit loads in Figure 13 for the vertical 
fin illustrate the influence of a lightly damped dutch roll mode 
on design limit loads. The SDG model of extreme turbulence 
accommodates the light damping of the dutch roll mode by 
virtue of the variable gust shape. Vertical fin design limit 
loads from the extreme turbulence analyses have the same 
probability of exceedance as those for the other aircraft 
components. The results from PSD analyses illustrate the 
effects of complex gust profiles and the large values of U n at 
low altitude (99 ft/scc for design speed V (1 ). The continuous 
turbulence intensity level for design speed V is 32 percent 
greater than that for speed V c „ If this increased intensity was 
eliminated for speed V B , the difference between the SDG and 
PSD limit loads in Figure 12 would be reduced significantly. 
Elimination of the increased intensity for speed V H does not 
imply the elimination of speed V g as a design condition. Very 
lightly damped rigid body modes may cause critical loads at 
speeds lower than V c „ This is an area that should receive ad¬ 
ditional studv before considering the elimination of sneed V, 
as a d 
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Figure 13 Normalized Fin Root Limit Loads Due to 
Lateral Turbulence 
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4.3 Design Gust Intensity 

The gust intensity factors used with the SDG model of ex¬ 
treme turbulence were derived from, U ref , the reference gust 
velocity of JAR 25.341 Amendment 91/1. The reference gust 
velocity was derived from an accumulation of extreme event 
data collected during normal commercial airline operations. 
This derivation included the design cruise speeds of aircraft 
that normally cruise in various altitude bands. The use of 
these design cruise speeds causes a deviation from the statis¬ 
tics of the measured data; but results in conservative values 
for design gust intensities. The limit loads due to SDG ex¬ 
treme turbulence presented in previous sections have been de¬ 
scribed as having equal probability of exceedance for all 
components of the aircraft. However, this is true only at any 
given altitude, since the gust intensity factors were derived 
from U ref . 

Design gust intensity factors would more accurately reflect the 
statistics of the measured data if the flight speeds and the time 
spent in a specific altitude range by the data collecting aircraft 
were used in the derivation. The design values would then re¬ 
flect the statistics of the recorded events. Thus, the altitude 
variation of design gust intensities would reflect the statistical 
nature of extreme turbulence as weighted by normal airline 
operations. 

As stated earlier, the probability factors used with the SDG 
model of extreme turbulence are provisional. These values 
should be re-evaluated and revised (if required) for use in an 
official SDG based design gust criterion. 

4.4 Nonlinear systems 

The preceding discussion was restricted to procedures and re¬ 
sults for analysis of linear systems. The process, however, is 
readily extended to nonlinear systems. SDG methodology for 
extreme turbulence is well suited for the development of gust 
load alleviation systems since it readily accommodates control 
system nonlinearities. The computational effort increases 
significantly because the search for worst case loads must 
span the applicable range of gust development lengths for each 
of the two possible ramp gusts and the distance between the 
ramps. However, computational effort can be reduced by 
combining a systematic (hill-climbing) search for the worst 
case with an element of randomization so the search will not 
stop at a local maxima, rather than finding the global max¬ 
ima 4 . 

Nonlinear automatic control systems designed to reduce loads 
induced by extreme (design level) turbulence must also oper¬ 
ate satisfactorily in continuous turbulence. Therefore, the ef¬ 
fects of any nonlinear systems must be evaluated for appropri¬ 
ate levels of continuous turbulence. SDG methodology for 
continuous turbulence, used with the search technique de¬ 
scribed in Reference 4, is a practical means of performing this 
evaluation. 

5. CONCLUSIONS 

While existing gust load criteria have yielded satisfactory lev¬ 
els of flight safety, it can be seen from the data presented 
above that SDG extreme turbulence methodology produces 
levels of dynamic response that are not provided by either the 
variable development length (1-cos) discrete gust or continu¬ 
ous turbulence methods. The degree of variation between 
loads predicted by the three methods is a function of the re¬ 
sponse characteristics of the individual aircraft. The results 


presented above may represent an extreme for commercial 
aircraft because of the relatively wide separation between rigid 
body and elastic mode frequencies. SDG modeling accommo¬ 
dates this frequency separation and yields equally probable 
loads for all components of the aircraft. Results from analyses 
of other aircraft may not show as much variation between 
methods. However, SDG methodology is the only one of the 
three analysis methods that produces consistent results for all 
aircraft. This is the result of maximizing both rigid body and 
elastic responses due to a statistical model of extreme turbu¬ 
lence, regardless of the frequency response characteristics of a 
particular aircraft. 

A design criterion based upon the SDG model of extreme tur¬ 
bulence has several advantages over the discrete gust and con¬ 
tinuous turbulence criteria presently in use. It uses a statisti¬ 
cal model of the type of turbulence that is likely to produce 
structural overload and failure. It is independent of the air¬ 
craft configuration. It provides equally probable limit loads, 
which provide the potential for improved structural efficiency 
of the aircraft. It would reduce the engineering effort required 
for certification by eliminating the need for multiple require¬ 
ments for design limit loads due to turbulence. 

The next step in the evolution of criteria for design limit loads 
due to turbulence should be the acceptance of SDG methodol¬ 
ogy for extreme turbulence. The process of acceptance should 
include a review of the design gust intensities and the prob¬ 
ability factors that maintain a constant probability of exceed¬ 
ance for critical gust profiles consisting of more than one 
ramp gust. 
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APPENDIX 

Statistical Discrete Gust Methodology for Design Limit 
Loads 

The SDG model of extreme (design level) turbulence repre¬ 
sents equally probable critical gust profiles with one (1-cos) 
ramp gust or two such gusts of alternating sign. The maxi¬ 
mum incremental velocity for a single ramp is a function of 
the gust intensity factor, U Q , and its development length, H, 
and is given by: 

U = U 0 ’H I/6 Ft/Scc (TAS) (A-l) 


For the present, a definition for the gust intensity factor based 
on the reference gust velocity of JAR 25.341, Amendment 
91/1 is suggested. 

Uo = U r ef/350' /6 /Va, Ft/Sec(TAS)/Ft l/6 (A-2) 


The velocity profile for each ramp gust, is given as a function 
of gust penetration distance, s, the maximum velocity, U, and 
the gust development length, H. 


.. U( TtsA 

u(s) = •— 1 - cos- 

2 V H ) 

u(s) = U 


for 0 < s < H (A-3) 
for s > H (A-4) 


To maintain a constant probability of exceedance for the more 
complex two-ramp critical profiles, a probability factor, P-,. is 
applied to the incremental velocities. A provisional value of 
0.75 is currently being used for P 2 . 

For the general case, worst case loads for extreme turbulence 
arc determined by varying the development length for a single 
ramp profile and the development lengths and spacing of ramp 
gusts of alternating sign in a two-ramp gust profile. 
Superposition of results from a family of single ramp gust is 
used to determine worst case loads and responses and the as¬ 
sociated critical gust profiles for linear systems. Development 
lengths for the family of ramp gusts must satisfy tw'o require¬ 
ments. First, the range of development lengths must be broad 
enough to assure maximum excitation of both rigid body and 
elastic modes of the aircraft. Second, the spacing of the de¬ 
velopment lengths must be sufficiently close to assure Finding 
the required extreme values. The following guidelines are 
suggested for use in initial analyses. However, the response 
characteristics of a particular aircraft will dictate the required 
range and spacing of the development lengths defining the 
family of gusts. 


Minimum development length, 15 feet 

Maximum development length, 

vertical gusts 400 feet 

lateral gusts 1,200 feet 

Incremental development length, 12 % 


Response analyses for each member of the family of ramp 
gusts will yield maximum and minimum values of opposite 
sign for desired loads and responses. Since the most complex 
SDG extreme turbulence profile consists of two ramp gusts of 
alternating sign; only the overall maximum and minimum ex¬ 
treme values (E| and E 2 ) are required to define worst case 
loads and responses. The design amplitude, y , for a given 


load or response is the maximum of the peak values, y k re¬ 
sulting from one and two ramp critical profiles. 

k 

Yk= p k£N for k= 1,2 (A-5) 

i=l 

P| and P 2 have values of 1.0 and 0.75, respectively. The de¬ 
sign amplitude is then: 

y =max(y,,y 2 ) (A-6) 

Associated with the extreme values, E, and E 2 , values arc gust 
development lengths, Hj and H-,, and times of occurrence, T, 
and T-, . These times and development lengths are used to de¬ 
fine critical profiles and (if desired) determined if any over¬ 
lapping of ramps in a two-ramp profile has occurred. Critical 
two-ramp profiles are defined by superimposing the ramp 
gusts with alternating signs such that the extreme values will 
add to form the worst case value for the load or response. 
The critical profile begins with the ramp, H 2 , causing the sec¬ 
ond extreme value, E-,. followed at the appropriate time by the 
ramp, H p (with opposite sign) that caused the first extreme 
value, E|. The time to start the second ramp, T , is given by: 

T s = T 2 -Ti (A-7) 

The possibility of overlapping ramps is checked by comparing 
the times at which the initial ramp ends, Tj, and the time at 
which the final ramp begins, Tf. For non-overlapping ramps: 

T f - T, > 0 (A-8) 

Experience has shown that overlap does not occur except 
W'hen one of the extreme values is the result of inertia loads 
caused by gust loading on other parts of the aircraft. In the 
case where local extreme values are the result of local gust 
loads, overlap has not been observed. However, one case of 
apparent overlap has been encountered. For a configuration 
with a highly swept wing and low wing loading the first ex¬ 
treme values of shear and bending moment near the wing tip 
were due to inertia loads. These inertia loads resulted from 
the inboard portion of the wing penetrating the gust and gen¬ 
erating significant rigid body acceleration before the outboard 
wing encountered significant gust loading. The result is iner¬ 
tia loads of larger magnitude than the loads associated with 
dynamic overshoot of loads due to local gust loading. This 
condition has not been encountered for a critical load case. 
However, if such loads were present for critical conditions, 
they would be used in the calculation of y , because this is a 
real loading condition. The overlap of gust ramps resulting 
from this kind of loading is not considered to be a violation of 
the SDG model because the initial extreme value is not the 
result of local gust loading. 
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A REVIEW OF GUST LOAD CALCULATION METHODS AT de HAVILLAND 

John Glaser 

de Havilland Inc., 123 Garratt Blvd., Downsview 
Ontario, CANADA, M3K 1Y5 


SUMMARY 

The development of an analysis system for the routine 
calculation of gust response loads is reviewed in some 
detail in this report. While the system provides adequate 
design strength margins, a more robust and effective 
system would reduce user workload, computer costs and 
analysis time. It is suggested that other analysis systems 
could benefit similarly, particularly when considering 
the demands imposed by highly nonlinear aircraft sys¬ 
tems, the trend toward full finite element structural 
dynamic models and the relentless quest for structural 
efficiency. It is proposed that improvements in analysis 
systems could evolve from the collective experience in 
gust loads methodologies acquired within the aeronauti¬ 
cal community. To capitalize on that collective experi¬ 
ence, it is recommended that a working group of gust 
load specialists be formed to assemble and evaluate cur¬ 
rent and promising methods for calculating gust loads 
and to recommend standardized airplane test cases, both 
rigid and elastic, for validating analysis methods and 
results. 


1 INTRODUCTION 

This paper is not about the modeling of atmospheric tur¬ 
bulence or the modeling of aircraft structures; these 
matters, although extremely important, have been and 
will continue to be investigated by such bodies as the 
AGARD Structures and Materials Panel (SMP), the 
FAA-hosted Gust Specialists’ Meetings, the AIAA/ 
ASME/ASCE/AHS/ASC Structures, Structural Dynam¬ 
ics and Materials Conferences, and the International 
Modal Analysis Conferences (IMAC) to mention a few. 

Rather, the purpose of this paper is to focus attention on 
the methods used for calculating response loads given 
that models for atmospheric turbulence and the airplane 
structure have been defined. Loads methods for highly 
nonlinear aircraft in continuous turbulence fall within 
this category but will not be discussed here; this prob¬ 
lem is currently receiving considerable attention by the 
FAA-hosted Gust Specialists’ meetings and other 
groups. 


Figure 1 provides a flow diagram for the method used at 
de Havilland. We begin by assuming linearity and set up 
the flight dynamic and structural dynamic equations of 
motion for the airplane in terms of its rigid and elastic 
modes of vibration. Frequency response functions 
(FRFs) for modal coordinates and hence for loads are 
calculated assuming flight through one dimensional sine 
wave gust fields. We cam then continue in the frequency 
domain to calculate loads for continuous turbulence 
(usually von Karman) or we can transform into the time 
domain to obtain loads for a step gust from which 
versed sine or (1-cos) gust responses can be calculated 
using the convolution integral. 

As we look more deeply into the process we find that a 
number of analytical assumptions and computational 
procedures me required to calculate gust response loads, 
each having their influence on user workload, computer 
requirements, analysis time and the results obtained. 
This paper will describe some of the approaches taken 
by de Havilland in developing its gust loads analysis 
system. This review is presented with the hope that it 
will stimulate renewed attention within the aeronautical 
community on this important aspect of the gust loads 
analysis problem. 

The basis of experience for this report is de Havilland’s 
family of turboprop commuter airplanes with maximum 
take-off weights ranging from 30,000 to 45,000 pounds 
(13,500 to 20,500 kg), with maximum design cruise 
speeds of about 250 kts eas and with maximum design 
altitudes of 25,000 feet (7,500 m). The analysis models 
for these aircraft were developed using stick models for 
the structure and unsteady strip theory aerodynamics 
with Yates corrections (Ref. 1). 


2 SETTING UP THE EQUATIONS OF MOTION 
- SIMULATING FLIGHT DYNAMICS 

Limited accountability was the first method used at de 
Havilland (some 30 years ago) to include elastic 
responses due to gust encounter. The gust model used 
was discrete with versed sine shape. The general 
approach was to compare elastic-to-rigid gust response 
ratios for the new design with those of an acceptable ref- 


Paper presented at a Workshop on “Aircraft Loads due to Turbulence and their Impact on 
Design and Certification” held in Lillehammer, Norway on 5 May 1994. 
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erence airplane and to augment the prevailing gust crite¬ 
ria loads by the increase in dynamic factor, if present. 
The aircraft model used included only rigid body heave 
and first wing bending modes; there was no attempt to 
include flight dynamics in the gust calculations in those 
days. 

Flight dynamics first became important at de Havilland 
with the introduction of the continuous turbulence crite¬ 
rion of FAR P;u1 25, Appendix G, because of the rela¬ 
tively high energy at low frequencies that characterize 
continuous turbulence models. Our first approach to 
simulate combined flight dynamic and structural 
dynamic effects was to write the equations of motion in 
a body fixed coordinate system. However, because we 
ran into some analytical difficulties which we felt could 
not be resolved in the time available, we quickly 
retreated to the customary space fixed coordinate sys¬ 
tem . 

Good flight stability characteristics (phugoid, short 
period, dutch roll) were achieved in the space fixed sys¬ 
tem by manually adjusting some of the aerodynamic 
terms in the flight dynamic matrix. This method was 
subsequently automated by introducing a so-called “sky 
hook” - a single aerodynamic panel - located at the air¬ 
craft eg to achieve the flight stability characteristics 
required. Consistent with the manual adjustment 
method, the sky hook was arranged to affect only the 
flight dynamic equations (ie, panel forces did no work in 
elastic modes). While the sky hook method was simple 
to implement, we felt uneasy with the approach because 
of the inconsistency introduced between rigid and elas¬ 
tic degrees of freedom. 

Our current approach for setting up the combined equa¬ 
tions is to distribute a set of aerodynamic panels over 
the length of the fuselage that affect rigid and elastic 
modes consistently and that provide the flight stability 
characteristics required. While clearly an improvement 
over the sky hook method, we are still not satisfied 
because the aerodynamic coefficient distribution 
required for flight stability does not agree well with 
fuselage aerodynamic test data. 

Some vertical gust loads results for continuous turbu¬ 
lence are presented in Table 1 for three flight dynamic 
simulation methods; no attempt to simulate stability, the 
sky hook method and the fuselage panel method. The 
results show that simulating flight stability is important 
and that it matters how that simulation is achieved. 


3 GUST LOAD RESPONSE CALCULATIONS 

The Mission Analysis Method and N 0 

As part of our evaluation of the Mission Analysis 
method of Appendix G. we ran into some difficulty in 
establishing tut acceptable method for defining N t) , the 
number of zero crossings in unit time. This parameter is 
calculated from a definite integral over frequency with 
convergence depending on the specific formulation of 
the Sears function used. With guidance from Refs. 2 anil 
3, we investigated the acceleration response for the case 
of a rigid airplane free in heave only. We found that con¬ 
verged values for N„ could be obtained using approxi- 
mat ions to the Sears function that decayed as k"“ or 
faster, for large values of k, where k is the reduced fre¬ 
quency. For example, the two-term Kussner approxima¬ 
tion for the incompressible Sears function (Ref. 4). gave 
a value for N 0 of 2.3 whereas the Houbolt compressible 
formulation (Ref. 3) gave 2.6. Reference 2 suggests that 
physically more realistic values of N u are obtained if the 
integration is truncated at a frequency corresponding to 
95% of the converged value of A, the response magni¬ 
tude due to unit root-mean-square (rms) gust intensity. 
Using this definition, the values obtained for N 0 were 
0.85 and 0.81, respectively. Due to the ambiguity in 
defining N 0 and for more compelling philosophical rea¬ 
sons regarding the use of Mission Analysis for defining 
design loads, we decided very early to adopt the Design 
Envelope method of Appendix G. 

The foregoing investigation on N u caused us to wonder 
what effect the various Sears function approximations 
(Fig. 2) might have on loads. This is shown in Table 2 
for an elastic airplane. The results have been normalized 
with respect to Fung’s formulation which gives a very 
close approximation to the exact, incompressible Se;us 
function (Ref. 5). Compared to Fung, the results are 
noted to vary from +37 to -6 7 . Houbolt’s compressible 
formulation (Ref. 3) for M=() is also observed to he in 
good agreement with Fung. Time did not permit our 
investigation into lifting surface formulations for Sears 
functions. 

The Force Summation Method and Other Methods 
for Calculating Loads 

As mentioned earlier, our first method for taking elastic 
effects into account required the calculation of elastic to 
rigid gust response ratios. This required the use of a 
Force Summation Method (Refs. 6, 7, and 8) as it is the 
only method that can be used to calculate response loads 
for a rigid aircraft. For a variety of reasons. Force Sum- 
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mation remains the preferred method for calculating 
loads at de Havilland as well as at most airplane manu¬ 
facturers (based on a survey conducted by the AGARD 
SMP). Furthennore, simplified rigid airplane analysis 
models can be used to provide important checks on 
results and it is always useful to keep a “calibrated” eye 
on elastic-to-rigid response ratios as a reality check. 

Another very attractive reason for using the Force Sum¬ 
mation method is that fewer modes need to be carried to 
obtain converged dynamic load results than for methods 
such as Mode Acceleration or Mode Displacement (Ref. 
6). This is illustrated in Fig. 3 for wing root and horizon¬ 
tal tail root shear, bending moment and torque loads. 
Note that in some cases, the Mode Displacement 
method has not converged fully even with 80 elastic 
modes and frequencies up to approximately 300 Hz. 
Furthennore, it is difficult to know whether convergence 
has actually been achieved when using that method. 


4 DISCRETE GUST LOADS AND BALANCED 
LOADS 

Transforming The Frequency Response Function 

For calculating discrete gust loads, the procedure used at 
de Havilland is to first calculate the response for a unit 
step gust from a Fourier transformation of the frequency 
response function and to use the convolution integral to 
obtain the load response for any arbitrary gust pattern of 
interest (Fig. 1). 

It can be shown (Ref. 9) that the response, A(t), to a unit 
step gust can be calculated from the following transfor¬ 
mation expression: 

2 f °° I (co) coswt 

A (t) = R (0) + - ——-dco (1 

« J 0 ® 

where H (to) = R(co) +/I(co) is the frequency re¬ 
sponse function for the particular load of interest, co is 
the circular frequency, t is time and i = J—\. 

Due to the rigid body modes, H(co) is not well defined at 
co = 0. This problem is addressed (Ref. 10) by setting 
R(0) = R(e) and I (0) = 0, and by assuming linear 
variations in their values over the frequency range 0 to 
e. Here, e is a small value of co, and is determined by 
computational experimentation so that expected behav¬ 
iours for A(t) are obtained at large values of time, t. For 
example, loads should converge to zero and velocities to 
constant values. 


The application of the convolution integral for calculat¬ 
ing the response to an arbitrary discrete gust profile is 
straightforward and requires no comment. However, it 
should be noted that some gust patterns, such as the 
familiar versed sine, can mask inaccuracies that may 
exist in the step response function at t = 0 and t —»°° 
before convolution is applied. This process for calculat¬ 
ing transient loads is subjective and tedious; a more 
robust methodology would improve effectiveness. 

Balanced Loads 

The need for balanced loads became important when the 
Stress Group at de Havilland requested time correlated, 
externally applied load distributions for their calcula¬ 
tions. In this context, balanced loads means that the 
externally applied load distribution, including inertia 
reactions, produces no residual force or moment when 
integrated over the whole airplane. 

This requirement exposed several new problems. 

(1) We first discovered an inconsistency in our applica¬ 
tion of the force summation method for calculating 
frequency response functions which, once realized, 
was easily resolved. 

(2) The sky hook method for achieving required flight 
stability characteristics (already discounted for 
other reasons as explained above) produced dis¬ 
crete and sometimes significant loads at the aircraft 
eg that were difficult to rationalize when calculat¬ 
ing internal loads. This reinforced our resolve to 
establish realistic aerodynamic distributions for 
both flight stability and loads. 

(3) The Stress Group’s method for calculating internal 
loads for the fuselage expected shear distributions 
only whereas the Loads Group provided both local 
shear and bending moment. This was a process 
problem that was easily resolved once realized. 

In summary, the degree of balance achieved for an air¬ 
craft in free flight provides a useful test on the loads 
methods used. 


5 CONTINUOUS TURBULENCE, BALANCED 
LOADS AND THE SDG METHOD 

When it came time to calculate Design Envelope loads, 
we were confronted with a new problem: how to define 
balanced load sets for flight in random turbulence for 
which there is no time basis for correlation. 
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By way of a solution, we undertook an investigation of 
J. G. Jones’ Statistical Discrete Gust (SDG) method 
(Ref. 11). This method gave loads approximating those 
obtained for the Design Envelope method, once the 
design gust intensity was calibrated (Chapter VI. Ref. 
12). Being a time domain method, time correlated find 
hence balanced load distributions could be defined to 
provide a solution to our problem. The SDG method had 
other attractions as well; for example, it could be used to 
represent other (more realistic) models of atmospheric 
turbulence (discrete or continuous) and, being a time- 
domain method, it could be used for calculating Design 
Envelope loads for aircraft with highly nonlinear sys¬ 
tems. 

In the end, however, the SDG method was not used at de 
Havilland for design because a parallel study based on 
Noback's work (Ref. 13) showed that balanced load dis¬ 
tributions for random turbulence could in fact be defined 
using correlation coefficients. Therefore, the effort to 
validate the SDG method for design and certification 
was avoided. With reference to Fig. 3, it should be noted 
that balanced load conditions are always obtained irre¬ 
spective of the number of modes carried. 

The SDG method, as proposed at the time, was gener¬ 
ally concluded to be overly complex for constructing 
critical gust patterns. While a simpler formulation has 
since been proposed (not formally documented to date), 
due to limited resources, we have not had an opportunity 
to evaluate this new proposal. Although we believe that 
the SDG method and other similar methods merit fur¬ 
ther evaluation, we wish to express our concern about 
the general level of complexity that these new methods 
require. Complex methods are costly to implement and 
validate and their standardized application in industry to 
achieve consistent results will be virtually impossible 
unless standtirdized algorithms are made generally 
available. 

Currently, a NASA-developed code for determining 
critical discrete gust patterns for von Karman continu¬ 
ous turbulence has been made available for evaluation 
(Ref. 14). The method is based on Matched Filter The¬ 
ory ruid is applicable to airplanes with linear and nonlin¬ 
ear systems. The effectiveness of this development and 
distribution approach should be monitored because it 
could prove to be an effective means for the standardiza¬ 
tion of complex computational methods within the aero¬ 
nautical community. 


6 AERODYNAMIC CORRECTION FACTORS 

Aerodynamic corrections are used to modify strip or 
doublet lattice aerodynamic formulations to account for 
the effects of three-dimensionality, viscosity, or tran¬ 
sonic flow. Generally, these corrections derive from 
wind tunnel data or more sophisticated aerodynamic 
codes. The correction method used at de Havilland 
applies factors to the derivatives and cross-derivatives 
for each streamwise strip. This approach was recently 
applied in showing compliance with the new require¬ 
ment for control surface loads, JAR 25.391, introduced 
by JAR Amendment 91/1. 

Two other aerodynamic correction methods ;ire summa¬ 
rized by T. Noll (Ref. 15) viz: applying correction fac¬ 
tors to airfoil section properties and applying correction 
factors to the pressures or downwashes on individual 
boxes. Due to limited resources, these alternate aerody¬ 
namic correction methods have not been evaluated at de 
Havilkmd to date. 


7 DISCUSSION 

This report has reviewed the development of an analysis 
system for calculating gust loads and has indicated some 
of the difficulties experienced in developing that system. 
Due to high development costs, the preferred approach 
of evaluating all methods and choosing the best could 
not be pursued: rather, we made what appeared to be 
rational decisions on methodology at the outset and built 
our system accordingly. The end result, although accept¬ 
able for design safety, requires considerable diligence 
on the part of the user with the obvious consequent 
impact on prognun costs and schedule. 

While our analysis system would benefit from revisions 
in the methods used, we believe other systems could 
benefit similarly, particularly in view of some emerging 
requirements and trends. Consider the following: 

• The design of aircraft with highly nonlinear systems 
(active controls) which will require more complex 
and computer costly analysis routines, 

• The trend toward full finite element structural dy¬ 
namic models to achieve more efficient structural 
designs and design processes. 

• The ever increasing need to reduce design cost and 
cycle time while developing more efficient structures. 
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In summary, we perceive a two-fold need; the need for 
more effective gust load analysis systems to design effi¬ 
cient structures within tight schedules, and the need to 
validate the results of the analysis systems used. 


8 RECOMMENDATIONS 

Especially over the last 10 years, considerable experi¬ 
ence has been acquired in the development and applica¬ 
tion of routine gust load calculation methods. While this 
experience resides mostly with the end users in industry, 
considerable contribution find insight have emanated 
from Research Institutes as well. It is therefore evident 
that, if further gains tire to be made in gust loads analy¬ 
sis methods, a united participation of the experts in the 
field will be required. 

It is therefore recommended that a working group of 
gust loads experts from Industry and Research Institutes 
be formed and that this working group take into consid¬ 
eration the following recommendations as a basis for its 
terms of reference. Before proceeding further, it should 
be emphasized that it is not the intention of this proposal 
that a single, standardized analysis system evolve from 
the deliberations of such a working group but that more 
effective analysis systems tailored to specific needs be 
defined and validated from a digest of the collective 
experience accumulated by the experts in the field. 

(1) As a first step toward achieving more effective 
analysis systems for routine gust load computa¬ 
tions, it is recommended that a compendium of cur¬ 
rently used and promising methods be assembled 
that would provide a detailed definition of each 
method giving its advantages, disadvantages, limi¬ 
tations and examples of results. This compendium 
should consider the following very different 
approaches: 

• The method for setting up and solving the equa¬ 
tions of motion for an elastic airplane in body- 
fixed coordinates as presented in Ref. 16 by 
Waszak and Schmidt of Purdue University. 

• The procedure presented by Lusebrink and Son- 
der of Deutsche Airbus (Chapter V, Ref. 10) 
which sets up and solves the equations of 
motion in the time domain using indicia! 
unsteady aerodynamic formulations. 

• The procedure presented in Ref. 17 by Anthony 
Pototzky and Boyd Perry III of NASA Langley 


which sets up and solves s-plane and state space 
formulations of the equations of motion. 

• The procedure presented in Ref. 18 by Ph. Nicot 
and C. Petiau of Dassault Aviation which sets 
up and solves the equations of motion using a 
loads basis method as opposed to the more com¬ 
mon modes basis method. 

(2) It is recommended that standard test cases be 
defined to provide a basis for validating newly 
developed analysis methodologies or to verify 
problem setups using established analysis systems. 
One standard case that has proven invaluable is the 
so-called Pratt gust formula (FAR 25.341) that 
gives the discrete gust load factor for a rigid air¬ 
craft with freedom in heave only. Additional stan¬ 
dard cases that include other rigid body freedoms 
and structural flexibility should also be considered 
in order to achieve more comprehensive valida¬ 
tions. 
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TABLE 1 EFFECT OF FLIGHT DYNAMIC SIMULATION ON LOADS 

CONTINUOUS TURBULENCE VERTICAL GUST RESULTS ELASTIC AIRPLANE 



FLIGHT DYNAMIC SIMULATION METHOD 

NO SIMULATION 
WING + TAIL 
AERO ONLY 

SKY HOOK 

FUSELAGE 

PANELS 


WING 

VERTICAL SHEAR 

1.00 

.94 

ROOT 

VERTICAL BM 

1.00 

.95 

RESPONSE 

TORQUE 

1.00 

.99 

HORIZONTAL 

VERTICAL SHEAR 

1.00 

1.10 

TAIL ROOT 
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.87 
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STABILITY 

DAMPING RATIO 

.33 

.50 



TABLE 2 EFFECT OF SEARS FUNCTION APPROXIMATIONS ON LOADS 
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FIG. 1 GUST RESPONSE CALCULATION FLOW DIAGRAM 
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SPECIAL EFFECTS OF GUST LOADS 
ON MILITARY AIRCRAFT 

by 

Dr. John C. Houbolt 
Chief Aeronautical Scientist (retired) 
NASA Langley Research Center 
Hampton Va. 23681 


ABSTRACT 

In the operation of airplanes, atmospheric turbulence.creates 
a broad spectrum of problems. The nature of these problems is 
presented in this paper. Those that are common to both the 
commercial carriers and to the military fleet are discussed 
first. Attention is then focused on the problems that are of 
special concern in military operations. An aim is to bring out 
the need for continued effort in the gust research area. 


GENERAL 

The flight of aircraft through rough air has been studied for 
many years, but is still of major importance. Upon examination, 
one finds the breadth of the problems created somewhat 
surprising. In general terms, severe turbulence has been of 
concern because of the control "upset" problem and because of its 
influence on the static design strength of the airplane; moderate 
turbulence is of concern because it leads to problems in 
achieving precise flight and is a major source of fatigue damage 
to the structure. Many reviews have been written to bring out the 
specific problems created by turbulence. The points made are 
valid even to this day. A look at some of these points serves 
well to illuminate the extent and nature of the flight problems 
that are created by atmospheric turbulence. Back in 1964 alone, 
there were 8,500 pilot reports of turbulence encounter of the 
unanticipated clear air variety. Of these 2,060 (about 25%) were 
of severe nature, 135 were extreme. 

Generally, turbulence encounter affects all aviation 
activities, including defense and civil operation, in the 
following ways: 

a) Mission cancellation and diversion around turbulence 


Paper presented at a Workshop on “Aircraft Loads due to Turbulence and their Impact on 
Design and Certification” held in Lillehammer, Norway on 5 May 1994. 
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b) Increased aircraft inspection 

c) Temporary loss of aircraft control; resultant injury to 
passengers and crew requiring compensation or medical 
treatment 

d) Shortening of aircraft life time due to metal fatigue 

e) Need for special aircrew training 

f) Shortage of airspace due to increased number of "baby 
jets" 

Certain aspects unique to civil and military operation are 
also present. In 1964, for civil aircraft alone, it is estimated 
that due to injury, diversions, training, information and 
inspection, a total expense of $18,000,000 was involved due to 
turbulence encounter. This figure does not include loss of the 
use of grounded aircraft, loss of time due to injuries, or 
expense of settling claims. Also the increase in number of 
general aviation jets, operating in the 20-40,000 foot 
environment, where restrictions vertically and horizontally may 
be imposed at times due to turbulence, result in increased 
economic concern to operators. 

Military aircraft have experienced numerous accidents 
and incidents because of turbulence. Additional difficulties 
presented to military operations include: 

a) Hazards in air-to-air refueling 

b) Reduced available airspace for conducting simulated 
missions 

c) Degradation of reconnaissance sensor and bombing 
accuracies 

d) (More as discussed in a later section) 

The cost of turbulence to DOD operations is difficult to 
evaluate, but estimates for the 3-year 1963-1965 period due to 
aircraft lost and damaged are placed at around $30,000,000. Costs 
of special inspections, minor repairs, minor injuries, and for 
special training are not known. In addition, value cannot be 
placed on loss of pilots or on the effect on actual combat 
missions due to lack of advanced information on the location of 
clear air turbulence. 
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GUST RESEARCH EFFORT 

We list here in brief fashion the nature of the research 
effort that has been underway for years to understand and solve 
the problems that arise from gust encounter. It is convenient to 
list the studies in four broad categories. 

Turbulence measurement: This effort has involved three main 
sources. Low altitude measurement on tall towers, specific 
research flights on well instrumented airplanes to measure the 
turbulence, and through means of special recorders in aircraft to 
measure turbulence encounter and severity during routine 
scheduled flights. 

Turbulence modeling; By analytical studies, and through use of 
the turbulence measurements, mathematical modeling of turbulence 
has been made. The modeling has been in terms of at least four 
related aspects: the concept of discrete or isolated type gust 
encounter, the concept of continuous (random) turbulence 
encounter form, the establishment of turbulence design severity 
levels, and the determination of the time spent in turbulence. 

Airplane response analysis due to turbulence: Here studies 
are made to determine how the airplane responds to turbulence. 
Modeling consists of representing the airplane by various 
structural degrees of freedom, and of using non-steady 
aerodynamics to account for the airplane motion and gust 
penetration effects. The results sought are the maximum loads, 
acceleration, stresses, and load encounter curves for airplane 
fatigue analysis and life prediction. In turn these results are 
used to design the structural members, particularly the wing, and 
to design special devices such as gust loads alleviation systems. 

Specifications and regulations: The research studies do not 
directly write the specifications or regulations. However the 
studies are fed as guidelines to the regulators so that they can 
use the findings to establish the design specifications and 
procedures for turbulence encounter that ultimately appear in the 
regulations. 


ADDED SPECIAL MILITARY PROBLEMS 

Some of the turbulence problems that are special for military 
aircraft have been mentioned in the preceding sections. By way of 
emphasis we list here certain problems of turbulence encounter 
that are of particular concern, some of which are unique to 
military type aircraft. 
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Turbulence encounter has a marked effect on degrading the 
structural life of the aircraft (this has become to be known as 
the aging aircraft problem). In fact, gust encounter accounts for 
40-50% of the damage hours of the structure. 

In low level missions, which is of particular significance to 
military operations, but where turbulence can be rather severe 
and prevalent, damage hours build up ten times more rapidly 
than in higher altitude flight. 

T-tailed aircraft take a rather rough beating during air¬ 
refueling operations due to turbulence generated by the tanker 
aircraft. 

Gust alleviation systems are often needed to control airplane 
motion, such as yaw dampers. These devices are also needed to 
improve the ride comfort for the pilot, and to improve sensor 
capability and bombing accuracy. 

We note that the items listed in this section are in need of 
greater research effort. 

We also note that, whereas turbulence causes the most damage 
to the wings, the ground to air cycle (fuselage pressurization) 
creates most of the damage hours to the fuselage. 


CLOSURE 

For possible use to obtain a quick overview of this report, 
or for presentation purposes, the contents of this report are 
summarized in three condensed listings in the three pages that 
follow. 
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ADDED SPECIAL MILITARY PROBLEMS 


MARKED EFFECT ON DEGRADING STRUCTURAL LIFE 

OF THE AIRCRAFT (THE AGING AIRCRAFT PROBLEM) 


GUSTS ACCOUNT FOR 40-50% OF THE DAMAGE HOURS 
OF THE STRUCTURE 


IN LOW LEVEL MISSIONS DAMAGE HOURS BUILD UP 
10 X MORE RAPIDLY 


T-TAILED AIRCRAFT TAKE A ROUGH BEATING DURING 
AIR-REFUELING (TURBULENCE FROM THE TANKER 
AND THE ATMOSPHERE) 


GUST ALLEVIATION IS NEEDED TO CONTROL AIRPLANE 

MOTION AND STRUCTURAL MODES, FOR RIDE COMFORT, 
AND FOR SENSOR AND BOMBING ACCURACY 


NOTE: GROUND TO AIR CYCLE (FUSELAGE PRESSURIZATION) 
CREATES MOST OF THE DAMAGE HOURS TO THE 
FUSELAGE 
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